Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,128,259 A→G V46A (GTG→GCG)  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,2590AG99.8% 6091.7 / ‑4.7 1619V46A (GTG→GCG) metJDNA‑binding transcriptional repressor MetJ
Reads supporting (aligned to +/- strand):  ref base A (0/1);  major base G (746/869);  minor base T (3/0);  total (749/870)
Fisher's exact test for biased strand distribution p-value = 9.90e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 1619 total aligned reads displayed.

CACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTT  >  NC_000913/4128113‑4128404
                                                                                                                                                  |                                                                                                                                                 
  cGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTg                                                                                                                                               >  2:2934428/1‑149 (MQ=255)
   gCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAATGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTgt                                                                                                                                              >  2:1859805/1‑148 (MQ=255)
   gCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTCCGCCTgt                                                                                                                                              <  2:1643954/149‑2 (MQ=255)
    cATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTgtc                                                                                                                                             >  1:81581/1‑147 (MQ=255)
     aTGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTgtct                                                                                                                                            >  1:2127477/1‑146 (MQ=255)
        aTCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGtctacg                                                                                                                                         <  2:1209198/149‑1 (MQ=255)
        aTCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGtctacg                                                                                                                                         <  1:1187183/149‑1 (MQ=255)
        aTCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGtctacg                                                                                                                                         >  1:1803882/1‑149 (MQ=255)
        aTCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGtctacg                                                                                                                                         <  2:2391599/149‑1 (MQ=255)
         tctcTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTAcgc                                                                                                                                        >  2:604310/1‑149 (MQ=255)
         tctcTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTAcgc                                                                                                                                        >  2:1578674/1‑149 (MQ=255)
             ttttGCCGCGTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACgcgtgc                                                                                                                                    <  1:330586/149‑1 (MQ=255)
              tttGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACgcgtgcg                                                                                                                                   >  1:1478519/1‑149 (MQ=255)
              tttGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACgcgtgcg                                                                                                                                   <  2:2437578/149‑1 (MQ=255)
                  ccGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCa                                                                                                                               >  1:376364/1‑149 (MQ=255)
                  ccGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCa                                                                                                                               >  1:1864012/1‑149 (MQ=255)
                    gCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATc                                                                                                                             >  2:2755865/1‑149 (MQ=255)
                       tCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGt                                                                                                                          <  1:2711170/149‑1 (MQ=255)
                           ggATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAgg                                                                                                                      >  2:477932/1‑149 (MQ=255)
                            gATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGa                                                                                                                     >  2:3094860/1‑149 (MQ=255)
                            gATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGa                                                                                                                     <  2:2809147/149‑1 (MQ=255)
                             aTTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGAt                                                                                                                    <  1:694403/149‑1 (MQ=255)
                              tttCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGAtt                                                                                                                   >  2:1769478/1‑149 (MQ=255)
                              tttCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGAtt                                                                                                                   <  1:1989473/149‑1 (MQ=255)
                               ttCGTCGCTGCGCTCTTTCCGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGAttt                                                                                                                  <  2:2060649/149‑1 (MQ=255)
                               ttCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGAttt                                                                                                                  <  2:1150401/149‑1 (MQ=255)
                                tcgtcgCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGAtttt                                                                                                                 <  2:1911177/149‑1 (MQ=255)
                                 cgtcgCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGAttttt                                                                                                                <  1:1528085/149‑1 (MQ=255)
                                  gtcgCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTa                                                                                                               <  2:1497880/149‑1 (MQ=255)
                                    cgCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAac                                                                                                             <  2:1292091/149‑1 (MQ=255)
                                                 aCGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGCTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGaa                                                                                                <  1:3023804/149‑1 (MQ=255)
                                                  cGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAAt                                                                                               >  2:769966/1‑149 (MQ=255)
                                                  cGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAAt                                                                                               <  2:2279605/149‑1 (MQ=255)
                                                   gCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATg                                                                                              <  1:2924122/149‑1 (MQ=255)
                                                    cAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATgg                                                                                             >  1:1394619/1‑149 (MQ=255)
                                                      gATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGaa                                                                                           <  2:233780/149‑1 (MQ=255)
                                                           gCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGt                                                                                      <  2:2235337/149‑1 (MQ=255)
                                                            catcatCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTa                                                                                     <  1:2569641/149‑1 (MQ=255)
                                                             atcatcCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTaa                                                                                    <  1:2404958/149‑1 (MQ=255)
                                                              tcatcCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAAt                                                                                   >  1:486573/1‑149 (MQ=255)
                                                                 tcCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCtt                                                                                >  1:971600/1‑149 (MQ=255)
                                                                 tcCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCtt                                                                                >  1:2195484/1‑149 (MQ=255)
                                                                  ccGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCttt                                                                               <  1:1043677/149‑1 (MQ=255)
                                                                   cGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCtttt                                                                              <  2:1984803/149‑1 (MQ=255)
                                                                   cGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCtttt                                                                              <  1:2763258/149‑1 (MQ=255)
                                                                   cGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCtttt                                                                              <  1:306302/149‑1 (MQ=255)
                                                                   cGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCtttt                                                                              <  2:43498/149‑1 (MQ=255)
                                                                   cGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCtttt                                                                              <  2:2977861/149‑1 (MQ=255)
                                                                   cGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCtttt                                                                              <  2:2666094/149‑1 (MQ=255)
                                                                     gCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTg                                                                            >  1:437476/1‑149 (MQ=255)
                                                                     gCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTg                                                                            <  2:1393014/149‑1 (MQ=255)
                                                                         aGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACtt                                                                        <  2:1462133/149‑1 (MQ=255)
                                                                           gTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGt                                                                      <  1:246435/149‑1 (MQ=255)
                                                                           gTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGt                                                                      <  1:2345194/149‑1 (MQ=255)
                                                                           gTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGt                                                                      <  2:2188136/149‑1 (MQ=255)
                                                                           gTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGt                                                                      <  1:1298021/149‑1 (MQ=255)
                                                                           gTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGt                                                                      <  2:826983/149‑1 (MQ=255)
                                                                             tGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTc                                                                    <  2:3045143/149‑1 (MQ=255)
                                                                             tGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTc                                                                    <  2:660300/149‑1 (MQ=255)
                                                                              gCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCa                                                                   <  1:1010813/149‑1 (MQ=255)
                                                                                ccGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCAct                                                                 <  1:747/149‑1 (MQ=255)
                                                                                ccGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCAct                                                                 >  1:625259/1‑149 (MQ=255)
                                                                                 cGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCActc                                                                >  1:2059195/1‑149 (MQ=255)
                                                                                 cGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCCTTTTGACTTGTTCActc                                                                >  2:1042878/1‑149 (MQ=255)
                                                                                  ggTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCActct                                                               >  2:2493882/1‑149 (MQ=255)
                                                                                     aaaGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACtcttct                                                            <  1:1605509/149‑1 (MQ=255)
                                                                                     aaaGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACtcttct                                                            <  1:1639848/149‑1 (MQ=255)
                                                                                        ggCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGc                                                         <  2:1090942/149‑1 (MQ=255)
                                                                                        ggCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGc                                                         <  1:2849368/149‑1 (MQ=255)
                                                                                          cATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCg                                                       >  1:933642/1‑149 (MQ=255)
                                                                                          cATGCAGAAACGCTTCGCACAGCAGCTCGCGGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCGTTTTGACTTGTTCACTCTTCTTGCCg                                                       >  1:1091301/1‑149 (MQ=255)
                                                                                               aGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTc                                                  <  2:952843/149‑1 (MQ=255)
                                                                                               aGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTc                                                  >  2:1334708/1‑149 (MQ=255)
                                                                                                gAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCa                                                 >  1:1348201/1‑149 (MQ=255)
                                                                                                 aaaCGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAg                                                <  1:2269505/149‑1 (MQ=255)
                                                                                                   aCGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCg                                              <  2:2876278/149‑1 (MQ=255)
                                                                                                     gCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGATGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTa                                            <  1:282753/149‑1 (MQ=255)
                                                                                                       ttCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATg                                          <  2:713393/149‑1 (MQ=255)
                                                                                                          gCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGc                                       <  1:881883/149‑1 (MQ=255)
                                                                                                           cacaGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCt                                      >  2:1686418/1‑149 (MQ=255)
                                                                                                            acaGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTg                                     <  1:311398/149‑1 (MQ=255)
                                                                                                            acaGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTg                                     <  2:2704974/149‑1 (MQ=255)
                                                                                                             cagcagCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGa                                    <  2:2136403/149‑1 (MQ=255)
                                                                                                               gcagcTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGata                                  >  1:1712259/1‑149 (MQ=255)
                                                                                                                cagcTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGatat                                 >  2:869248/1‑149 (MQ=255)
                                                                                                                   cTCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATAtt                              >  1:1919344/1‑149 (MQ=255)
                                                                                                                    tCGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTc                             >  1:1752409/1‑149 (MQ=255)
                                                                                                                     cGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTcg                            >  2:325885/1‑149 (MQ=255)
                                                                                                                     cGCTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTcg                            >  2:262734/1‑149 (MQ=255)
                                                                                                                       cTGTTGGTAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATCcgcc                          <  2:537234/149‑1 (MQ=255)
                                                                                                                              tAGCGTGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCAt                   <  1:540670/149‑1 (MQ=255)
                                                                                                                                   tGACGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGc              >  2:2563818/1‑149 (MQ=255)
                                                                                                                                     aCGCAGGTTGTTCGCCTGTCTACGTGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCa            <  1:2642110/149‑1 (MQ=255)
                                                                                                                                     aCGCAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCa            <  2:1814499/149‑1 (MQ=255)
                                                                                                                                      cGCAGGTTGTTCGCCTGTCTGCGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCAt           >  1:1139556/1‑149 (MQ=255)
                                                                                                                                        cAGGTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATga         >  1:2603827/1‑149 (MQ=255)
                                                                                                                                          ggTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATgaga       >  2:374543/1‑149 (MQ=255)
                                                                                                                                          ggTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATgaga       >  1:2961646/1‑149 (MQ=255)
                                                                                                                                          ggTTGTTCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATgaga       >  1:223597/1‑149 (MQ=255)
                                                                                                                                            ttgttCGCCTGTCTACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATa     >  1:1250679/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTT  >  NC_000913/4128113‑4128404

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: