Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,781,022 T→G 7.6% intergenic (+5/‑193) 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,0220TG 7.6% 2304.8 / 45.7 807intergenic (+5/‑193)lldD/trmLL‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (176/549);  new base G (75/7);  total (251/556)
Fisher's exact test for biased strand distribution p-value = 7.79e-34
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 835 total aligned reads displayed.

ATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGC  >  NC_000913/3780872‑3781171
                                                                                                                                                      |                                                                                                                                                      
  cggATCGACAAAGAGATGACCGTGGCGATGACGCTGACTGGCGCGAAATCGAGCAGCGAACTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGttt                                                                                                                                                       <  2:2462899/148‑1 (MQ=255)
  cTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGttt                                                                                                                                                       <  2:8238591/150‑1 (MQ=255)
   tGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTg                                                                                                                                                      >  2:5307295/1‑150 (MQ=255)
    gATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTgc                                                                                                                                                     >  1:7653325/1‑150 (MQ=255)
    gATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTgc                                                                                                                                                     <  2:1241707/150‑1 (MQ=255)
      tCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTgcgc                                                                                                                                                   >  1:1142086/1‑150 (MQ=255)
        gAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGccc                                                                                                                                                 >  2:702617/1‑150 (MQ=255)
         aaaaaGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGcccc                                                                                                                                                >  1:839129/1‑150 (MQ=255)
           aaaGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCt                                                                                                                                              <  2:5916963/150‑1 (MQ=255)
             agagATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGcccccacc                                                                                                                                            >  1:4068714/1‑147 (MQ=255)
              gagaTGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCAc                                                                                                                                           <  2:7378283/150‑1 (MQ=255)
              gagaTGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCAc                                                                                                                                           <  1:8175349/150‑1 (MQ=255)
               agaTGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCAcc                                                                                                                                          <  2:6852904/150‑1 (MQ=255)
                  tGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCCCCCTa                                                                                                                                       >  1:5049233/1‑150 (MQ=255)
                  tGAAAGTGGCGATGACGCTGACTGGCGCGAAACCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCCCAcccaa                                                                                                                                       >  1:1689849/1‑148 (MQ=255)
                    acaGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAAc                                                                                                                                     <  1:8953713/148‑1 (MQ=255)
                     aaGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAAcc                                                                                                                                    <  1:5909113/150‑1 (MQ=255)
                        tGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCCCCCCCAAAccccc                                                                                                                                 >  2:6568165/1‑148 (MQ=255)
                          gCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCtctc                                                                                                                               <  2:9215026/150‑1 (MQ=255)
                           cGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTcc                                                                                                                              <  2:4535093/150‑1 (MQ=255)
                            gATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGACTGCGGCATAGTCGTGTGCGCCCCCCCCCCCCAACcccccccc                                                                                                                             >  1:2035869/1‑150 (MQ=255)
                             aTGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCt                                                                                                                            <  1:7875831/150‑1 (MQ=255)
                             aTGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCt                                                                                                                            <  1:2282764/150‑1 (MQ=255)
                             aTGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCt                                                                                                                            <  2:4867258/150‑1 (MQ=255)
                               gACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCa                                                                                                                          <  1:8842720/150‑1 (MQ=255)
                               gACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCa                                                                                                                          <  1:7356937/150‑1 (MQ=255)
                               gACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCa                                                                                                                          <  2:2059130/150‑1 (MQ=255)
                               gACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCa                                                                                                                          <  2:7508268/150‑1 (MQ=255)
                                 cGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAgg                                                                                                                        <  2:5800286/150‑1 (MQ=255)
                                     gACTGGCGCGAAATCGATCAGCGAAAGTACGCAAGAGTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGgaga                                                                                                                   <  2:8486417/150‑1 (MQ=255)
                                       cTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCTGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGATTGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAgg                                                                                                                 <  1:1843528/150‑1 (MQ=255)
                                       cTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTGAGGGAGAgg                                                                                                                 <  2:2667115/150‑1 (MQ=255)
                                        tGGCGCGATATCGATCTGCGTAATTACTCAAGTTTCTCTGGTGTAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAATGGATTTCGGCTTAGTCGTTTGCGCACCCTCACCATAACCCTCTCCCTGAGGtagttgg                                                                                                                <  1:8709552/150‑5 (MQ=255)
                                        tGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAggg                                                                                                                <  1:2170310/150‑1 (MQ=255)
                                        tGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCCCACCCAAAACCACAACCACACGGGGAggg                                                                                                                >  2:8710937/1‑150 (MQ=255)
                                          gcgcGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGAACCCAAACCCACACCCCCAGGGGGAggggg                                                                                                              >  1:5684146/1‑149 (MQ=255)
                                           cgcgAAATCGATCAGCGACATTACGAAATATTCGCTGTTGCAGTGTCTGGGTACAGATTTTCCTGCGGCACTGTATCCCATGGCGAAAGGGAATGCGGCATATTCGTTTGCGACCACTCAAACTAACCCTATACCTCATGGAGAggttac                                                                                                             <  2:2828161/150‑5 (MQ=255)
                                           cgcgAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAAcc                                                                                                                                    <  1:8621938/128‑1 (MQ=255)
                                           cgcgAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGAc                                                                                                             <  2:3298925/150‑1 (MQ=255)
                                           cgcgAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCCCACCCAAACCCAAACCCACAGGGCGAgggggc                                                                                                             >  1:7569070/1‑148 (MQ=255)
                                           cgcgAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCCCACCCCAACCCCCGCCCGCCGGGAgggggggc                                                                                                             >  2:3010153/1‑148 (MQ=255)
                                                   cGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGgc                                                                                                     <  1:7990728/150‑1 (MQ=255)
                                                      tCAGCGAACTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCt                                                                                                  <  2:5232620/150‑1 (MQ=255)
                                                      tCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGTTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCt                                                                                                  <  1:2864871/150‑1 (MQ=255)
                                                      tCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGTGGGGACCGTTCGGCGCt                                                                                                  <  1:4499988/150‑1 (MQ=255)
                                                      tCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCt                                                                                                  <  2:3724789/150‑1 (MQ=255)
                                                      tCAGCGAAAGTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCt                                                                                                  <  2:8417418/150‑1 (MQ=255)
                                                       cAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtg                                                                                                 <  2:5014488/150‑1 (MQ=255)
                                                         gCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgta                                                                                               <  1:3885575/150‑1 (MQ=255)
                                                          cGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCACACCCAAAACCACTCCCGAAGGGGGAGGGGGCCGATCGGGGCagtat                                                                                              >  2:1504381/1‑150 (MQ=255)
                                                           gAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCtgtatg                                                                                             <  2:1734101/150‑1 (MQ=255)
                                                             aaTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCACACCCAAACCCACAACCGCAGGGAGAGGGGACCGTTCGGCGCtgtatgta                                                                                           >  2:7714469/1‑150 (MQ=255)
                                                                  cGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCACCCGAACCCACACCCGCAGGGCGAGGGGACCGGGCGGCGCGGTATGTACTccc                                                                                      >  1:8256407/1‑150 (MQ=255)
                                                                    cAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCCCCCCCAAAACCACAACCACAGGGGGTGGGGGCCCGTCGGGGCGGTATTTAATccccc                                                                                    >  1:3471651/1‑148 (MQ=255)
                                                                      aGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCACAACCAAACCCACTCCCAAAGGGAGGGGGGGACCGTCGGCGCAGTATTTTCTACCacac                                                                                  >  1:6095760/1‑150 (MQ=255)
                                                                            gCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGaa                                                                            <  1:7113643/150‑1 (MQ=255)
                                                                             cTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGaaa                                                                           <  2:1102036/150‑1 (MQ=255)
                                                                             cTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGaaa                                                                           <  2:6200181/150‑1 (MQ=255)
                                                                             cTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGaaa                                                                           <  1:6250599/150‑1 (MQ=255)
                                                                              tGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCCCCCCCAAACCCACACCCACAGGGGGAGGGGGCCCGTGGGCGCTGTATTTACAACCGCACTCAGAAAc                                                                          >  1:554621/1‑150 (MQ=255)
                                                                                 tGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGac                                                                       <  2:8836845/150‑1 (MQ=255)
                                                                                     ggggCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCACCACCAAACCCACGACCCCAGGGAGAGGGGGCCCGTCGGCGCGGTTTGTACTCCCACACACTGTAACGAAAccc                                                                   >  1:6633721/1‑149 (MQ=255)
                                                                                       ggCTGGGTAAAGAGTTGCCTGCGGCAGTGGCTCCCTTGGCTAATGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCTGGGAGTGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCa                                                                 <  1:6736291/150‑1 (MQ=255)
                                                                                       ggCTGGGTAAAGAGGTGCCTGCGGCACTGGCTCCCAGGGCGAAAGGGAATGCGGCAGAGTCGTTTGCGCCCCCTCACCCTACCCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCCCTCTGAAACGACACCGCa                                                                 <  2:1019562/150‑1 (MQ=255)
                                                                                        gCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCCCCCCCAAAACCACAACCACAGGGACGGGGGACCCGTCGGCGCGGTGTTTACAACCACACACAGAAAACACAcccccc                                                                >  1:4797223/1‑146 (MQ=255)
                                                                                          tGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCActc                                                              <  2:3929343/150‑1 (MQ=255)
                                                                                               aaaGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCttttt                                                         <  2:6697884/150‑1 (MQ=255)
                                                                                                    gTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtctc                                                    <  1:1318254/150‑1 (MQ=255)
                                                                                                     ttGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCCCACCCGAAACCACACCCACACGGGGGGGGGGCCCGTGGGCGCGGGATGTTATACCACACCCAGAAACGAAACCCCCCACATTTTTTATcc                                                   >  1:6918612/1‑150 (MQ=255)
                                                                                                        ccTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTc                                                <  2:4679299/150‑1 (MQ=255)
                                                                                                        ccTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCGCCCGAAACCACTCCCACAGGGGGTGGGGGCCGTGGGGGGCGGTGTGGGATCCCGACCTCAGAAACGGCACCGGCCGCgttttttatgccgc                                                >  2:7929394/1‑143 (MQ=255)
                                                                                                         cTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCg                                               <  1:48522/150‑1 (MQ=255)
                                                                                                               cACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCCCCCCCAAAACCCCACCCGCCCGGGGGGGGGGCCCGTCGGCCCGGGTTGTACAACCGCACACGTAAACAACACCCCACACAGTTTTTGTGCCACGcccccc                                         >  2:5476727/1‑148 (MQ=255)
                                                                                                                  tGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCgg                                      <  1:1450988/150‑1 (MQ=255)
                                                                                                                  tGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCgg                                      <  1:2548898/150‑1 (MQ=255)
                                                                                                                     cTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTTTTTTTTTCTCCCTCGCCCCTCCGGGga                                   <  1:1005444/150‑1 (MQ=255)
                                                                                                                     cTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGga                                   <  1:3096824/150‑1 (MQ=255)
                                                                                                                     cTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGga                                   <  2:6380895/150‑1 (MQ=255)
                                                                                                                     cTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGga                                   <  1:3266004/150‑1 (MQ=255)
                                                                                                                     cTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGga                                   <  1:6487754/150‑1 (MQ=255)
                                                                                                                       cccATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGCAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGgaga                                 <  1:3023228/150‑1 (MQ=255)
                                                                                                                       cccATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGgaga                                 <  2:8953855/150‑1 (MQ=255)
                                                                                                                        ccATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGagag                                <  2:5110290/150‑1 (MQ=255)
                                                                                                                          aTGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAggg                              <  2:156297/150‑1 (MQ=255)
                                                                                                                          aTGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAggg                              <  1:5408677/150‑1 (MQ=255)
                                                                                                                             gCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCCCACCCAAAACCAAACCCACCGGGAGGGGGGACCCTTGGGCGCGGTTTGTTCAACCGAACACGGAAACCGCCCCCCAATAATTTTTTGACCCACGGCCCTCCGGGGAGAGGGCCg                           >  2:5621414/1‑150 (MQ=255)
                                                                                                                              cGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCCCCCTACCCCGCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCCCTCTGACCCGCCCCCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCgg                          <  1:4640946/150‑1 (MQ=255)
                                                                                                                               gAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGACCggg                         <  2:2173946/150‑1 (MQ=255)
                                                                                                                                aaaGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCgggg                        <  1:6373501/150‑1 (MQ=255)
                                                                                                                                 aaGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGt                       <  1:5204041/150‑1 (MQ=255)
                                                                                                                                   gggAATGCGGCATAGTCGTTTGCGCCCCCTCACCCAAACCCTCCCCCTCTGGGAGAGGGGACCTTTCGGCGCTGTTTGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGa                     <  1:2452912/150‑1 (MQ=255)
                                                                                                                                    ggAATGCGGCATAGTCGTTTGCGCCCCCTCTCCCTAACCCTCTCCCTCAGGGAGAGGTGTCCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGCCACCGCACTCTTTTTTTCTCCCTCGCCCCTCCTGGGAGAGGGCCGGGGTGAg                    <  2:3494598/150‑1 (MQ=255)
                                                                                                                                    ggAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAg                    <  2:7638906/150‑1 (MQ=255)
                                                                                                                                     gattGCGGCATAGTCGTTGGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGTGGGGACCGTTCGGCGCTGTTTGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgg                   <  1:4947525/147‑1 (MQ=255)
                                                                                                                                     gAATGCGGCATAGTCGTGTGCGCCCCCTCACCCTACCCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgg                   <  2:7198566/150‑1 (MQ=255)
                                                                                                                                       aTGCGGCATAGTCGTTTGCGCCCCCCCACCCAAAACCACTCCCACAGGGAGAGGGGTCCCGTGGGCGCGGTTTTTAAACACCCACACAGAAACCCCACCCCACAAATTTTTTATACCAAGCCCCACCCGGGAGAGGGCCCgggggggggg                 >  2:402814/1‑150 (MQ=255)
                                                                                                                                           ggCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCTCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaaa             <  2:9111082/150‑1 (MQ=255)
                                                                                                                                              attttcttttGCGCCCCCACACCATAACCCTCACCATCTGGTTTTGGGTACCTTTCGGCGGTTTTTTTACTCCATCACTCTTAAACGACACGGAAATTTTTTTTTCTCCCTCGCCCCTCGGGGTAGTGGGCGGGGTTGTGGGTAAAAGGc          <  2:3090739/143‑1 (MQ=255)
                                                                                                                                                    gTTTGCGCCCCCTCGCCCTAACCCGCTCCCGCCGGGAGCGGGGCCCGTTCGGCGCGGTTTGTACGCCCTCACTCGGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGGGGGCCGGGGTGCGGGGAAAAGGCCGCACt    >  2:3623968/1‑150 (MQ=255)
                                                                                                                                                     tttGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCCCTg   >  2:4327250/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
ATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGC  >  NC_000913/3780872‑3781171

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: