Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,018,983 A→T 7.4% intergenic (+46/+30) matP → / ← ompA Ter macrodomain organizer matS‑binding protein/outer membrane protein A (3a;II*;G;d)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,018,9830AT 7.4% 809.2 / 35.7 256intergenic (+46/+30)matP/ompATer macrodomain organizer matS‑binding protein/outer membrane protein A (3a;II*;G;d)
Reads supporting (aligned to +/- strand):  ref base A (131/101);  new base T (15/4);  total (150/106)
Fisher's exact test for biased strand distribution p-value = 8.81e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 268 total aligned reads displayed.

AAACGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCC  >  NC_000913/1018834‑1019132
                                                                                                                                                     |                                                                                                                                                     
aanCGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGGTCAATTTACTAAAGGCAAAAAAAACCCCgca                                                                                                                                                       >  1:2703121/1‑150 (MQ=255)
aaaCGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCgca                                                                                                                                                       >  1:641984/1‑150 (MQ=255)
    gATTGTTCAGCTGATTGAAGATTCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATAGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGTCAAAAAAAACCCAgcagcac                                                                                                                                                   >  1:3523244/1‑149 (MQ=255)
       tGTTCAGCTGTTTTAAGATGCGTAAAACAAAGAGAAATACGCTAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCTTTGCTGGATAAGGAATAACTGACGAAAGTCTTTTAAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCgg                                                                                                                                                <  1:946504/150‑1 (MQ=255)
       tGTTCAGCTGAATGAAGATGCTGAAAACAAACAGAAATACGCGAATAAAATGAATTCTCTGAAGCAGGATCTGCAGGCATTGCTGCGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCgg                                                                                                                                                <  2:3883193/150‑1 (MQ=255)
        gTTCAGCTGATTGAATATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCggg                                                                                                                                               <  1:3351619/150‑1 (MQ=255)
        gTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCTAATAAATTTTCTTCTCTTAAGCAGGATCTGCAGGCATTGCTGGTTAAGTAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCggg                                                                                                                                               <  1:2936602/150‑1 (MQ=255)
        gTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCggg                                                                                                                                               <  1:2575228/150‑1 (MQ=255)
        gTTCAGCTGATTGAAGATGCGGAAAAAAAAGAGAAATACGCTAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGTAATAACTGACTAAAGTCAGTTAAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCggg                                                                                                                                               <  1:2166333/150‑1 (MQ=255)
           cAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTTGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAGAAAAACCCCGCTGCTGCGGGGtt                                                                                                                                            <  2:3220288/150‑1 (MQ=255)
             gCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCTATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGtttt                                                                                                                                          >  1:593608/1‑150 (MQ=255)
              cTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGttttt                                                                                                                                         >  2:4124328/1‑150 (MQ=255)
              cTGATTGAAGAGGCGGAAAACAAAGAGAAATACGCGACTAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCAGTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCTGCTGCGGGGttttt                                                                                                                                         <  2:1553453/150‑1 (MQ=255)
                  ttGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTAc                                                                                                                                     <  2:2333133/150‑1 (MQ=255)
                   ttaaGATGCGTAAAACAAAGATAAATACGCTAATAAAATTTCTTCTCTTAAGCAGGATCTGCAGGAATTGCTGGGTAAGGAATAACTGACGAAAGTCTTTTCAATTTACTAACGGAAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTAcc                                                                                                                                    <  1:2114373/148‑1 (MQ=255)
                            ggAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCCGCAGCGGGGTTTTTCTACCAGACGAGaa                                                                                                                           >  2:3460817/1‑150 (MQ=255)
                                 aCAAAGAGAATTCCGCGAATAACATGTCTTCACTGAAGCAGGAGCTGCAGGCATTGCTGGTTAAGGAATAACTGACGAATGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTaa                                                                                                                      <  1:121925/150‑1 (MQ=255)
                                 aCAAAGAGAAATCCGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTaa                                                                                                                      <  1:3974039/150‑1 (MQ=255)
                                 aCAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCTGGCATTGCTGGGTAAGTAATAACTGACGAAAGTAAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTaa                                                                                                                      <  2:378276/150‑1 (MQ=255)
                                 aCAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCCGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAGAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTaa                                                                                                                      <  1:778024/150‑1 (MQ=255)
                                 aCAAAGACAAAAACCCGAAAAAAATTTCTACTCTTAAGCAGGATCTGCAGGCATTGCTGGGTAAGCAATAACTGACGAAAGACAGTTCAAATTACTAAAGCCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTaa                                                                                                                      <  1:3397812/150‑1 (MQ=255)
                                   aaaGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGc                                                                                                                    >  1:502354/1‑150 (MQ=255)
                                   aaaGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTGTCTACCAGACGAGAACTTAAGc                                                                                                                    >  2:4060428/1‑150 (MQ=255)
                                       agaAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGc                                                                                                                <  2:901862/150‑1 (MQ=255)
                                           aTACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCt                                                                                                            >  1:1839792/1‑150 (MQ=255)
                                           aTACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCt                                                                                                            >  1:11566/1‑150 (MQ=255)
                                           aTACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGAAGCAGCGGGGTTTTTCTACAAGACGAGAACTTAAGCCTGCGGCt                                                                                                            >  2:799902/1‑150 (MQ=255)
                                               gcgAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGt                                                                                                        >  2:1478005/1‑150 (MQ=255)
                                                cgAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGtt                                                                                                       >  1:2332358/1‑150 (MQ=255)
                                                 gAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAATAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTa                                                                                                      >  1:136084/1‑150 (MQ=255)
                                                 gAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTa                                                                                                      >  2:3659403/1‑150 (MQ=255)
                                                 gAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTa                                                                                                      >  1:1943823/1‑150 (MQ=255)
                                                   aTAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACa                                                                                                    <  1:2461395/150‑1 (MQ=255)
                                                   aTAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACa                                                                                                    <  1:2828489/150‑1 (MQ=255)
                                                   aTAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACa                                                                                                    <  2:2585720/150‑1 (MQ=255)
                                                   aTAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACa                                                                                                    >  2:3127425/1‑150 (MQ=255)
                                                   aTAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACa                                                                                                    <  2:3355225/150‑1 (MQ=255)
                                                   aTAAAATCTCTACTCTCAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTCACGAAACTCATTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACa                                                                                                    <  2:1756806/150‑1 (MQ=255)
                                                             ttctCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGa                                                                                          <  1:3207539/150‑1 (MQ=255)
                                                             ttctCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGa                                                                                          >  2:1450882/1‑150 (MQ=255)
                                                                tctGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATAc                                                                                       >  2:2056425/1‑150 (MQ=255)
                                                                tctGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATAc                                                                                       >  1:1402662/1‑150 (MQ=255)
                                                                tctGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTGTGATAc                                                                                       >  1:2491464/1‑150 (MQ=255)
                                                                 ctGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATAcc                                                                                      <  1:3846206/150‑1 (MQ=255)
                                                                      gCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTaa                                                                                 >  2:2238426/1‑150 (MQ=255)
                                                                      gCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTCTGATACCTTTaa                                                                                 >  1:160365/1‑150 (MQ=255)
                                                                               gCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATCCATTTAACTTCGAtct                                                                        >  1:1665679/1‑150 (MQ=255)
                                                                               gCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGAtct                                                                        <  2:3913820/150‑1 (MQ=255)
                                                                                  ggCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATctctc                                                                     >  2:2158934/1‑149 (MQ=255)
                                                                                   gCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTAc                                                                    <  2:1640401/150‑1 (MQ=255)
                                                                                   gCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTGTGAGACCTTTAACTTCGAGCGCTAc                                                                    >  2:2991281/1‑150 (MQ=255)
                                                                                   gCATGGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCACTTTACTAAAGGCAAAAAAAACCCCGCCGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCGAc                                                                    <  2:3691579/150‑1 (MQ=255)
                                                                                    cATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTAcg                                                                   >  1:1874812/1‑150 (MQ=255)
                                                                                    cATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTAcg                                                                   <  1:2960927/150‑1 (MQ=255)
                                                                                     aTTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAACCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTAcgc                                                                  <  1:3314498/150‑1 (MQ=255)
                                                                                       tGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTATCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGcga                                                                >  2:2490178/1‑150 (MQ=255)
                                                                                         cTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTATCTTCGATCTCTACGcgacg                                                              >  2:3307413/1‑150 (MQ=255)
                                                                                         cTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGcgacg                                                              >  2:3394240/1‑150 (MQ=255)
                                                                                          tGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGcgacga                                                             >  1:3947139/1‑150 (MQ=255)
                                                                                          tGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGcgacga                                                             >  1:3067588/1‑150 (MQ=255)
                                                                                             gTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATcc                                                          <  2:3108669/150‑1 (MQ=255)
                                                                                                aGGAATAACTGACGACAGTCAGTTCACTTTACTACCGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACACCGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGa                                                       <  2:2559967/150‑1 (MQ=255)
                                                                                                  cactaccTGCCGAAAGTCAGTCCAATCTACAAACGCCAAAAAAAACCCCGCAGCCGCGGGGATTTTCTACCAGACGAGAACTCAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGc                                                     <  2:2370933/144‑1 (MQ=255)
                                                                                                     tAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAg                                                  <  1:1775926/150‑1 (MQ=255)
                                                                                                       aCTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGc                                                >  2:1965309/1‑150 (MQ=255)
                                                                                                         tGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAg                                              >  1:1123056/1‑150 (MQ=255)
                                                                                                         tGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTGTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAg                                              >  1:1448736/1‑150 (MQ=255)
                                                                                                           aCGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTc                                            >  1:719592/1‑150 (MQ=255)
                                                                                                             gAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGa                                          >  1:2743713/1‑150 (MQ=255)
                                                                                                                aGTCAGTTCAATTTACGAACGGCAAACAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCa                                       <  2:2703121/150‑1 (MQ=255)
                                                                                                                   cAGTTCAATTTACTAAACGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTg                                    <  1:3156368/150‑1 (MQ=255)
                                                                                                                    aGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTgc                                   >  1:2209789/1‑150 (MQ=255)
                                                                                                                      ttCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACGTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCCGGCCGTCGAGCCGTgcag                                 >  2:1098280/1‑150 (MQ=255)
                                                                                                                             tACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTg                          <  1:1470691/150‑1 (MQ=255)
                                                                                                                                tAAAGGCAAAAAAATCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGttt                       >  1:3084319/1‑150 (MQ=255)
                                                                                                                                tAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGttt                       >  1:2075281/1‑150 (MQ=255)
                                                                                                                                tAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGttt                       >  1:2766492/1‑150 (MQ=255)
                                                                                                                                  aaGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCa                     <  2:1715020/150‑1 (MQ=255)
                                                                                                                                  aaGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCa                     >  1:3674969/1‑150 (MQ=255)
                                                                                                                                   aGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCAc                    <  1:3935014/150‑1 (MQ=255)
                                                                                                                                      cAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGtt                 <  1:1066563/150‑1 (MQ=255)
                                                                                                                                      cAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGtt                 >  1:829884/1‑150 (MQ=255)
                                                                                                                                        aaaaaaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGt               >  2:575672/1‑150 (MQ=255)
                                                                                                                                        aaaaaaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGt               <  1:3507968/150‑1 (MQ=255)
                                                                                                                                        aaaaaaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGt               <  2:3623893/150‑1 (MQ=255)
                                                                                                                                        aaaaaaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCGGCGGCTGTGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTCCGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGt               >  2:2620397/1‑150 (MQ=255)
                                                                                                                                         gaaaaaCCCCGCTGCTGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTc              >  2:1942769/2‑150 (MQ=255)
                                                                                                                                         gaaaaaCCCCGCTGCTGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTc              >  1:2369416/2‑150 (MQ=255)
                                                                                                                                         aaaaaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTATCTTCGATCTCTACGCGACGATCCGGAGCCAGGCCGTCGATCAGTGCTGCACGCTGTTTCTCGTTGTc              >  2:3709319/1‑150 (MQ=255)
                                                                                                                                         aaaaaaCCCCGCAGCAGCGGGGTTTTTCTAACAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGATGTTTCACGTTGTc              >  1:3877123/1‑150 (MQ=255)
                                                                                                                                          aaaaaCCCCGCTGCTGCGGGGTTTTTTTACCAGCCGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCCGCACGCTGTTTCACGTTGTca             >  1:2380391/1‑150 (MQ=255)
                                                                                                                                            aaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGTTCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTcaca           >  2:1045828/1‑150 (MQ=255)
                                                                                                                                             aaCCCCGCAGCAGCGGGTTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTATTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCTCGCTGTTTCACGTTGTCACAg          <  2:1185687/150‑1 (MQ=255)
                                                                                                                                             aaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAg          <  2:757659/150‑1 (MQ=255)
                                                                                                                                                cccGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCCGCACGCTGTTTCACGTTGGCACAGgtg       >  1:3624785/1‑150 (MQ=255)
                                                                                                                                                  cGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGtt     >  1:2034594/1‑150 (MQ=255)
                                                                                                                                                  cGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGATGTCACAGGTGtt     >  2:44181/1‑150 (MQ=255)
                                                                                                                                                    cagcagCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGACTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGc   <  2:136102/150‑1 (MQ=255)
                                                                                                                                                    cagcagCGGGGTTTTTCTACCAGACGAGAACATAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATACGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTAACGTTGTCACAGGTGTTGc   >  1:1527750/1‑150 (MQ=255)
                                                                                                                                                     agcagcGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGcc  <  2:3320425/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAACGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCC  >  NC_000913/1018834‑1019132

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: