Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,128,316 A→G V27A (GTT→GCT)  metJ ← transcriptional repressor, S‑adenosylmethionine‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,128,3160AG99.3% 528.6 / NA 140V27A (GTT→GCT) metJtranscriptional repressor, S‑adenosylmethionine‑binding
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (74/65);  minor base C (1/0);  total (75/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 140 total aligned reads displayed.

GGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTA  >  NC_000913/4128171‑4128463
                                                                                                                                                 |                                                                                                                                                   
ggCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACTCCTTAAGAGGCATGGAAGCCGTa                                                                                                                                                <  4:188465/151‑1 (MQ=255)
ggCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTa                                                                                                                                                <  2:861507/151‑1 (MQ=255)
   atcatcCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATc                                                                                                                                             >  1:1007709/1‑151 (MQ=255)
   atcatcCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAGCCGTAATc                                                                                                                                             <  1:53312/151‑1 (MQ=255)
     catcCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCtt                                                                                                                                           >  1:705164/1‑151 (MQ=255)
     catcCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCtt                                                                                                                                           <  1:1233725/151‑1 (MQ=255)
      atcCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCttt                                                                                                                                          >  1:687815/1‑151 (MQ=255)
       tcCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCtttt                                                                                                                                         >  2:333928/1‑151 (MQ=255)
       tcCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGAGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCtttt                                                                                                                                         >  1:537631/1‑151 (MQ=255)
        ccGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCttttt                                                                                                                                        >  1:605903/1‑151 (MQ=255)
        ccGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCttttt                                                                                                                                        <  1:1699030/151‑1 (MQ=255)
         cGGCAAAGGTTGCCCGGTAAAGGCATGCGGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTg                                                                                                                                       >  2:1007830/1‑151 (MQ=255)
         cGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAATGCTTCGCACAGCAGCTCGCTGTTGGCAGCGAGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCtttt                                                                                                                                         >  1:73762/1‑149 (MQ=255)
         cGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTg                                                                                                                                       <  2:1393141/151‑1 (MQ=255)
            cAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACt                                                                                                                                    >  2:1525866/1‑151 (MQ=255)
            cAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACt                                                                                                                                    <  4:116177/151‑1 (MQ=255)
                ggTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGtt                                                                                                                                <  2:299309/151‑1 (MQ=255)
                   tGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCa                                                                                                                              <  1:1147979/150‑1 (MQ=255)
                        ggTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACtctt                                                                                                                         <  2:304741/150‑1 (MQ=255)
                           aaaGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTcttctt                                                                                                                      <  1:1240762/150‑1 (MQ=255)
                           aaaGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTcttctt                                                                                                                      <  2:1415802/150‑1 (MQ=255)
                            aaGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGc                                                                                                                    >  4:60810/1‑151 (MQ=255)
                             aGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGcc                                                                                                                   >  3:70292/1‑151 (MQ=255)
                                cATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTg                                                                                                                >  2:1622673/1‑151 (MQ=255)
                                 aTGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACAACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGc                                                                                                               >  1:1356511/1‑151 (MQ=255)
                                   gCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTc                                                                                                             <  2:1523801/151‑1 (MQ=255)
                                          cGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTACGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTAt                                                                                                      >  1:1695625/1‑151 (MQ=255)
                                             ttCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATggg                                                                                                   >  1:1357259/1‑151 (MQ=255)
                                               cGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCt                                                                                                 <  1:926434/151‑1 (MQ=255)
                                                  acaGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGat                                                                                              <  2:531070/151‑1 (MQ=255)
                                                    agcagcTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGatat                                                                                            <  1:561892/151‑1 (MQ=255)
                                                    agcagcTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGatat                                                                                            >  1:1276415/1‑151 (MQ=255)
                                                      cagcTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGatatat                                                                                          >  2:430212/1‑151 (MQ=255)
                                                        gcTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTc                                                                                        <  1:468982/151‑1 (MQ=255)
                                                        gcTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTc                                                                                        <  1:1520479/151‑1 (MQ=255)
                                                            gCTGTTGGAAGCGTGACGCAGGATGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTcgccg                                                                                    >  1:433522/1‑151 (MQ=255)
                                                             cTGTTGGTAGCGTGACGCGGGTTGTTCACCTGAAGACGCGTGCGTTCATCGGGGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTcgccg                                                                                    >  1:900983/1‑150 (MQ=255)
                                                             cTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTcgccgc                                                                                   >  2:964646/1‑151 (MQ=255)
                                                              tGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCt                                                                                  >  2:1347120/1‑151 (MQ=255)
                                                              tGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCt                                                                                  >  1:341279/1‑151 (MQ=255)
                                                                  ggTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCAt                                                                              <  4:459950/151‑1 (MQ=255)
                                                                      gCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAg                                                                          <  1:1057387/151‑1 (MQ=255)
                                                                          gACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCa                                                                       <  1:1185969/150‑1 (MQ=255)
                                                                           aCGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATg                                                                     <  2:851837/151‑1 (MQ=255)
                                                                            cGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATga                                                                    <  2:399236/151‑1 (MQ=255)
                                                                             gCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATgag                                                                   <  1:5698/151‑1 (MQ=255)
                                                                                ggTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATa                                                                <  2:1634001/151‑1 (MQ=255)
                                                                                ggTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATa                                                                <  1:1544668/151‑1 (MQ=255)
                                                                                  ttgttCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACt                                                              >  3:243725/1‑151 (MQ=255)
                                                                                   tgttCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACtt                                                             >  3:544308/1‑151 (MQ=255)
                                                                                   tgttCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACtt                                                             >  4:162523/1‑151 (MQ=255)
                                                                                   tgttCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACt                                                              <  2:643023/150‑1 (MQ=255)
                                                                                    gttCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACtt                                                             >  2:1052524/1‑150 (MQ=255)
                                                                                       cACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAAt                                                          <  2:822476/150‑1 (MQ=255)
                                                                                        aCCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATcc                                                        <  1:1639289/151‑1 (MQ=255)
                                                                                         ccTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCct                                                       >  1:771807/1‑151 (MQ=255)
                                                                                             acgacgCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTc                                                   <  1:639244/151‑1 (MQ=255)
                                                                                             acgacgCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTc                                                   >  1:438214/1‑151 (MQ=255)
                                                                                                acgCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTc                                                >  2:199938/1‑151 (MQ=255)
                                                                                                acgCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTc                                                <  1:1201660/151‑1 (MQ=255)
                                                                                                acgCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTc                                                <  3:401941/151‑1 (MQ=255)
                                                                                                  gcgtgcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCa                                               <  1:1426000/150‑1 (MQ=255)
                                                                                                    gtgcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAAt                                             >  1:1418215/1‑150 (MQ=255)
                                                                                                    gtgcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATa                                            <  2:100750/151‑1 (MQ=255)
                                                                                                      gcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATaa                                           <  2:532864/150‑1 (MQ=255)
                                                                                                      gcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTGGACGTGTGCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATaa                                           >  2:393260/1‑150 (MQ=255)
                                                                                                       cgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAAt                                         <  3:15302/151‑1 (MQ=255)
                                                                                                          tCATCGGCGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAAt                                         >  4:80010/1‑148 (MQ=255)
                                                                                                            aTCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTg                                       >  2:1356358/1‑148 (MQ=255)
                                                                                                             tCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagac                                   >  1:1130160/1‑151 (MQ=255)
                                                                                                             tCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagac                                   >  2:880209/1‑151 (MQ=255)
                                                                                                                gTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccag                                <  1:1654750/151‑1 (MQ=255)
                                                                                                                gTGAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagacc                                  >  2:484035/1‑149 (MQ=255)
                                                                                                                  gAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccagac                              >  3:178883/1‑151 (MQ=255)
                                                                                                                  gAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccagac                              <  4:424318/151‑1 (MQ=255)
                                                                                                                  gAGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccagac                              >  2:1256222/1‑151 (MQ=255)
                                                                                                                   aGGATTTTTAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccagacc                             >  3:135869/1‑151 (MQ=255)
                                                                                                                       tttttAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACcacac                         >  3:379908/1‑151 (MQ=255)
                                                                                                                        ttttAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACcacac                         >  1:1193909/1‑150 (MQ=255)
                                                                                                                         tttAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACcacac                         >  1:942583/1‑149 (MQ=255)
                                                                                                                          ttAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAg                       <  1:651596/150‑1 (MQ=255)
                                                                                                                          ttAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGt                      >  2:1263168/1‑151 (MQ=255)
                                                                                                                           tAACACCTTCAGAGGAATGGAAGCCGTAATCTTTTTCACTTGTTCACTCTTCTTGCCTTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCCGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGTCCACACAGt                      <  1:1578108/150‑1 (MQ=255)
                                                                                                                           tAACACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGtt                     >  2:469857/1‑151 (MQ=255)
                                                                                                                            aaCACCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGtt                     <  2:1207137/150‑1 (MQ=255)
                                                                                                                               acCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGAtg                 <  2:664581/151‑1 (MQ=255)
                                                                                                                               acCTTAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGAtg                 >  2:815511/1‑151 (MQ=255)
                                                                                                                                   tAAGAGGACTGGAACCCGCAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTggg             >  4:472645/1‑151 (MQ=255)
                                                                                                                                   tAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGAtgtg               <  1:80537/149‑1 (MQ=255)
                                                                                                                                   tAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTggg             >  3:162562/1‑151 (MQ=255)
                                                                                                                                   tAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTggg             >  1:1009356/1‑151 (MQ=255)
                                                                                                                                   tAAGAGGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTgg              <  2:705267/150‑1 (MQ=255)
                                                                                                                                       agGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACt         >  4:320235/1‑151 (MQ=255)
                                                                                                                                       agGAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACt         <  2:149608/151‑1 (MQ=255)
                                                                                                                                        ggAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCGTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACt         <  1:284553/150‑1 (MQ=255)
                                                                                                                                         gAATGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTg        >  2:329186/1‑150 (MQ=255)
                                                                                                                                          aaTGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGAc      >  1:1694515/1‑151 (MQ=255)
                                                                                                                                           aTGGAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGAcc     <  2:115968/151‑1 (MQ=255)
                                                                                                                                             ggAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGt   >  4:374370/1‑151 (MQ=255)
                                                                                                                                              gAAGCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTa  <  2:1356659/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTA  >  NC_000913/4128171‑4128463

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