Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,781,088 T→A 6.6% intergenic (+71/‑127) 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,0880TA 6.6% 809.7 / 3.4 255intergenic (+71/‑127)lldD/trmLL‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (104/120);  new base A (15/4);  total (130/125)
Fisher's exact test for biased strand distribution p-value = 7.81e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 259 total aligned reads displayed.

AAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGA  >  NC_000913/3780941‑3781237
                                                                                                                                                   |                                                                                                                                                     
aaGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCa                                                                                                                                                     <  1:1943626/150‑1 (MQ=255)
aaGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCACACCCAAAACCACAACCACAGGGGGAGGGGGCCGGTAGGCGCGGTATTTACTCCCGCa                                                                                                                                                     >  1:2206200/1‑150 (MQ=255)
    ttCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCActct                                                                                                                                                 <  2:397707/150‑1 (MQ=255)
      cGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGa                                                                                                                                               <  1:1546089/150‑1 (MQ=255)
       gCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCACACCCAAAACCACTACCACAGGGAGAGGGGGACCTTAGGCGCTGTATGTAAAACCACACTCAGaa                                                                                                                                              >  1:2766189/1‑150 (MQ=255)
        cTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGaaa                                                                                                                                             <  1:1129453/150‑1 (MQ=255)
        cTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGaaa                                                                                                                                             <  1:3506701/150‑1 (MQ=255)
        cTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGaaa                                                                                                                                             <  2:2747814/150‑1 (MQ=255)
        cTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGaaa                                                                                                                                             <  1:3975523/150‑1 (MQ=255)
        cTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGaaa                                                                                                                                             <  1:1631082/150‑1 (MQ=255)
         tGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAAc                                                                                                                                            <  2:166797/150‑1 (MQ=255)
          ggTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCGCACCCAAAACCAAAACCACAGGGAGAGGGGGCCCTTAGGCGCGGTATTGACTCCCACACAAAGAAACg                                                                                                                                           >  2:2868531/1‑150 (MQ=255)
           gggctggggGTGGGTTTAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGATTGCGGCATTGTCGTTTGCGCCCCCTCAGGCTAACCCTCTCCCTCAGGGTGTGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGa                                                                                                                                          <  2:2185211/145‑1 (MQ=255)
           gTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTTCGCCCCCCCACCCAAACCCACACCCGCACGGAGAGGGGCCCCTTCGGCGCGGTATGTAAACCCACACACTGAAACGa                                                                                                                                          >  2:704159/1‑150 (MQ=255)
              cAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGacac                                                                                                                                       <  2:3743798/150‑1 (MQ=255)
              cAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCCCACCCCAACCCGCTCCCCCAGGGGGAGGGGCCCCTTCGGCGCGGTGTGTACACCCGCACACAGCAACCacac                                                                                                                                       >  2:4124418/1‑150 (MQ=255)
                   gCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCACCCAAAACCACACCCACAGGGCGCGGGGCCCCTTCGGCGCGGTATGTACACCCACACGCTGAAACGACACCCCAc                                                                                                                                  >  2:2410162/1‑150 (MQ=255)
                         tAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCCCACCCCAAACCACACCCACAGGGGGAGGGGGCCGGTCGGCGCATTATGTACAACCACACACTGAAACAACACCCCACACattt                                                                                                                            >  1:2858133/1‑150 (MQ=255)
                         tAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCCCACCCAAAACCACAACCACACGGGGTGGGGACCCGTAGGCGCGGTATGTTATGACCCACACAGAAACCACCCCGCAcaaattt                                                                                                                            >  1:240737/1‑144 (MQ=255)
                          aaaGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCttttt                                                                                                                           <  2:1399058/150‑1 (MQ=255)
                            agagTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTTTGTACTCCCTCACTCTGAAACGACACCGCACTCttttttt                                                                                                                         <  1:2622841/150‑1 (MQ=255)
                            agagTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCttttttt                                                                                                                         <  2:2094746/150‑1 (MQ=255)
                              agTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGTGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtct                                                                                                                       <  1:1836841/150‑1 (MQ=255)
                                ttGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTcc                                                                                                                     <  2:4027919/150‑1 (MQ=255)
                                ttGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCCCCCAAAACCACACACGCCGGGCGGGGGGGCCCGTCGGGGGGGGTTGTTCACCCGCACACTGTAACCACACCCCCCACTTTTTTTGTcc                                                                                                                     >  2:212573/1‑150 (MQ=255)
                                   ccTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGTGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTc                                                                                                                  <  1:3289342/150‑1 (MQ=255)
                                    cTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCg                                                                                                                 <  2:1206064/150‑1 (MQ=255)
                                     tGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGc                                                                                                                <  1:4073195/150‑1 (MQ=255)
                                             tGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCgg                                                                                                        <  1:2331537/150‑1 (MQ=255)
                                              ggCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCggg                                                                                                       <  1:3650891/150‑1 (MQ=255)
                                               gCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTTCGCCCCCGCACCCAAACCCACACCCGCAGGGCGAGGGGTCCCGTCGGCGCGGTAAGTACAGCCACACGAGGAAACCACACCCCCAACTTTTTTTCTGCCAAGCCCCACCcggg                                                                                                      >  2:2534968/1‑150 (MQ=255)
                                                cTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGga                                                                                                     <  2:3826074/150‑1 (MQ=255)
                                                cTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGga                                                                                                     <  1:1594331/150‑1 (MQ=255)
                                                 tCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGgag                                                                                                    <  2:2597683/150‑1 (MQ=255)
                                                 tCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCCAACCCTCTCCCTCCGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGgag                                                                                                    <  1:3953304/150‑1 (MQ=255)
                                                  cccATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGgaga                                                                                                   <  1:2762775/150‑1 (MQ=255)
                                                  cccATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGgaga                                                                                                   <  1:1180760/150‑1 (MQ=255)
                                                  cccATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGgaga                                                                                                   <  1:2276589/150‑1 (MQ=255)
                                                   ccctGGCGACAGGGACTGCGGCATCGGCGTTGGCGCCCCCTCCCCCTCCCCCGCGCCCGCCGGGAGAGGGGACCGGGCGGCGCGGTATGTACGCCCTCCCTCTGAACCGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGagag                                                                                                  <  2:43643/147‑1 (MQ=255)
                                                    cATGGCGAACGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTACCCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAgg                                                                                                 <  2:2720007/150‑1 (MQ=255)
                                                      tGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGc                                                                                               <  2:3084167/150‑1 (MQ=255)
                                                      tGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCACACCCTAACCCGCACCCTCAGGGGGAGGGGGACCTTGGGGGCGGTAATTAAAACCACCCACAGTAAACACCCCCCCCTCTTTTTTTAACCCTAGCCCCACCCCGGAGAgggg                                                                                               >  2:418727/1‑149 (MQ=255)
                                                       ggCGACAGGGAAGGCGGCATAGTCGTTTGCGCCCCCTCACCCTACCCCTCTCCCTCCGGGAGAGGGGACCGTTCGGCGCTGTCTGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCTCCCCGGAGGGGCGAgggag                                                                                              <  1:2586506/150‑3 (MQ=255)
                                                          gAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCggg                                                                                           <  2:543344/150‑1 (MQ=255)
                                                          gAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCACCCGAACCCACAACCACAGGGGGGGGGGGCCGGTCGGGGCGGTATGGTCTCCCACACACAGGAAACAAACCCCCCACTTTTTTTGTACCACGCCCCACCCGGGAGGGGGCCcgg                                                                                           >  2:3398921/1‑150 (MQ=255)
                                                          gAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCACACCCAAAACCAAAACCCCCCGGAGAGGGGCCCCTTCCGCGCGGGATGTAAGCCCCCACACAGAAAACACACCCCACAAATTTTTTAAACCCACCCCCACCCGGGCGCGGGCCggg                                                                                           >  2:457144/1‑150 (MQ=255)
                                                          gAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCGCACCCGAAACCACTACCACAGGGCGAGGGGACCCTTCGGCGGTGTATGTTCTCCCGCACAATGTAAAGACACCCCACTAGTTTTTTCTCCCACGCCCCACCCGGGGGGGGGCCggg                                                                                           >  2:1443893/1‑150 (MQ=255)
                                                          gAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCACACCCGAACCCACTCCCGAAGGGGGAGGGGGCCCTTCGGCGCGGGATGGACTACCCAAATCTGAAACCACACCCCACACGTTTTTTCTCCCACGCCCCACCCGGGCGAGGGGCggg                                                                                           >  2:3067525/1‑150 (MQ=255)
                                                            aaGGGACTGCGTCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTTAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGt                                                                                         <  1:11987/150‑1 (MQ=255)
                                                                gAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgg                                                                                     <  2:2916835/150‑1 (MQ=255)
                                                                      ggCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCCGGGAGAGGGGACCGTTCGGCGCTGTATGTACGCCCTCCCTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaaa                                                                               <  2:3979909/150‑1 (MQ=255)
                                                                      ggCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaaa                                                                               <  2:1174254/150‑1 (MQ=255)
                                                                      ggCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaaa                                                                               <  2:3221185/150‑1 (MQ=255)
                                                                      ggCATAGTCGTTTGCGCCCCCACACCCTAACCCACTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaaa                                                                               <  2:3532999/150‑1 (MQ=255)
                                                                                tttGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGACACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGCGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTg                                                                     >  2:3145569/1‑150 (MQ=255)
                                                                                         ccTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCCCTGCTCACCAAt                                                            >  2:1549625/1‑150 (MQ=255)
                                                                                               ccTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGGGCGGGCCGGGGGGAGGGGAACAGGCCGCGCTGCGCACCAATTGTcgg                                                      >  1:2345347/1‑148 (MQ=255)
                                                                                               ccTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTg                                                      >  2:4109785/1‑150 (MQ=255)
                                                                                                cTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCACTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGc                                                     >  1:3220992/1‑150 (MQ=255)
                                                                                                    cccTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCTCCCCTCCGGGGAGAGGGCCGGGGTGCGGGGACAAGGCCGAACGGCTAACAAATGTTATGCGATc                                                 >  2:2766991/1‑150 (MQ=255)
                                                                                                     ccTCTCCCTCGGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACGCTGAAACGACACCGCACTCTGTTTTTCTCCCTCGCCCCTCCGGGGCGCGGGCCGGGGTGAGGGGACAAGGCCGCACTGCTCACCCCTTTTCTGCTCTcc                                                >  2:2221833/1‑150 (MQ=255)
                                                                                                     ccTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATcc                                                >  2:2539472/1‑150 (MQ=255)
                                                                                                     ccTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATcc                                                >  1:977312/1‑150 (MQ=255)
                                                                                                         tcCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGcc                                            >  1:1053983/1‑150 (MQ=255)
                                                                                                           ccTCAGGGAGAGGGGACAGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGcccc                                          >  2:947859/1‑150 (MQ=255)
                                                                                                             tCAGGGAGAGGGGACCGTTCGGCGCTGCATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCg                                        >  1:1605121/1‑150 (MQ=255)
                                                                                                               aGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCa                                      >  2:60111/1‑150 (MQ=255)
                                                                                                                       gggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCATTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAACGCCGCACTCCACACCAATTTTCTGCTATCCTGCCACCGCACTAAgggg                              >  2:746959/1‑150 (MQ=255)
                                                                                                                         gACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGc                            >  1:698655/1‑150 (MQ=255)
                                                                                                                          aCCGTTCGGCGTTTTATTTACTCCCAAACTCTCAAACCACACAAAAATATTTTTTCCCCCCCCGCCCCTCGGGGTAGCGGGCGGGGGTGCGGGCAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCa                           <  2:4241130/150‑1 (MQ=255)
                                                                                                                          aCCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCa                           >  1:1309538/1‑150 (MQ=255)
                                                                                                                          aCCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCa                           >  1:2936693/1‑150 (MQ=255)
                                                                                                                            cGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGt                         >  1:4212470/1‑150 (MQ=255)
                                                                                                                              ttCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTAt                       >  1:1577753/1‑150 (MQ=255)
                                                                                                                              ttCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGGGAGGGGACAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTCAGGGGGCAGTAt                       >  2:2987073/1‑150 (MQ=255)
                                                                                                                                gggCGCTTTTTGCACACCCACACGCTAAACCGACACCGCTCTCTTTTTTTCACCCCCGCCCCGGGGGGGCGGGGGGGGGGGTGCGGGCCCAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGc                     <  2:1759551/149‑1 (MQ=255)
                                                                                                                                cGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGc                     >  2:1937164/1‑150 (MQ=255)
                                                                                                                                cGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATAATCTGCTATCCTGCCCCCGCACTAAGGGGGCCGTATGc                     >  2:3615765/1‑150 (MQ=255)
                                                                                                                                cGGCGATGTATGTACTCCCTCACTCTGAAACGAAAACGAACTCTTTTATTCTCCCTCGCCACACCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACAGCTCACCAATTTTCTGATATCCTCACCCCGCACTAAGGGGGCAGTATGc                     >  2:2632090/1‑150 (MQ=255)
                                                                                                                                 ggCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCt                    >  2:3620082/1‑150 (MQ=255)
                                                                                                                                    gcTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaaa                 >  1:611296/1‑150 (MQ=255)
                                                                                                                                     cTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCCGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAAc                >  2:198202/1‑150 (MQ=255)
                                                                                                                                            taCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTAc         >  1:269457/1‑150 (MQ=255)
                                                                                                                                             aCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACt        >  2:582851/1‑150 (MQ=255)
                                                                                                                                             aCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACt        >  1:440446/1‑150 (MQ=255)
                                                                                                                                               tCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGTGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTCAGGGGGCAGTATGCTAAACAGCGTACttt      >  2:3626646/1‑150 (MQ=255)
                                                                                                                                               tCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTACACATCGTACttt      >  2:3367386/1‑150 (MQ=255)
                                                                                                                                               tCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACttt      >  1:4215390/1‑150 (MQ=255)
                                                                                                                                               tCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACttt      >  1:241276/1‑150 (MQ=255)
                                                                                                                                               tCCCTCACTCTGAAACGACACCGCACTCTATTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACttt      >  1:1785825/1‑150 (MQ=255)
                                                                                                                                                cccACACTCCCACACCACACCCCATTCTTTTTGTCCCCCCCGCCCGGGGGGGGTGGGGGGGGGGGTGTGGCCACAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTa     <  1:3897537/150‑1 (MQ=255)
                                                                                                                                                 ccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAACAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTAc    >  1:1256997/1‑150 (MQ=255)
                                                                                                                                                 ccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTCTGCTCAACATCGTACTTTAc    >  1:2932471/1‑150 (MQ=255)
                                                                                                                                                 ccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTAc    >  1:755487/1‑150 (MQ=255)
                                                                                                                                                 ccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTAc    >  1:209210/1‑150 (MQ=255)
                                                                                                                                                 ccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTAc    >  1:100167/1‑150 (MQ=255)
                                                                                                                                                  cTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACg   >  1:1956990/1‑150 (MQ=255)
                                                                                                                                                  cTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACg   >  1:984164/1‑150 (MQ=255)
                                                                                                                                                   tCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGa  >  2:1666797/1‑150 (MQ=255)
                                                                                                                                                   tCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCACCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCACCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGa  >  1:3120209/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGA  >  NC_000913/3780941‑3781237

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: