Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,781,038 T→A 13.4% intergenic (+21/‑177) lldD → / → trmL L‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,781,0380TA13.4% 589.4 / 10.7 219intergenic (+21/‑177)lldD/trmLL‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (16/154);  new base A (22/2);  total (62/157)
Fisher's exact test for biased strand distribution p-value = 9.19e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 219 total aligned reads displayed.

TGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTG  >  NC_000913/3780890‑3781185
                                                                                                                                                    |                                                                                                                                                     
tGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTa                                                                                                                                                      >  1:3221427/1‑150 (MQ=255)
 gAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTaa                                                                                                                                                     <  1:1175953/150‑1 (MQ=255)
  aaaGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAAc                                                                                                                                                    <  2:1785670/150‑1 (MQ=255)
   aaGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAAcc                                                                                                                                                   <  1:4924095/150‑1 (MQ=255)
   aaGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAAcc                                                                                                                                                   <  1:2196269/150‑1 (MQ=255)
     gTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCct                                                                                                                                                 <  1:3291638/150‑1 (MQ=255)
        gCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCtctc                                                                                                                                              <  2:5091661/150‑1 (MQ=255)
          gATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCCCACCCCAAaccacaacc                                                                                                                                            >  2:4811438/1‑144 (MQ=255)
            tGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTc                                                                                                                                          <  1:4993895/150‑1 (MQ=255)
             gACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAATGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCCCCCGCACCCACAcccccc                                                                                                                                         >  2:2792068/1‑147 (MQ=255)
             gACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCa                                                                                                                                         <  2:5357051/150‑1 (MQ=255)
             gACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCa                                                                                                                                         <  2:3894859/150‑1 (MQ=255)
              aCGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑g                                                                                                                                       <  1:3883238/150‑1 (MQ=255)
               cGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTGA‑gg                                                                                                                                      <  1:2976930/150‑1 (MQ=255)
                 cTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGGTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGGGTGCCCCCCCCCCCCCCCACCCACTcccccccggc                                                                                                                                     >  2:3483348/1‑143 (MQ=255)
                   gACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGgaga                                                                                                                                  <  1:553606/150‑1 (MQ=255)
                   gACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCACCCCCGAACCCACACCCACA‑GGgagt                                                                                                                                  >  1:4504247/1‑149 (MQ=255)
                   gACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCACCCCCCAAACCACAACCCCA‑gggggt                                                                                                                                  >  2:4230222/1‑147 (MQ=255)
                    aCTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGagag                                                                                                                                 <  1:1920774/150‑1 (MQ=255)
                     cTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTGA‑GGGAGAgg                                                                                                                                <  1:2498364/150‑1 (MQ=255)
                     cTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTGA‑GGGAGAgg                                                                                                                                <  2:2254349/150‑1 (MQ=255)
                     cTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAgg                                                                                                                                <  1:4055032/150‑1 (MQ=255)
                      tGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAggg                                                                                                                               <  1:5402464/150‑1 (MQ=255)
                        gcgcGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGa                                                                                                                             <  1:1646695/150‑1 (MQ=255)
                        gcgcGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGa                                                                                                                             <  1:3371994/150‑1 (MQ=255)
                        gcgcGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCCCCCCCAAAAACACACCCACA‑CGGAGAggggg                                                                                                                             >  1:2065403/1‑149 (MQ=255)
                         cgcgAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGAc                                                                                                                            <  2:1252504/150‑1 (MQ=255)
                         cgcgAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCGCACCCAAACCCACACCCGCA‑gggggggggggc                                                                                                                            >  2:5283223/1‑148 (MQ=255)
                          gcgAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTGA‑GGGAGAgggttag                                                                                                                           <  2:794925/150‑5 (MQ=255)
                          gagAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTGA‑GGGAGAgggttag                                                                                                                           <  2:3221427/148‑5 (MQ=255)
                           cgAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCCCCCCCAAAACCACACCCAAC‑GGGCGGGGGGAccc                                                                                                                          >  1:1048209/1‑149 (MQ=255)
                           cgAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCGCACCCAAAACCCCACCCACA‑CGGCGAGGGGACCg                                                                                                                          >  1:4411104/1‑150 (MQ=255)
                             aaaTCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGtt                                                                                                                        <  2:3061837/150‑1 (MQ=255)
                               aTCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCg                                                                                                                      <  2:1192791/150‑1 (MQ=255)
                                tCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGTGTTGCCTGCGGCACTGGCTCCCATGGCTAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCgg                                                                                                                     <  1:106790/150‑1 (MQ=255)
                                tCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCgg                                                                                                                     <  1:1378163/150‑1 (MQ=255)
                                   aTCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCAAACCCAAACCCACTACCTCA‑GGGAGAGGGGACCGTTCGgcgc                                                                                                                  >  2:1379039/1‑150 (MQ=255)
                                      aGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTTCGCCCCCCCCCCCCAAACCCCAACCCCA‑GGGGGGGGGGGCCCGTCGgcgccgg                                                                                                               >  2:1168956/1‑147 (MQ=255)
                                      aGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTTCGCCCCCACACCCAAAACCACAACCACA‑GGGAGAGGGGACCCTTCGGCGCagt                                                                                                               >  2:3467654/1‑150 (MQ=255)
                                      aGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCGCACCCGAACCCACACCCGCA‑GGGAGAGGGGCCCCTTCGGCGCggt                                                                                                               >  2:4063004/1‑150 (MQ=255)
                                       gCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCACACCCAAACCCACTCCCACA‑GGGAGAGGGGGCCGTTCGGCGCtgta                                                                                                              >  1:1837414/1‑150 (MQ=255)
                                        cGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCtgtat                                                                                                             <  2:411534/150‑1 (MQ=255)
                                           aaTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCtgtatgta                                                                                                          <  2:3785217/150‑1 (MQ=255)
                                              tACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTc                                                                                                       <  1:5011218/150‑1 (MQ=255)
                                                cGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCACACCCAAACACACACCCACA‑CGGGGGGGGGGCCCGTCGGCGCGTTATGTAcaacc                                                                                                     >  1:3192942/1‑146 (MQ=255)
                                                 gCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTGA‑g                                                                                                                                       <  2:4370468/115‑1 (MQ=255)
                                                 gCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCt                                                                                                    <  1:5176908/150‑1 (MQ=255)
                                                   aaGATTCGATGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGTCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCAACCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTTCTCCCTCa                                                                                                  <  2:459222/150‑1 (MQ=255)
                                                    aGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCTCCCCCTAACCCTAACCCTCA‑GGGGGAGGGGACCCTTCGGCGCTGTATGTTCTCCCTCAc                                                                                                 >  1:3048118/1‑150 (MQ=255)
                                                     gATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCAct                                                                                                <  1:1778339/150‑1 (MQ=255)
                                                       ttCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCACACCCAAACCCTCGCCCGCA‑GGGAGAGGGGACCCTTCGGCGCTGTATGTAATCCCACAcact                                                                                              >  2:3661031/1‑150 (MQ=255)
                                                         cGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGa                                                                                            <  2:5383375/150‑1 (MQ=255)
                                                                gCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGaca                                                                                     <  2:1282730/150‑1 (MQ=255)
                                                                 cAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGacac                                                                                    <  2:4668808/150‑1 (MQ=255)
                                                                    gggCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCAAAACCAAAACCAAAACCAAA‑GGGAGAGGGGAACCTTAGGCGCAGTATGTACACCCACACACTGAAACGACAcccc                                                                                 >  2:415574/1‑148 (MQ=255)
                                                                     gggtGGGTAATGTGTTGCCTGCGGCACTGGCTCCCTTGGCTATTGGGTTTGCGGCATTGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCT‑GGGTGTGGGGACCGTTCGGCGGTTTTTTTACTCCCTCACTCTGAAACGACACCGCa                                                                                <  1:3748319/147‑1 (MQ=255)
                                                                     ggCTGGGTAATGAGTTGCCTGCGGCACTGGCTCCCATGGCGACAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCACCCAAAACCGCTACCTCC‑AGGGGTGGGGCCCGTTCGGCGCAGTATGTACTACCGCCCACAGCAACGACACCCCa                                                                                >  1:5355804/1‑150 (MQ=255)
                                                                        tGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCCCCCCCGAAACCCAAACCACA‑GGGAGAGGGGACCCCTAGGCGCGGTATTTAAACACACACACAGAAACCACAccccccac                                                                             >  2:3941882/1‑144 (MQ=255)
                                                                          ggTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCtt                                                                           <  2:297046/150‑1 (MQ=255)
                                                                          ggTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCCCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCtt                                                                           <  2:2578880/150‑1 (MQ=255)
                                                                            tAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCtttt                                                                         <  2:4462518/150‑1 (MQ=255)
                                                                             aaaGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCttttt                                                                        <  1:1518058/150‑1 (MQ=255)
                                                                              aaGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCACCCGAAACCCCACCCACA‑CGGGGAGGGGGCCCGGCGGCGCGGGTTGTTGAACCACACGCTGGAAAGACACCCCCCACtttttt                                                                       >  1:394599/1‑150 (MQ=255)
                                                                               agagTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCttttttt                                                                      <  1:143419/150‑1 (MQ=255)
                                                                               agagTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCttttttt                                                                      <  1:857508/150‑1 (MQ=255)
                                                                               agagTTGCCTGCGGCACTGGCGCCCAGGGCGACAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTAGGTACTCCCTCACTCTGAAACGACACCGCACTCttttttt                                                                      <  1:2202964/150‑1 (MQ=255)
                                                                                  gTTGCCTGCGGCAGTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtctc                                                                   <  1:4160315/150‑1 (MQ=255)
                                                                                  gTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtctc                                                                   <  1:4811452/150‑1 (MQ=255)
                                                                                  gTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtctc                                                                   <  1:765853/150‑1 (MQ=255)
                                                                                     gCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCt                                                                <  1:2396102/150‑1 (MQ=255)
                                                                                      ccTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTc                                                               <  2:2761367/150‑1 (MQ=255)
                                                                                       cTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCg                                                              <  1:3548565/150‑1 (MQ=255)
                                                                                            gCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑t                                                        <  1:4665806/150‑1 (MQ=255)
                                                                                                tGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCCgg                                                    <  2:4852780/150‑1 (MQ=255)
                                                                                                 ggCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGACACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCCggg                                                   <  2:2347481/150‑1 (MQ=255)
                                                                                                 ggCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCCggg                                                   <  1:421941/150‑1 (MQ=255)
                                                                                                    tCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCACCCCCAAACCCAAGACCGCC‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACAACCACACAAAGAAACGACACCCCACACATTTTTTAACCCACGCCCC‑ACCCGGgag                                                >  2:5000626/1‑150 (MQ=255)
                                                                                                     cccATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCCGGGgaga                                               <  1:2433239/150‑1 (MQ=255)
                                                                                                      ccATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCCGGGGagag                                              <  2:5381310/150‑1 (MQ=255)
                                                                                                       cagggCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCCGGGGAGAgg                                             <  1:1884528/147‑1 (MQ=255)
                                                                                                       cATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCCGGGGAGAgg                                             <  1:2580567/150‑1 (MQ=255)
                                                                                                        aTGGCGAAAGGGAATGCGGCATAGTCATTTGCGCCCCCTCACCCAAACCCTCTCCCTCA‑GGGTGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCCGGGGAGAggg                                            <  1:4148204/150‑1 (MQ=255)
                                                                                                           gCGAAAGGGAATGCGGCATCGTCGTTTGCGCCCCCGCCCCCTCACCCTCTCCCTCC‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAACCGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCCGGGGAGAGGGCCg                                         <  2:3212107/150‑1 (MQ=255)
                                                                                                           gCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCCGGGGAGAGGGCCg                                         <  1:1165315/150‑1 (MQ=255)
                                                                                                             gAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCTCCCC‑GGAGGGGCGAgggagaaa                                       <  2:576631/150‑6 (MQ=255)
                                                                                                              aaaGGGAATGCGGCATAGTCGTTTGCGCCCCCCCACCCCAAACCACAACCAAACGGGGGGGGGGCCCGTAC‑GCGCGGTTTGTTCGCCCACACACAGAAACGAAACCCCACAAATTTTTTAGACCACCCCCC‑CCCCGGGGGAGGGGCgggg                                      >  2:2147661/1‑150 (MQ=255)
                                                                                                               aaGGGAATGCGGCATAGTCGTTTGCGCCCCCACACCCGAAACCACTCCCGCA‑GGGAGAGGGGGCCCGTGGGCGCTGTGTGTAATCCCGCACGATGAAACGACACCGCCCACTTTTTTTTTACCCCGCCCCGCCCGGGGGGCGGGCGgggg                                      >  2:4901423/1‑150 (MQ=255)
                                                                                                                  ggAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTACCCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCACACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCCGGGGAGAGGGCCGGGGTGAg                                  <  1:2775131/150‑1 (MQ=255)
                                                                                                                  ggAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCCGGGGAGAGGGCCGGGGTGAg                                  <  2:2915797/150‑1 (MQ=255)
                                                                                                                  ggAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCCGGGGAGAGGGCCGGGGTGAg                                  <  1:3630399/150‑1 (MQ=255)
                                                                                                                   gAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCCGGGGAGAGGGCCGGGGTGAgg                                 <  2:4691986/150‑1 (MQ=255)
                                                                                                                    aaTGCGGCATCGTCGTTTGCGCCCCCTCACCCTCCCCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCCGGGGAGAGGGCCGGGGTGAggg                                <  2:3069699/150‑1 (MQ=255)
                                                                                                                       gCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCCGGGGAGAGGGCCGGGGTGAGGGGaa                             <  2:1977163/150‑1 (MQ=255)
                                                                                                                        cGGCATAGTCGTTTGCGCCCCCTCACCCTCACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCCCTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCCGGGGAGAGGGCCGGGGTGAGGGGaaa                            <  2:2399066/150‑1 (MQ=255)
                                                                                                                                          ccccTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCa          >  2:4780783/1‑150 (MQ=255)
                                                                                                                                             cTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAtt       >  1:3304708/1‑150 (MQ=255)
                                                                                                                                                aCCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTc    >  1:3130488/1‑150 (MQ=255)
                                                                                                                                                 cccTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCt   >  1:5223571/1‑150 (MQ=255)
                                                                                                                                                  ccTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTg  >  1:1095656/1‑150 (MQ=255)
                                                                                                                                                  ccTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTg  >  1:4322308/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCC‑TCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTG  >  NC_000913/3780890‑3781185

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: ‑