Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 518,334 T→G 7.2% intergenic (‑55/+6) ybbN ← / ← ybbO DnaK co‑chaperone, thioredoxin‑like protein/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913518,3340TG 7.2% 773.1 / 5.6 275intergenic (‑55/+6)ybbN/ybbODnaK co‑chaperone, thioredoxin‑like protein/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (116/129);  new base G (25/2);  total (143/132)
Fisher's exact test for biased strand distribution p-value = 2.89e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 285 total aligned reads displayed.

ACCAAAAATAGAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCTTC  >  NC_000913/518186‑518483
                                                                                                                                                    |                                                                                                                                                     
  cAAAAATAGAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCa                                                                                                                                                    <  1:1414950/150‑1 (MQ=255)
    aaaaTAGAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCCGTTGTCGATATTTTCTTTGAAATGGGGGCTTAAGCGCGCGCTTCAAc                                                                                                                                                  >  1:2495473/1‑150 (MQ=255)
      aaTAGAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTc                                                                                                                                                >  2:383504/1‑150 (MQ=255)
       aTAGAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCa                                                                                                                                               >  2:2965965/1‑150 (MQ=255)
        tAGAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCAc                                                                                                                                              <  1:3559511/150‑1 (MQ=255)
        tAGAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCAc                                                                                                                                              <  1:159936/150‑1 (MQ=255)
        tAGAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCAc                                                                                                                                              <  1:498074/150‑1 (MQ=255)
           aaCAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTTACATGGGGGCTTTAGCGCGCGCTTCAACGCAcccc                                                                                                                                           >  1:1029256/1‑150 (MQ=255)
           aaCAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGGTGTCGATATTTTCTTTGAAAAGGGGGCTTAATCGCGCGCGTGAACTCAcccc                                                                                                                                           >  1:1782302/1‑150 (MQ=255)
            aCAGCACCGGAGTGGTCATCGACTGTTCAATAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCt                                                                                                                                          <  1:3818294/150‑1 (MQ=255)
            aCAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCt                                                                                                                                          <  2:3365421/150‑1 (MQ=255)
              aGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTGGATATTTTCTTTGACATGGGGGCTGAAGCGCGCGCGTCGACTCACCCCGGc                                                                                                                                        >  1:906600/1‑150 (MQ=255)
               gCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCa                                                                                                                                       >  1:2353578/1‑150 (MQ=255)
                  ccGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAata                                                                                                                                    <  1:4068301/150‑1 (MQ=255)
                    ggAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATAtt                                                                                                                                  <  1:1675129/150‑1 (MQ=255)
                     gAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATAttt                                                                                                                                 <  2:3465251/150‑1 (MQ=255)
                     gAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATAttt                                                                                                                                 <  2:3569959/150‑1 (MQ=255)
                      aGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATAtttt                                                                                                                                <  1:660779/150‑1 (MQ=255)
                      aGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATAtttt                                                                                                                                <  1:936748/150‑1 (MQ=255)
                      aGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATAtttt                                                                                                                                <  2:609979/150‑1 (MQ=255)
                       gTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTg                                                                                                                               <  1:1953240/150‑1 (MQ=255)
                        tGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGCTATTTTCTTTGACATGGGGGCTTAATCGCGCGCTTCAACTCACCCCGGCAATATTTTGt                                                                                                                              >  1:1656434/1‑150 (MQ=255)
                         ggTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTc                                                                                                                             <  1:3963937/150‑1 (MQ=255)
                           tCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCGGGGTGTTGGGGGTATTTTCTTTTAAATTGGGGGGTTAGCGCGCGCGTAAAAGAAACCCAGCAAAAttttttgca                                                                                                                           >  1:1831332/1‑145 (MQ=255)
                           tCATCGACTGTTCAAGAACCTGTTGCAGGATAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCGCTGGTGTCGATATATTATGTGACAAGGGGGCGTAAGCGCGCGCGTCAACACACCCCGGCAAAAATTTGTCCa                                                                                                                           >  2:2829016/1‑150 (MQ=255)
                            cATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCAt                                                                                                                          <  2:1003777/150‑1 (MQ=255)
                             aTCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATc                                                                                                                         <  1:1287274/150‑1 (MQ=255)
                              tCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGGTGTCGAGATTTTCTTTGAAATGGGGGGGTAAGCGCGCGCGTCAACACACCCCGGCAAAAATTTTTCCATCa                                                                                                                        >  1:1175672/1‑150 (MQ=255)
                               cGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCAc                                                                                                                       <  1:1053358/150‑1 (MQ=255)
                                gACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTAACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCAcg                                                                                                                      >  2:3444894/1‑150 (MQ=255)
                                   tGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACgcgc                                                                                                                   <  1:2780728/150‑1 (MQ=255)
                                     ttCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAAATCACCCCTGCAATATTTTGTCCATCACGCGccc                                                                                                                 >  2:566230/1‑150 (MQ=255)
                                     ttCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAAATCACCCCAGCAATATTTTGTCCATCACGCGccc                                                                                                                 >  2:659926/1‑150 (MQ=255)
                                        aaGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAAAATTTTGTCCATCACGCGCCCCgg                                                                                                              >  2:256829/1‑150 (MQ=255)
                                        aaGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACAAGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCgg                                                                                                              >  2:1705901/1‑150 (MQ=255)
                                        aaGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGGTGTCGATATTTTATTTTAAATGGGGGCTTAATCCCCCGCTTAAACACACCCCCGCAATAATTTTTCCATCACCCCCCCCgg                                                                                                              >  2:3798215/1‑150 (MQ=255)
                                         aGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTTTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCATCAACTCACCCCTGCAATAATTTGTCCATCACGCGCCCCGgc                                                                                                             >  2:2720886/1‑150 (MQ=255)
                                         aGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTGGTTATTTTCTTTTACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGgc                                                                                                             >  2:1727303/1‑150 (MQ=255)
                                         aGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCAGTTGTCGATATTTTCTTTGGCATGGGGGCGTAAGCGCGCGCTTCAACACACCCCAGCAAAATTTTGTCCAACACCCGCCCCGgc                                                                                                             >  2:3941840/1‑150 (MQ=255)
                                          gAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTATTTGAAATTGGGGCTTAAGCGCGCGCTTCAACTCACCCCGGCAATATTTTGTCCATCACGCGCCCCGgca                                                                                                            >  2:4039386/1‑150 (MQ=255)
                                           aaCCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGgcag                                                                                                           <  1:3197491/150‑1 (MQ=255)
                                           aaCCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGgcag                                                                                                           <  1:2149312/150‑1 (MQ=255)
                                                 ttGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGgcg                                                                                                     <  1:3051230/150‑1 (MQ=255)
                                                  tGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGgcgc                                                                                                    <  2:1854091/150‑1 (MQ=255)
                                                       gTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAg                                                                                               >  2:3706411/1‑150 (MQ=255)
                                                               cGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTa                                                                                       <  1:3671623/150‑1 (MQ=255)
                                                                       ttGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAg                                                                               <  2:2689121/150‑1 (MQ=255)
                                                                          aCAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTAGATATTTTCTTTGTCATGGGGGCTTAATCGCGCGCTTCAAAACACCCCTGCAATATTTTGTCCATCACCCGCCCCGGCAGCAGGCGCTTAAGCACCAATACCGCCCCGGTc                                                                            >  2:1920281/1‑150 (MQ=255)
                                                                           cAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCa                                                                           <  2:2316590/150‑1 (MQ=255)
                                                                                       cGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTTACAAGGGGGCTTAAGCGCGCGCTTCAACTCACCCCGGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCGTAAGCACCCATACCGCCCAGGTCACCCCCGGCACCg                                                               >  2:3901217/1‑150 (MQ=255)
                                                                                       cGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCAACCCTGCACTATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCg                                                               >  2:4064995/1‑150 (MQ=255)
                                                                                        ggACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACACAACCCTGCAATATTTTGTACATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCAATACCGCCCAGGTCACCAGCGTCACCgg                                                              >  2:1289635/1‑150 (MQ=255)
                                                                                        ggACATGGAGTCGCTCTCTGTTGTCGATATTTTCATTGACATGGGGGCTTAATCGCGCGCTTCAACTCAACCCAGCAAAAATTTGTACAAAACCCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCATCGTCACCgg                                                              >  2:1591552/1‑150 (MQ=255)
                                                                                               gAGTCGCTGTCTGTTGTCGATATTTTATTTGACATTGGGGCTTAAGCGCGCGCTTCAACTCACCCCAGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACACAAACCCCCCAGGTCACCCGCGTCACCGGATAgcgc                                                       >  2:3224164/1‑150 (MQ=255)
                                                                                                 gTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCAt                                                     <  1:208140/150‑1 (MQ=255)
                                                                                                    gCTCTCTGTTGTCGATATTTTCTTTGACATTGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATATTGTCCATCACCCGCCCCGGCAGCAGGCGCATAACCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATcat                                                  >  1:202271/1‑148 (MQ=255)
                                                                                                    gCTCTCTGTTGTCGATATTTTATTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCAGCAATATTTTGTCCATCACGCGCCCCCGCAGCAGGCGCTTAAGCACCAATACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCtt                                                  >  1:549665/1‑150 (MQ=255)
                                                                                                     ctctctGTTGTCGATATTTTCTTTGACATGGGGGCTTTAGCGCGCGCTTCAACTCACCCCTGCAATATATTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTc                                                 <  1:747554/150‑1 (MQ=255)
                                                                                                     ctctctGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTc                                                 <  1:1562038/150‑1 (MQ=255)
                                                                                                     ctctctGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTc                                                 >  2:2915830/1‑150 (MQ=255)
                                                                                                     ctctctGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTc                                                 <  2:1831332/150‑1 (MQ=255)
                                                                                                     ctctctGTTGTCGATATTTTATTTGACATGGGGGCTTAAGCGCGCGCTTCAAATCACCCCTGCAATATTTTGTACATAACGCGCCCCGGCAGCAGGCGCTTAAGCACCAATACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTc                                                 >  1:3869876/1‑150 (MQ=255)
                                                                                                       ctctGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCgg                                               <  1:2042873/150‑1 (MQ=255)
                                                                                                             tgtCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTtctc                                         <  1:2521083/150‑1 (MQ=255)
                                                                                                              gtCGATATTTTCTTTGACATGGGGGCTTAATCGCGCGCGTCAACACACCCCTGCAATATTTTTTACATCAAGCGCCCCGGCAGCAGGCGCTTAAGCACCAATACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCg                                        >  1:1491793/1‑150 (MQ=255)
                                                                                                              gtCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCCGCAGCAGGCGCTTAAGCACCAAAACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCg                                        >  1:4078826/1‑150 (MQ=255)
                                                                                                                 gATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCta                                     >  2:1089682/1‑150 (MQ=255)
                                                                                                                 gATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCta                                     >  2:1748316/1‑150 (MQ=255)
                                                                                                                  ataggttCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCtaa                                    <  1:3480194/145‑1 (MQ=255)
                                                                                                                   tatTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCtaat                                   <  2:214583/150‑1 (MQ=255)
                                                                                                                   tatTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCtaat                                   >  1:3273437/1‑150 (MQ=255)
                                                                                                                   tatTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCtaat                                   >  1:2104922/1‑150 (MQ=255)
                                                                                                                     ttttCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGGAATATTTTGTCCATCTCGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCAGCAGCGTCTCCGGATAGCGCATGTTCGGCTTCTCGCtaataa                                 >  2:1202180/1‑150 (MQ=255)
                                                                                                                     ttttCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCtaataa                                 >  2:858496/1‑150 (MQ=255)
                                                                                                                     ttttCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCtaataa                                 >  2:817085/1‑150 (MQ=255)
                                                                                                                     ttttCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCtaataa                                 >  2:1460790/1‑150 (MQ=255)
                                                                                                                     ttttCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCtaataa                                 >  2:1036730/1‑150 (MQ=255)
                                                                                                                       ttCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATaaaa                               >  2:2063913/1‑150 (MQ=255)
                                                                                                                         cTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGc                             <  2:3089871/150‑1 (MQ=255)
                                                                                                                           ttGACATGGGGGCTTAAGCGCGCGGTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCAt                           >  1:2258066/1‑150 (MQ=255)
                                                                                                                           ttGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCAt                           >  2:1653185/1‑150 (MQ=255)
                                                                                                                           ttGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCAt                           >  1:3639090/1‑150 (MQ=255)
                                                                                                                            tGACATGGGGTCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATg                          <  2:11132/150‑1 (MQ=255)
                                                                                                                               cATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCg                       >  2:2483044/1‑150 (MQ=255)
                                                                                                                                 tGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTa                     >  1:3029601/1‑150 (MQ=255)
                                                                                                                                 tGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCAATACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTa                     >  1:3719228/1‑150 (MQ=255)
                                                                                                                                 tGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATAATTTGTCCATCACGCGCCCCTGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTa                     <  2:685699/150‑1 (MQ=255)
                                                                                                                                  gggggCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTAc                    <  1:2730616/150‑1 (MQ=255)
                                                                                                                                   ggggCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACt                   >  1:3144509/1‑150 (MQ=255)
                                                                                                                                   ggggCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCACTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACt                   >  2:1681268/1‑150 (MQ=255)
                                                                                                                                      gCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTg                >  1:3749757/1‑150 (MQ=255)
                                                                                                                                        ttAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTGTc              <  1:389874/150‑1 (MQ=255)
                                                                                                                                         tAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTGTcc             <  2:1287229/150‑1 (MQ=255)
                                                                                                                                          aaGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTGTcca            >  2:435764/1‑150 (MQ=255)
                                                                                                                                                 cgcTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGc     >  1:418573/1‑150 (MQ=255)
                                                                                                                                                  gcTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCt    >  2:2990888/1‑150 (MQ=255)
                                                                                                                                                  gcTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCt    >  2:919141/1‑150 (MQ=255)
                                                                                                                                                    ttCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCTTc  >  2:2096919/1‑150 (MQ=255)
                                                                                                                                                    ttCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCTTc  <  2:2011690/150‑1 (MQ=255)
                                                                                                                                                    ttCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACTTAGTCCACCACCGCTTc  >  2:2168010/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ACCAAAAATAGAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGATAGCGCATCTTCGGCTTCTCGCTAATAAAAGCATGGCGTACTTTGTCCACCACCGCTTC  >  NC_000913/518186‑518483

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: