Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,959,457 T→G 8.0% intergenic (+12/+75) ilvC → / ← ppiC ketol‑acid reductoisomerase, NAD(P)‑binding/peptidyl‑prolyl cis‑trans isomerase C (rotamase C)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,959,4570TG 8.0% 310.5 / 5.2 110intergenic (+12/+75)ilvC/ppiCketol‑acid reductoisomerase, NAD(P)‑binding/peptidyl‑prolyl cis‑trans isomerase C (rotamase C)
Reads supporting (aligned to +/- strand):  ref base T (34/65);  new base G (9/2);  total (43/67)
Fisher's exact test for biased strand distribution p-value = 3.24e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 114 total aligned reads displayed.

GCTATTCCGGAAGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCA  >  NC_000913/3959308‑3959599
                                                                                                                                                     |                                                                                                                                              
gCTATTCCGGAAGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGTTGTTGCGGGTTAAGTGCGCGCTGAt                                                                                                                                                <  1:329358/150‑1 (MQ=255)
gCTATTCCGGAAGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGTTGTTGCGGGTTAAGTGCGCGCTGAt                                                                                                                                                <  1:1141686/150‑1 (MQ=255)
gCTATTCCGGAAGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAt                                                                                                                                                <  2:423117/150‑1 (MQ=255)
gCTATTCCGGAAGGCGCGGTAGATAACGGGCAACTGCGTGAAGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAt                                                                                                                                                >  2:1134355/1‑150 (MQ=255)
 ctcttCCGGAAGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATg                                                                                                                                               >  1:314210/4‑150 (MQ=255)
   aTTCCGGAAGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGcc                                                                                                                                             <  2:2191069/150‑1 (MQ=255)
   aTTCCGGAAGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGcc                                                                                                                                             <  1:97647/150‑1 (MQ=255)
      ccGGAAGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAATCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATgggcgg                                                                                                                                          >  2:1493/1‑145 (MQ=255)
       cGGAAGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATCTGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGGCcaca                                                                                                                                         >  2:911928/1‑150 (MQ=255)
       cGGAAGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCa                                                                                                                                         >  2:1921792/1‑150 (MQ=255)
       cGGAAGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCa                                                                                                                                         >  2:1575095/1‑150 (MQ=255)
       cGGAAGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTGAAGTGCGCGCTGAGGCCCGCa                                                                                                                                         >  2:1477736/1‑150 (MQ=255)
        ggAAGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTAGATGACAGATATGAAACGTATTGCTGTTGCGGGTTCAGTGCGCGCTGAGGcccgccc                                                                                                                                        >  1:476265/1‑146 (MQ=255)
             gcgcGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGa                                                                                                                                   <  2:1842742/150‑1 (MQ=255)
             gcgcGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGa                                                                                                                                   <  1:1063548/150‑1 (MQ=255)
               gcgGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTCTATGACAGATATGAAACGTATTGCTGTTGCGGGTGAAGTGCGCGCTGATGCCCGCACCCCGAcc                                                                                                                                 >  2:1688830/1‑150 (MQ=255)
               gcgGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGAcc                                                                                                                                 <  1:1230118/150‑1 (MQ=255)
               gcgGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGAcc                                                                                                                                 <  2:356955/150‑1 (MQ=255)
               gcgGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGAcc                                                                                                                                 <  1:702050/150‑1 (MQ=255)
               gcgGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGAcc                                                                                                                                 <  1:2398517/150‑1 (MQ=255)
               gcgGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGAcc                                                                                                                                 <  1:1230242/150‑1 (MQ=255)
                   tAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCctct                                                                                                                             <  1:1745973/150‑1 (MQ=255)
                   tAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCctct                                                                                                                             <  1:636901/150‑1 (MQ=255)
                    aGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCtctc                                                                                                                            <  1:2160725/150‑1 (MQ=255)
                     gATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTcc                                                                                                                           <  1:1666602/150‑1 (MQ=255)
                         aCGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCcaca                                                                                                                       <  2:1894653/150‑1 (MQ=255)
                         aCGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCcaca                                                                                                                       <  2:60625/150‑1 (MQ=255)
                           gggCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGGCCGCACCCCCAcccacacccacacg                                                                                                                     >  2:435719/1‑148 (MQ=255)
                           gggCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagg                                                                                                                     >  2:1785097/1‑150 (MQ=255)
                           gggCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagg                                                                                                                     <  1:3650/150‑1 (MQ=255)
                             gCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTccctgtggga                                                                                                                   <  2:2088449/150‑8 (MQ=255)
                             gCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTccctgtggga                                                                                                                   <  2:2354230/150‑8 (MQ=255)
                             gCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCACTCCCACaggga                                                                                                                   >  1:1988992/1‑150 (MQ=255)
                             gCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCGCACCCCCCCCCACACCCACagggg                                                                                                                   >  2:504340/1‑149 (MQ=255)
                               taCTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggaga                                                                                                                 <  2:2121341/149‑1 (MQ=255)
                                aCTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggagag                                                                                                                <  1:2194613/150‑1 (MQ=255)
                                 cTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggagagg                                                                                                               <  2:1584197/150‑1 (MQ=255)
                                 cTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggagagg                                                                                                               <  1:2133567/150‑1 (MQ=255)
                                     gTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAgggagagggaga                                                                                                           <  1:1216653/150‑1 (MQ=255)
                                          gtgAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCcaca                                                                                                                       <  1:303585/133‑1 (MQ=255)
                                          gtgAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCcaca                                                                                                                       >  2:303585/1‑133 (MQ=255)
                                            gAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATCTGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCGCACCCCGACCCACTCCCCCCCGGCGAGGGGGAAAACACt                                                                                                    >  2:82429/1‑150 (MQ=255)
                                            gAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCCACCCACTCCCCCAGGGCGGGGGGGAAAacaca                                                                                                    >  2:1321332/1‑149 (MQ=255)
                                            gAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCTCACCCCGACCCACTCCCACAGGGCGAGGGGGAAAACACt                                                                                                    >  1:1175115/1‑150 (MQ=255)
                                             aaCGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAACCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCAACCCCACCCAAAACCCCACGGAGAGGGAGgaaaaaaaa                                                                                                   >  2:66383/1‑145 (MQ=255)
                                             aaCGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTc                                                                                                   <  1:1252063/150‑1 (MQ=255)
                                             aaCGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTc                                                                                                   <  2:1210151/150‑1 (MQ=255)
                                             aaCGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCACCCCCCCACCCACACCCACCCGGAGCGGGGGAAAacacac                                                                                                   >  1:2032030/1‑148 (MQ=255)
                                               tgAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCaa                                                                                                 <  1:1076055/149‑1 (MQ=255)
                                               cGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCaa                                                                                                 <  2:145057/150‑1 (MQ=255)
                                                 aaGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGGCCCCACCCCCCCCCACACCCCCCCGGCGGGGGGGGAAACACACaaca                                                                                               >  2:1076055/1‑148 (MQ=255)
                                                          gCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCt                                                                                      <  1:816408/150‑1 (MQ=255)
                                                          gCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCt                                                                                      <  1:1951599/150‑1 (MQ=255)
                                                          gCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCACACCCCCACCCTCTACCACAGGGAGAGGGAGAAAACACACAATGCCATCTccg                                                                                      >  1:1390468/1‑149 (MQ=38)
                                                          gCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGAGCGCGCTGATGCCCACACCCCGACCCACACCCACAGGGAGAGGGAGAAAACACAAAAGGCCGtcaaca                                                                                      >  1:1390502/1‑146 (MQ=38)
                                                           cAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTAACCCCGACCCACTCCCACAGGGAGAGGGAGAAAAAAATCAAGGCCATCTACTg                                                                                     >  2:1354617/1‑150 (MQ=38)
                                                             gCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGAGGCCCGCACCCCCACCCACACCCCCACGGAGAGGGCGAAAACACGCAACGCCGTCTCCgggg                                                                                   >  1:1934983/1‑149 (MQ=255)
                                                                aTGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                                                                <  1:1098922/150‑1 (MQ=255)
                                                                 tGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCACTCCCACAGGGAGAGGGAGAAAACACACAAGGCCTTCTCCTGGAGAAg                                                                               >  1:2214801/1‑150 (MQ=255)
                                                                     aTTGAGCAGGTAGGTTAGAAACTGCGCGGCTATCTGACAGATATGAAACGGATTGCTGTTGCGGGGTAAGTGCGCGCTGATGCCCTCCCCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt                                                                           <  1:2189248/150‑1 (MQ=38)
                                                                     aTTGAGCAGGTAGGTACGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt                                                                           <  1:2189250/150‑1 (MQ=255)
                                                                     aTTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACCGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCGCCTGGAGAAGGCCt                                                                           <  1:1127710/150‑1 (MQ=255)
                                                                     aTTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt                                                                           <  1:671726/150‑1 (MQ=255)
                                                                         aGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCATCTCCTGGAGAAGGCCTTGCt                                                                       >  2:1214382/1‑150 (MQ=255)
                                                                           cAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTAt                                                                     <  1:306621/150‑1 (MQ=255)
                                                                            aggtaggtAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCCCCCCACACCCCCAGGGGGGGGGCGAAAACACTCAAAGCCATCTCCAGGCGAAGGGCCGGCTAtt                                                                    >  1:168739/1‑150 (MQ=255)
                                                                             ggtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGACAACCCTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTa                                                                   <  2:794556/150‑1 (MQ=255)
                                                                                   taAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCg                                                             >  2:2185608/1‑150 (MQ=255)
                                                                                   taAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCACTCCCACAGGGCGAGGGAGAAAACACTCAAGGCCGTCTCCTGGAGAAGGCCGTGCTATTAGTTGCg                                                             >  2:1464666/1‑150 (MQ=255)
                                                                                        aaCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACa                                                        <  2:2169754/150‑1 (MQ=255)
                                                                                                ggatataTGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATTGCCTCAACCCGACCCAATACCACAGGGATTTGGGGAAAAAAATCAAAGCCTTATGCTGGAGTATGCCTTTGAATTATTTTAGGTAAAGCACCATa                                                >  2:1253460/4‑150 (MQ=255)
                                                                                                ggCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCCACCCACACCCACACGGAGAGGGGGAAAAAACACAAGGCCATCTCCAGGGGAAAGCCGTGCTATAAGTTGCGGTACAGCAccaaa                                                >  2:946457/1‑147 (MQ=255)
                                                                                                  cTATATGACAGATATGACACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCACGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAat                                              <  1:1298643/150‑1 (MQ=255)
                                                                                                  cTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAat                                              <  1:1401160/150‑1 (MQ=255)
                                                                                                      atGACAGCTATGCAACGGAGTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCCAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAatgatg                                          <  1:2192339/150‑1 (MQ=38)
                                                                                                      atGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAatgatg                                          <  1:1519374/150‑1 (MQ=255)
                                                                                                      atGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAatgatg                                          <  1:546770/150‑1 (MQ=255)
                                                                                                      atGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAatgatg                                          <  1:301892/150‑1 (MQ=255)
                                                                                                      atGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACACGGAGAGGGAGAAAACACTCAACGCCTTCTCCTCGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAatgatg                                          <  1:766517/150‑1 (MQ=255)
                                                                                                      atGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCACTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTAGCTATTAGTTGCGGTACAGCACCTTAatgatg                                          <  1:549765/150‑1 (MQ=255)
                                                                                                       tGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGAtgt                                         <  2:164060/150‑1 (MQ=255)
                                                                                                        gACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGAtgtg                                        <  1:658045/150‑1 (MQ=255)
                                                                                                          cAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGat                                      <  2:967225/150‑1 (MQ=255)
                                                                                                             atatGAAACGTATTGCTGTTGCGGGTTAAGTGCACGCTGAGGCCCGCACCCCCACCCACTACCACAGGGAGAGGGAGAAAAAACACAAGGCCATATCCGGGAGGATGCCGTGGTAATAGTTGGGGGTCAGCCCCATAATTATTTTatata                                   >  1:1496636/1‑149 (MQ=255)
                                                                                                                  aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACTGGGAGAGGGAGAAATCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACTGCTCCTTTATGATGTGATATCCGAAc                              <  1:1993812/150‑1 (MQ=38)
                                                                                                                  aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAAc                              <  2:2316134/150‑1 (MQ=255)
                                                                                                                    aCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCACACCCCCACCCACTCCCACAGGGCGAGGGAGAAAACAATCAAGGCCGTGACCAGGAGAAGGCCGTGCTATTAGTTGCGGTACAGCACCATAATGATGTTTTATCAGAAcag                            >  2:1642148/1‑148 (MQ=255)
                                                                                                                       tATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCAACCCGACCCAATCCCAAAGGTAGTGGGGGGAAACACTCAAAGCCTTCGCCTGGAGTAGGCCGTGCTAATAGTTGCGGTGCAGCCCCCTAAAGATGTGTTATCCGAACggggg                         >  2:674826/1‑149 (MQ=255)
                                                                                                                       tATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCAAACCCCGACCCGCAACCACAGGGATAGGGAGAAAAAAAAAAAGGCCTTCAACATGAGAAGGCCTTTCTAAAAGTTTCGGTCAAGCAACATAAAGATTTGATATACGAACAGGgt                         >  2:481176/1‑150 (MQ=255)
                                                                                                                             tGTTGCGGTTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACTGGGTGTGGGAGAAAACACTCACGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGTTTTGATATCCGAACTGGGTGTGccg                   <  2:168739/150‑3 (MQ=38)
                                                                                                                                tGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGATGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCgg                <  1:1575095/150‑1 (MQ=255)
                                                                                                                                tGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCgg                <  1:911928/150‑1 (MQ=255)
                                                                                                                                tGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCgg                <  1:1921792/150‑1 (MQ=255)
                                                                                                                                   gggTTAAGTGCGCGCTGATGCCCACACCCCCACCCACTCCCACAGGGAGAGGGAGAAAACAAGCAAGGCCCTCTCCGGGCGAACGCCGTGCGATTAGTTGCGGTACAGCACCATAATGATGTGATAACCCAACAGGGGGTGCAGCGGGcc             >  2:1714069/1‑150 (MQ=25)
                                                                                                                                     gTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCgg           <  1:2001813/150‑1 (MQ=255)
                                                                                                                                     gTTAAGTGCGCGCTGATGCCCGCACCCCCACCCGCACCCACAGGGAGAGGGAGAAAAAACACAACGCCATAACCAGGAGAAGGCCGTGCTATAAGTTGCGGTACACCACCATAAAGATGTGATATACGAACGGGGTGTGCCCCGGGCCgg           >  2:1915397/1‑150 (MQ=19)
                                                                                                                                       tAAGTGCGCGCTGAGGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTc         <  2:630214/150‑1 (MQ=255)
                                                                                                                                         aGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGTGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCgg       <  1:1493/150‑1 (MQ=255)
                                                                                                                                           tGCGCGCTGATGCCCTCACCCCGACCCTCTACCACAGGGAGTTGGAGAAAACACTCAAGGCCTTCTACTGGAGAAGGCCTTGCTAATAGTTTCGGTACAGCACCATAATGATGTGAAATACGAACAGGGTGTGCAGCGGGCCCGTCGGCt     >  2:1963235/1‑150 (MQ=37)
                                                                                                                                              gcgcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCa  <  1:1069276/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
GCTATTCCGGAAGGCGCGGTAGATAACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCA  >  NC_000913/3959308‑3959599

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: