Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,442,078 A→G 3.9% K499K (AAA→AAG thiC → thiamin (pyrimidine moiety) biosynthesis protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077793,442,0780A→G 3.9% 13.1 ‑0.3 ‑0.0206K499K (AAA→AAGthiCthiamin (pyrimidine moiety) biosynthesis protein
Reads supporting (aligned to +/- strand):  new base (4/4):  ref base (87/111):  total (91/115)
Fisher's exact test strand distribution p-value = 7.34e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.45e-01

Only 100 of 250 total aligned reads displayed.

CTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGCGAAAGGGCATCCGGGCG  >  NC_007779/3441978‑3442178
                                                                                                     |                                                                                                    
cTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCaaaa                                                                                                      <  1:9434003/101‑1 (MQ=255)
 tACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAag                                                                                                     <  1:7395214/101‑1 (MQ=255)
 tACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAag                                                                                                     <  1:7199490/101‑1 (MQ=255)
 tACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAag                                                                                                     <  1:10803723/101‑1 (MQ=255)
  aCCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAaga                                                                                                    >  1:1537459/1‑101 (MQ=255)
   ccGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAagag                                                                                                   <  1:7241198/101‑1 (MQ=255)
   ccGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAagag                                                                                                   <  1:14204352/101‑1 (MQ=255)
     gacatcgcgcCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTGCGCAGGCTGGTACCGCCATGCTCTGCTACGTCACGCCGAAGGAGCa                                                                                                 >  1:12888273/7‑101 (MQ=255)
       tatTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATc                                                                                               >  1:2825191/1‑101 (MQ=255)
        atcgcgcCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTGCGCAGGCTGGTACCGCCATGCTCTGCTACGTCACGCCGAAGGAGCACCt                                                                                              <  1:5969454/98‑1 (MQ=255)
        atTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCt                                                                                              <  1:10447992/101‑1 (MQ=255)
        atTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCt                                                                                              <  1:8130984/101‑1 (MQ=255)
         ttGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTg                                                                                             <  1:13141453/101‑1 (MQ=255)
          tGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCTGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTgg                                                                                            <  1:3691954/101‑1 (MQ=255)
          cgcgcCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTGCGCAGGCTGGTACCGCCATGCTCTGCTACGTCACGCCGAAGGAGCACCTgg                                                                                            >  1:13294079/2‑101 (MQ=255)
           gcgcCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTGCGCAGGCTGGTACCGCCATGCTCTGCTACGTCACGCCGAAGGAGCTCCTggg                                                                                           <  1:1756912/101‑1 (MQ=255)
             gcCGGGCTATGACCACTTCACGTCGGGGATTGGAGCGGCGTTGTTTGGCTGGTTTGGCTGCGCGACGCTCTGTTTCGTAACGCCAAAAGAGCATCTgggcc                                                                                         >  1:1385426/1‑99 (MQ=255)
               cGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGGTGCTCTGTTACGCAACGACAAAAGAGCATCTgggagtg                                                                                       >  1:9298916/1‑97 (MQ=255)
               cGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTg                                                                                       >  1:8530109/1‑101 (MQ=255)
                gggCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGc                                                                                      <  1:11083446/101‑1 (MQ=255)
                gggCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGc                                                                                      <  1:2231772/101‑1 (MQ=255)
                  gCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGccc                                                                                    >  1:12103995/1‑101 (MQ=255)
                   cTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCa                                                                                   <  1:13576131/101‑1 (MQ=255)
                    tATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCaa                                                                                  >  1:12167644/1‑101 (MQ=255)
                     aTGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAAt                                                                                 >  1:12506214/1‑101 (MQ=255)
                      tGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATa                                                                                <  1:1049386/101‑1 (MQ=255)
                      tGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATa                                                                                <  1:9679065/101‑1 (MQ=255)
                      tGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCCCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATa                                                                                <  1:8840914/101‑1 (MQ=255)
                      tGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCCCGATGCCCTGTTACGTAACGCCAAAAAAGAAACGGGGTCTGCCCAATa                                                                                >  1:6347068/1‑101 (MQ=255)
                       gACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGGTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATaa                                                                               >  1:2204692/1‑101 (MQ=255)
                        aCCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATaaa                                                                              <  1:8186025/101‑1 (MQ=255)
                         ccACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCAAAAAGAGCATCTGGGTCTGCCCAATAaag                                                                             <  1:12191374/101‑1 (MQ=255)
                          cACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAaaga                                                                            <  1:53333/101‑1 (MQ=255)
                             ttCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAagaaga                                                                         <  1:8591628/101‑1 (MQ=255)
                               cACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATg                                                                       >  1:3347185/1‑101 (MQ=255)
                                aCGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGt                                                                      <  1:5627652/101‑1 (MQ=255)
                                aCGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGt                                                                      <  1:4174154/101‑1 (MQ=255)
                                aCGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGt                                                                      >  1:7637627/1‑101 (MQ=255)
                                 cGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGtt                                                                     >  1:6715266/1‑101 (MQ=255)
                                   tCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTaa                                                                   <  1:5795671/101‑1 (MQ=255)
                                    cGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAg                                                                  <  1:8931611/101‑1 (MQ=255)
                                     ggggATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGc                                                                 <  1:9592434/101‑1 (MQ=255)
                                      gggATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCa                                                                <  1:2871870/101‑1 (MQ=255)
                                        gATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTGCGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAgg                                                              <  1:4920775/101‑1 (MQ=255)
                                         aTTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAggg                                                             <  1:13659244/101‑1 (MQ=255)
                                          ttGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAgggg                                                            <  1:10026019/101‑1 (MQ=255)
                                          ttGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAgggg                                                            <  1:7663755/101‑1 (MQ=255)
                                          ttGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAgggg                                                            >  1:3135513/1‑101 (MQ=255)
                                           tGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGc                                                           >  1:1585312/1‑101 (MQ=255)
                                             gTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCtt                                                         <  1:14284579/101‑1 (MQ=255)
                                              tGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTa                                                        <  1:11839611/101‑1 (MQ=255)
                                              tGCGGCGATGATTGGCTGGTTTGGCTGCGCGACGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTa                                                        >  1:1704315/1‑101 (MQ=255)
                                               gcggcgATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTAt                                                       >  1:12668301/1‑101 (MQ=255)
                                               gcggcgATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTAt                                                       <  1:8717016/101‑1 (MQ=255)
                                               gcggcgATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTAt                                                       <  1:5876096/101‑1 (MQ=255)
                                                 ggcgATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCa                                                     >  1:10929483/1‑101 (MQ=255)
                                                   cgATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCAcc                                                   >  1:5568237/1‑101 (MQ=255)
                                                       gatTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCtata                                               <  1:8359349/101‑1 (MQ=255)
                                                        atTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATaa                                              <  1:61021/101‑1 (MQ=255)
                                                              cTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCAAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATtg                                         >  1:8224246/1‑101 (MQ=255)
                                                               tGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATtgc                                        <  1:9730655/101‑1 (MQ=255)
                                                                  tttGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATtgctgc                                     >  1:9030338/1‑101 (MQ=255)
                                                                  tttGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATtgctgc                                     >  1:6556076/1‑101 (MQ=255)
                                                                     ggCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCa                                  <  1:5153306/101‑1 (MQ=255)
                                                                     ggCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCa                                  <  1:3642712/101‑1 (MQ=255)
                                                                       cTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCAcg                                >  1:188496/1‑101 (MQ=255)
                                                                        tGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGGCCAATAAAGAAGATGTTAAGCAGGGGTTTATCTCCTATAAGATTGCTGCCCActc                               <  1:9083139/101‑3 (MQ=255)
                                                                        tGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCAcgc                               <  1:10120562/101‑1 (MQ=255)
                                                                        cgcgcGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCAcgc                               >  1:9244847/2‑101 (MQ=255)
                                                                         gcgcGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCAcgcc                              <  1:7158203/101‑1 (MQ=255)
                                                                          cgcgATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCAcgccg                             >  1:4638076/1‑101 (MQ=255)
                                                                           gcgATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCTCAcgccgc                            >  1:12945977/1‑101 (MQ=255)
                                                                           gcgATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCAcgccgc                            >  1:805527/1‑101 (MQ=255)
                                                                           gcgATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCAcgccgc                            >  1:2566400/1‑101 (MQ=255)
                                                                            cgATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCACGCCGCt                           >  1:12978690/1‑101 (MQ=255)
                                                                             gATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCACGCCGCTg                          >  1:792150/1‑101 (MQ=255)
                                                                                gCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCACGCCGCTGAcc                       >  1:6954667/1‑101 (MQ=255)
                                                                                 ctctGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCACGCCGCTGACCt                      >  1:3158145/1‑101 (MQ=255)
                                                                                 ctctGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTGAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCACGCCGCTGACCt                      <  1:4412817/101‑1 (MQ=255)
                                                                                  tctGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCACGCCGCTGACCTg                     >  1:8151819/1‑101 (MQ=255)
                                                                                  tctGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCACGCCGCTGACCTg                     >  1:11785811/1‑101 (MQ=255)
                                                                                   ctGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCACGCCGCTGACCTgg                    <  1:782528/101‑1 (MQ=255)
                                                                                    tGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGc                   <  1:460571/101‑1 (MQ=255)
                                                                                    tGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGc                   >  1:13076141/1‑101 (MQ=255)
                                                                                    tGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGACGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGc                   <  1:4284005/101‑1 (MQ=255)
                                                                                     gTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGCg                  <  1:10325192/101‑1 (MQ=255)
                                                                                     gTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGCg                  <  1:11459147/101‑1 (MQ=255)
                                                                                     gTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGCg                  <  1:14113160/101‑1 (MQ=255)
                                                                                      ttACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGCGa                 <  1:4892930/101‑1 (MQ=255)
                                                                                      ctACGTCACGCCGAAGGAGCACCTGGGTCTGCCGAACCGCGACGACGTCAAGGTCGGCGTGATTACTTACAAGATTTCCGCGCACGCCGCTGACCTGGcaa                 >  1:3518407/2‑99 (MQ=255)
                                                                                       tACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGCGaa                >  1:1279039/1‑101 (MQ=255)
                                                                                       tACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGCGaa                >  1:5293178/1‑101 (MQ=255)
                                                                                        aCGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGCGaaa               >  1:11742968/1‑101 (MQ=255)
                                                                                         cGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGCGAAAg              <  1:11647688/101‑1 (MQ=255)
                                                                                           tAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGCGAAAggg            >  1:5226243/1‑101 (MQ=255)
                                                                                              cGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGCGAAAGGGCAt         <  1:2540590/101‑1 (MQ=255)
                                                                                                ccAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGCGAAAGGGCATcc       <  1:13363460/101‑1 (MQ=255)
                                                                                                ccAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGCGAAAGGGCATcc       <  1:2010024/101‑1 (MQ=255)
                                                                                                 cAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGCGAAAGGGCATCCg      <  1:8731541/101‑1 (MQ=255)
                                                                                                  aaaaGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGCGAAAGGGCATCCgg     >  1:4025958/1‑101 (MQ=255)
                                                                                                     |                                                                                                    
CTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCACGCCGCTGACCTGGCGAAAGGGCATCCGGGCG  >  NC_007779/3441978‑3442178

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.