Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,128,264 A→T R44R (CGT→AGA)  metJ ← DNA‑binding transcriptional repressor MetJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,128,2640AT99.9% 5958.3 / NA 1581R44R (CGT→AGA) metJDNA‑binding transcriptional repressor MetJ
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base T (707/871);  total (710/871)
Fisher's exact test for biased strand distribution p-value = 2.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 1594 total aligned reads displayed.

ATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATC  >  NC_000913/4128121‑4128408
                                                                                                                                               |                                                                                                                                                
aTCTCTTTTGCTGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGtctacg                                                                                                                                             >  2:53863/1‑149 (MQ=255)
aTCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGtctacg                                                                                                                                             <  1:1187183/149‑1 (MQ=255)
aTCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGtctacg                                                                                                                                             >  2:1890198/1‑149 (MQ=255)
  ctctTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCTGCTCGCTGTGGGTAGCGTGACGCAGGTTGTTCGCCTGTCTAcgcg                                                                                                                                           >  2:1507549/1‑149 (MQ=255)
   tctTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACgcgt                                                                                                                                          >  2:2941950/1‑149 (MQ=255)
    ctTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACgcgtg                                                                                                                                         <  1:901557/149‑1 (MQ=255)
    ctTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACgcgtg                                                                                                                                         >  1:1808975/1‑149 (MQ=255)
    ctTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACgcgtg                                                                                                                                         <  1:2388923/149‑1 (MQ=255)
    ctTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACgcgtg                                                                                                                                         <  2:2494235/149‑1 (MQ=255)
    ctTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACgcgtg                                                                                                                                         >  1:1090942/1‑149 (MQ=255)
      tttGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACgcgtgcg                                                                                                                                       >  1:1478519/1‑149 (MQ=255)
      tttGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACgcgtgcg                                                                                                                                       >  1:2909878/1‑149 (MQ=255)
        tGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGtt                                                                                                                                     >  2:775610/1‑149 (MQ=255)
        tGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCATGTCTACGCGTGCGtt                                                                                                                                     <  2:14064/149‑1 (MQ=255)
            gCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATc                                                                                                                                 <  1:719300/149‑1 (MQ=255)
            gCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATc                                                                                                                                 <  1:1992759/149‑1 (MQ=255)
              ttCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCgg                                                                                                                               <  1:851735/149‑1 (MQ=255)
              ttCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCgg                                                                                                                               <  2:834198/149‑1 (MQ=255)
                ccGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTg                                                                                                                             <  1:640859/149‑1 (MQ=255)
                   ggATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCTCCTGTCTACGCGTGCGTTCATCGGTGAgg                                                                                                                          >  2:2338473/1‑149 (MQ=255)
                   ggATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAgg                                                                                                                          <  2:2770445/149‑1 (MQ=255)
                   ggATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAgg                                                                                                                          >  1:1129803/1‑149 (MQ=255)
                       ttCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGAttt                                                                                                                      <  2:1311165/149‑1 (MQ=255)
                         cgtcgCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGAttttt                                                                                                                    <  1:1528085/149‑1 (MQ=255)
                           tcgCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTaa                                                                                                                  >  2:1661095/1‑149 (MQ=255)
                            cgCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAac                                                                                                                 >  2:178386/1‑149 (MQ=255)
                               tGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACAcc                                                                                                              <  2:37485/149‑1 (MQ=255)
                                gcgcTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCt                                                                                                             <  2:2711403/149‑1 (MQ=255)
                                gcgcTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCt                                                                                                             >  1:2203276/1‑149 (MQ=255)
                                    tctTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAag                                                                                                         <  2:2568379/149‑1 (MQ=255)
                                     ctgtACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAaga                                                                                                        <  2:1362070/146‑1 (MQ=255)
                                       ttACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTCAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAgg                                                                                                      >  2:1180304/1‑149 (MQ=255)
                                       ttACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAgg                                                                                                      <  1:2843826/149‑1 (MQ=255)
                                        tACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGa                                                                                                     >  1:1350800/1‑149 (MQ=255)
                                         aCGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGaa                                                                                                    <  1:2706305/149‑1 (MQ=255)
                                           gCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATg                                                                                                  >  2:288981/1‑149 (MQ=255)
                                            cAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATgg                                                                                                 >  1:294791/1‑149 (MQ=255)
                                            cAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATgg                                                                                                 >  2:2408991/1‑149 (MQ=255)
                                            cAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATgg                                                                                                 >  2:869830/1‑149 (MQ=255)
                                              gATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGaa                                                                                               <  1:2139410/149‑1 (MQ=255)
                                               aTCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGaaa                                                                                              <  1:1756326/149‑1 (MQ=255)
                                                tCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAAc                                                                                             <  2:1001399/149‑1 (MQ=255)
                                                tCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAAc                                                                                             >  1:2634408/1‑149 (MQ=255)
                                                tCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAAc                                                                                             <  2:2627694/149‑1 (MQ=255)
                                                   gCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGt                                                                                          <  2:2984593/149‑1 (MQ=255)
                                                      tcatcCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAAt                                                                                       <  2:615137/149‑1 (MQ=255)
                                                         tcCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCtt                                                                                    <  1:3042142/149‑1 (MQ=255)
                                                          ccGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCttt                                                                                   <  1:1043661/149‑1 (MQ=255)
                                                              cAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGa                                                                               <  2:2011118/149‑1 (MQ=255)
                                                                   gTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGt                                                                          <  1:1298021/149‑1 (MQ=255)
                                                                   gTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGt                                                                          >  1:413447/1‑149 (MQ=255)
                                                                     tGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTc                                                                        <  2:3045143/149‑1 (MQ=255)
                                                                     tGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTc                                                                        <  2:1705284/149‑1 (MQ=255)
                                                                      gCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCa                                                                       <  1:1010813/149‑1 (MQ=255)
                                                                        ccGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCAct                                                                     <  1:953643/149‑1 (MQ=255)
                                                                        ccGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCAct                                                                     <  2:483391/149‑1 (MQ=255)
                                                                        ccGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCAGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCAct                                                                     <  1:2740401/149‑1 (MQ=255)
                                                                         cGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCActc                                                                    >  2:2667724/1‑149 (MQ=255)
                                                                         cGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCActc                                                                    >  2:2069337/1‑149 (MQ=255)
                                                                           gTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACtctt                                                                  >  1:2586370/1‑149 (MQ=255)
                                                                           gTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACtctt                                                                  <  1:2798172/149‑1 (MQ=255)
                                                                             aaaGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACtcttct                                                                <  2:956357/149‑1 (MQ=255)
                                                                                  cATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATTGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCg                                                           >  2:1823330/1‑149 (MQ=255)
                                                                                  cATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCg                                                           <  2:1122116/149‑1 (MQ=255)
                                                                                       aGAAACGCTTCGCACAGCACCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCAACGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTATTGCCGTGCTc                                                      >  2:1059042/1‑149 (MQ=255)
                                                                                         aaaCGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAg                                                    <  1:3127940/149‑1 (MQ=255)
                                                                                          aaCGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGc                                                   >  1:451256/1‑149 (MQ=255)
                                                                                             gCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTa                                                <  2:123644/149‑1 (MQ=255)
                                                                                              cTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTAt                                               <  2:2817833/149‑1 (MQ=255)
                                                                                              cTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTAt                                               <  1:771364/149‑1 (MQ=255)
                                                                                                 cGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATggg                                            <  2:1939359/149‑1 (MQ=255)
                                                                                                 cGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATggg                                            <  1:342878/149‑1 (MQ=255)
                                                                                                 cGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATggg                                            <  2:2440407/149‑1 (MQ=255)
                                                                                                   cacaGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCAACGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGAATGGGCt                                          >  1:3082783/1‑149 (MQ=255)
                                                                                                    acaGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTg                                         >  1:2286592/1‑149 (MQ=255)
                                                                                                    acaGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTg                                         <  2:1465307/149‑1 (MQ=255)
                                                                                                     cagcagCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGa                                        <  1:1397485/149‑1 (MQ=255)
                                                                                                     cagcagCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGa                                        <  2:2136403/149‑1 (MQ=255)
                                                                                                     cagcagCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGa                                        <  1:1636091/149‑1 (MQ=255)
                                                                                                         agcTCGCTGTTTGTAGCGTGACGCAGGTTGTTCGCCTGTCTTCGCGTTCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGatata                                    <  1:1716016/149‑1 (MQ=255)
                                                                                                           cTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATAtt                                  >  1:1919344/1‑149 (MQ=255)
                                                                                                           cTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATAtt                                  <  2:1581530/149‑1 (MQ=255)
                                                                                                           cTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATAtt                                  >  2:2194602/1‑149 (MQ=255)
                                                                                                            tCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTc                                 >  1:567260/1‑149 (MQ=255)
                                                                                                             cGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTcg                                >  2:325885/1‑149 (MQ=255)
                                                                                                               cTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTcgcc                              >  2:2243300/1‑149 (MQ=255)
                                                                                                                  ttGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCt                           <  2:2289590/149‑1 (MQ=255)
                                                                                                                   tGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTc                          >  1:2052757/1‑149 (MQ=255)
                                                                                                                   tGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTc                          >  1:496153/1‑149 (MQ=255)
                                                                                                                          gTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAg                   >  2:43637/1‑149 (MQ=255)
                                                                                                                               gCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATg              <  2:3098234/149‑1 (MQ=255)
                                                                                                                                   gttgttCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGAt          >  1:2527668/1‑149 (MQ=255)
                                                                                                                                    ttgttCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATa         >  1:3009343/1‑149 (MQ=255)
                                                                                                                                    ttgttCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATa         >  1:1250610/1‑149 (MQ=255)
                                                                                                                                     tgttCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATAc        <  1:2250891/149‑1 (MQ=255)
                                                                                                                                        tcgccTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTa     <  1:2440195/146‑1 (MQ=255)
                                                                                                                                          gccTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAAt   <  2:72508/148‑1 (MQ=255)
                                                                                                                                          gccTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAAt   <  1:1870408/148‑1 (MQ=255)
                                                                                                                                           ccTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTATCACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATc  >  2:661640/1‑149 (MQ=255)
                                                                                                                                           ccTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATc  >  1:1602288/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                
ATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATC  >  NC_000913/4128121‑4128408

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑