Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,018,986 A→T 5.7% intergenic (+49/+27) 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,9860AT 5.7% 831.4 / 30.1 237intergenic (+49/+27)matP/ompATer macrodomain organizer matS‑binding protein/outer membrane protein A (3a;II*;G;d)
Reads supporting (aligned to +/- strand):  ref base A (111/111);  new base T (12/2);  total (123/114)
Fisher's exact test for biased strand distribution p-value = 1.15e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 254 total aligned reads displayed.

GATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGT  >  NC_000913/1018838‑1019135
                                                                                                                                                     |                                                                                                                                                     
gATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGcagcag                                                                                                                                                      <  1:2552244/150‑1 (MQ=255)
   tGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCgg                                                                                                                                                   <  2:1624455/150‑1 (MQ=255)
     ttCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCgggg                                                                                                                                                 >  2:1751651/1‑150 (MQ=255)
       cAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGtt                                                                                                                                               <  1:970487/150‑1 (MQ=255)
       cAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAACCCCGCAGCAGCGGGGttt                                                                                                                                              >  1:2168459/1‑150 (MQ=255)
        aGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGttt                                                                                                                                              >  1:3687091/1‑150 (MQ=255)
          cTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGttttt                                                                                                                                            <  1:3597243/150‑1 (MQ=255)
             aTTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGACGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTa                                                                                                                                         <  2:1486210/150‑1 (MQ=255)
             aTTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTa                                                                                                                                         <  2:3165475/150‑1 (MQ=255)
              ttGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTAc                                                                                                                                        >  1:765188/1‑150 (MQ=255)
               tGACGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTAcc                                                                                                                                       <  2:2833015/150‑1 (MQ=255)
                     tGCGTAAAACAAAGATAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCATTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACga                                                                                                                                 <  1:2313918/150‑1 (MQ=255)
                     tGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACga                                                                                                                                 <  1:2412455/150‑1 (MQ=255)
                        ggAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGaa                                                                                                                              <  2:3098393/150‑1 (MQ=255)
                        ggAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGaa                                                                                                                              >  1:1739334/1‑150 (MQ=255)
                        ccaaaaaaaaCACAAATACCCTAAAAAAATGTCTACACTCAACCAGCATCTCCAGGCATTGCTGGGTAACCAATAACTGACCAAACTCACTACAATTTACTAAAGCCAAAAAAAACCCCGCAGCAGCGGGGATTTTCTACCAGACGAGaa                                                                                                                              <  1:2780095/148‑1 (MQ=255)
                         gAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACgagacc                                                                                                                             >  2:2689730/1‑148 (MQ=255)
                         gAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAAc                                                                                                                             <  2:1655584/150‑1 (MQ=255)
                             aCAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTaa                                                                                                                         <  2:2549147/150‑1 (MQ=255)
                             aCAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTaa                                                                                                                         <  1:2579772/150‑1 (MQ=255)
                               aaaGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGc                                                                                                                       >  1:1982643/1‑150 (MQ=255)
                               aaaGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGc                                                                                                                       >  1:3497464/1‑150 (MQ=255)
                               aaaGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGc                                                                                                                       >  1:3033318/1‑150 (MQ=255)
                                    gaAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCg                                                                                                                  >  1:870264/1‑150 (MQ=255)
                                       aTACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCt                                                                                                               <  1:3864626/150‑1 (MQ=255)
                                        tACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGa                                                                                                             <  2:3803742/150‑1 (MQ=255)
                                         aCGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGa                                                                                                             >  1:3098690/1‑150 (MQ=255)
                                           gcgAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGt                                                                                                           >  2:3838388/1‑150 (MQ=255)
                                                 aaaaTGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAAc                                                                                                     <  2:1739334/150‑1 (MQ=255)
                                                       tcttctCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCttt                                                                                               <  2:1664445/150‑1 (MQ=255)
                                                            tctGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATAc                                                                                          <  1:492185/150‑1 (MQ=255)
                                                            tctGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATAc                                                                                          >  1:1676393/1‑150 (MQ=255)
                                                            tctGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATAc                                                                                          >  1:1550049/1‑150 (MQ=255)
                                                            tctGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATAc                                                                                          >  1:2024591/1‑150 (MQ=255)
                                                             ctGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAATTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATAcc                                                                                         <  2:2230990/150‑1 (MQ=255)
                                                             ctGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATAcc                                                                                         >  2:2349109/1‑150 (MQ=255)
                                                               gAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCtt                                                                                       >  1:79907/1‑150 (MQ=255)
                                                                   cAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAAc                                                                                   <  2:3912101/150‑1 (MQ=255)
                                                                   cAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTGACAACGTCTTTGATACCTGTAAc                                                                                   >  2:3684612/1‑150 (MQ=255)
                                                                      gATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTc                                                                                <  1:2450587/150‑1 (MQ=255)
                                                                         cTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGAt                                                                             <  1:2734716/150‑1 (MQ=255)
                                                                          tGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGAtc                                                                            >  1:1369377/1‑150 (MQ=255)
                                                                          tGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGAtc                                                                            >  1:10504/1‑150 (MQ=255)
                                                                          tGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGAtc                                                                            >  1:460537/1‑150 (MQ=255)
                                                                          tGCAGGCATTGCTGGGTAAGGAATAACAGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGAtc                                                                            >  2:3690648/1‑150 (MQ=255)
                                                                             aGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATctct                                                                         >  2:3628214/1‑150 (MQ=255)
                                                                               gCATTGCTGGGTAAGGAATAACTGACGAAACTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTAc                                                                       <  1:1644736/150‑1 (MQ=255)
                                                                                cATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTAcg                                                                      <  1:138359/150‑1 (MQ=255)
                                                                                cATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTAcg                                                                      <  2:746698/150‑1 (MQ=255)
                                                                                     cTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTGACAACGTCTTTGATACCGTGAACTTCGATCTCTACGcgacg                                                                 >  1:1014366/1‑150 (MQ=255)
                                                                                         gTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATcc                                                             <  2:2595973/150‑1 (MQ=255)
                                                                                          tAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCg                                                            >  1:1161131/1‑150 (MQ=255)
                                                                                             ggAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGATTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAg                                                         <  2:224968/150‑1 (MQ=255)
                                                                                              gAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGc                                                        >  1:1460664/1‑150 (MQ=255)
                                                                                               aaTAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGcc                                                       <  2:3497464/150‑1 (MQ=255)
                                                                                                    cTGACGAAAGTCAGTTCAATTTACTACAGGCAACAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCa                                                  <  2:1228914/150‑1 (MQ=255)
                                                                                                     tGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAg                                                 >  2:592703/1‑150 (MQ=255)
                                                                                                       aCGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTc                                               >  1:2275058/1‑150 (MQ=255)
                                                                                                       aCGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAAATTAATCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACTCGACGATCCGGAGCCAGGCAGTc                                               >  2:3545823/1‑150 (MQ=255)
                                                                                                         gAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGa                                             >  1:312642/1‑150 (MQ=255)
                                                                                                               cAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCCGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTg                                       <  1:326914/150‑1 (MQ=255)
                                                                                                                aGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTgc                                      <  1:3943364/150‑1 (MQ=255)
                                                                                                                 ggacAATTTACTAAAGGAAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCTATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTgca                                     <  2:3432986/147‑1 (MQ=255)
                                                                                                                 gTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTgca                                     <  1:760955/150‑1 (MQ=255)
                                                                                                                 gTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTgca                                     <  1:1239288/150‑1 (MQ=255)
                                                                                                                 gTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTgca                                     <  1:230259/150‑1 (MQ=255)
                                                                                                                  ttCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTgcag                                    >  2:1275109/1‑150 (MQ=255)
                                                                                                                  ttCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTgcag                                    >  1:3397608/1‑150 (MQ=255)
                                                                                                                  ttCAATTTACTAAAGGCAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTgcagc                                   >  2:1455029/1‑150 (MQ=255)
                                                                                                                          aCTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGt                            >  2:3381281/1‑150 (MQ=255)
                                                                                                                           cTAAAGGCAAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGt                            >  1:1893593/1‑150 (MQ=255)
                                                                                                                            tAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGttt                          >  2:1815331/1‑150 (MQ=255)
                                                                                                                            tAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGGTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGttg                          >  2:2477690/1‑149 (MQ=255)
                                                                                                                             aaaGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTc                         <  1:3338278/150‑1 (MQ=255)
                                                                                                                              caGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCa                        <  2:1682135/149‑1 (MQ=255)
                                                                                                                              aaGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCa                        <  2:56706/150‑1 (MQ=255)
                                                                                                                                 gCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGACCTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGt                     <  2:860169/150‑1 (MQ=255)
                                                                                                                                  cAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGtt                    <  1:2689730/150‑1 (MQ=255)
                                                                                                                                   gcaaaaaaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGtt                    <  2:2168459/148‑1 (MQ=255)
                                                                                                                                      aaaaaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCGGCACGCTGTGTCGCGTTGTc                 >  2:1092526/1‑150 (MQ=255)
                                                                                                                                       aaaaaCCCCGCAGCAGCGGGGTTTTTCTCCCAGACGAGAACTTAAGCCTGCGGCTGAGGTACCACGTCGTTGATACCTTTACCTTCGATCTCTCCGCGACGATCCGGCGCCAGGCAGTCGATCCGTGCAGCACGCTGTTTCACGGTGTcc                >  2:141464/1‑149 (MQ=255)
                                                                                                                                       aaaaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTca                >  2:3096094/1‑150 (MQ=255)
                                                                                                                                       aaaaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTca                <  2:243438/150‑1 (MQ=255)
                                                                                                                                       aaaaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTca                <  1:2965255/150‑1 (MQ=255)
                                                                                                                                         aaaCCCCGCTGCTGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTcaca              >  1:3946042/1‑150 (MQ=255)
                                                                                                                                         aaaCCCCGCTGCTGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTcaca              >  1:3772550/1‑150 (MQ=255)
                                                                                                                                         aaaCCCCGCTGCTGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTcaca              >  1:1699903/1‑150 (MQ=255)
                                                                                                                                         aaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTcaca              >  1:2886712/1‑150 (MQ=255)
                                                                                                                                           aCCCCGCTGCTGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGGTGTCACAgg            >  1:3876009/1‑150 (MQ=255)
                                                                                                                                           aCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAgg            <  1:601427/150‑1 (MQ=255)
                                                                                                                                             cccGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGgtg          >  2:2282382/1‑150 (MQ=255)
                                                                                                                                             cccGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGgtg          >  2:2017214/1‑150 (MQ=255)
                                                                                                                                              ccGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCCGTGCAGCACGCTGTTTCACGTTGTCACAGgtgt         >  2:729125/1‑150 (MQ=255)
                                                                                                                                              ccGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGgtgt         <  2:1594686/150‑1 (MQ=255)
                                                                                                                                                gcagcaGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTg       >  2:1028474/1‑150 (MQ=255)
                                                                                                                                                 cagcagCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGc      >  1:3163112/1‑150 (MQ=255)
                                                                                                                                                 cagcagCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGc      <  1:2586249/150‑1 (MQ=255)
                                                                                                                                                 cagcagCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGc      <  1:3744007/150‑1 (MQ=255)
                                                                                                                                                   gcagcGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACATGTGTTGCCa    <  1:2345739/150‑1 (MQ=255)
                                                                                                                                                   gcagcGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCa    <  1:2207638/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGT  >  NC_000913/1018838‑1019135

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: