Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,708,597 T→C 5.8% intergenic (‑892/+19) dppA ← / ← proK dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor/tRNA‑Pro

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,708,5970TC 5.8% 401.7 / 15.5 179intergenic (‑892/+19)dppA/proKdipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor/tRNA‑Pro
Reads supporting (aligned to +/- strand):  ref base T (88/71);  new base C (5/5);  total (100/79)
Fisher's exact test for biased strand distribution p-value = 7.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-02

Only 100 of 182 total aligned reads displayed.

TCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGC  >  NC_000913/3708451‑3708745
                                                                                                                                                  |                                                                                                                                                    
      atatAGGATTACTCGCCAATTCGTGGCTCGACATGCGGGCCGTTGCTTGCGCAACGTTATCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCAACGAGCAGGc                                                                                                                                              <  2:16150/149‑1 (MQ=21)
      aTAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGACCCTCCCTGCCCCCGGCGGGGCGCATCCTCATGAGCTGCTGGTGCAAACAAGCCTACGGGTGGAGCAGGc                                                                                                                                              >  1:898747/1‑149 (MQ=255)
      aTAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCAGGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGc                                                                                                                                              >  1:975811/1‑149 (MQ=255)
        agagGATTCCTCGCCTCTGCGTGGCTCGCCCTGCGGGCCGTTGAGGGCGCAACGTTCTCTCGCTTCGCTCGCGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCATATGCAAAAAAGCCTGCTCGTTGAGCAGGCtt                                                                                                                                            <  2:373607/149‑1 (MQ=21)
          agGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGCGCCTCCCTGCCCCCGGCGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCGCGGTGAGCCGGCtttt                                                                                                                                          >  1:175807/1‑149 (MQ=25)
              ttACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGCGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGaa                                                                                                                                      >  2:559824/1‑149 (MQ=255)
                aCTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAAtt                                                                                                                                    <  1:986178/149‑1 (MQ=255)
                 cTCGCCACTTCGTGGCTCTCCCTGCGGTCCGGTGCTGGCGCAATGTTCTCTGGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAAttt                                                                                                                                   <  2:645396/149‑1 (MQ=18)
                 cTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCGGCCCCCGGAGGTTCGCAGCCTCATGAGCTGCAGCTGGAAAAAAGCCTGCGCGTTGGGCAGGGGGTTCGAAttt                                                                                                                                   >  1:618090/1‑149 (MQ=255)
                  tCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAGGCCTGCTCGGGGAGCAGGCGTTTCGAATTTg                                                                                                                                  >  1:509297/1‑149 (MQ=255)
                  tCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTg                                                                                                                                  >  2:71155/1‑149 (MQ=255)
                  tCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCAGCCTCATGAGCTGCCGCTGCAAAAAAGCCTGCTCGTGGCGCAGGCTTTTCGAATTTg                                                                                                                                  >  1:331066/1‑149 (MQ=25)
                    gCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGCGGTTCTCATCCTCATGAGCTGGAGATGCAAAAAAGCCTGCTCGTAGAGCAGGCTTTTCGAATTTGGt                                                                                                                                >  1:627273/1‑149 (MQ=255)
                    gCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGt                                                                                                                                >  2:597529/1‑149 (MQ=255)
                    gCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGt                                                                                                                                >  1:373641/1‑149 (MQ=255)
                    gCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGGTGAGCAGGCTTTTCGAATTTGGt                                                                                                                                >  2:417261/1‑149 (MQ=255)
                    gCCAATTCGGGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGt                                                                                                                                <  2:777970/149‑1 (MQ=255)
                     ccACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCCCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTc                                                                                                                               >  1:435666/1‑149 (MQ=255)
                        attCGTGGCTCGCCCGGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGt                                                                                                                            <  1:801998/148‑1 (MQ=255)
                           cGTGGCTCGCCCGGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGCGGTTCTCATCCTCATGAGCGGCAGCTGCAAAAAAGCCTGCTCGTTGAGCAGGCGTTTCGAATTTGGTCGGTGAt                                                                                                                         >  2:730052/1‑149 (MQ=21)
                             tGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGCGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATaa                                                                                                                       >  2:209916/1‑148 (MQ=255)
                                 tCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATGCGCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCCGGCTTGTCGAATTTGGTCGGTGATAGAGGa                                                                                                                   >  1:415858/1‑149 (MQ=255)
                                 tCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTAGAGCCGGCTTTTCGAATTTGGTCGGTGATAGAGGa                                                                                                                   >  2:1002494/1‑149 (MQ=255)
                                  cGCTCTTCGGGGCGTTGCTGTAGTATGGTTCTCTCGCTTCGGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGAt                                                                                                                  <  2:243551/149‑1 (MQ=14)
                                    cccTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGAtt                                                                                                                 >  1:726025/1‑148 (MQ=255)
                                     ccTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCGCCTCCTCAGGAGCGGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGGATTTGGTCGGTGATAGAGGATTCg                                                                                                               >  1:511999/1‑149 (MQ=18)
                                      cTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGCACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGGAGATGCAAAAAAGCCTGCTCGTTGAGCAGGGGTTTCGATTTTGGGGGGGGATTGAGGATTGGa                                                                                                              >  1:242205/1‑149 (MQ=18)
                                        gcgggcggTTCCTGGAGCCTCGTTATCTGGCATCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAAc                                                                                                            <  2:819236/149‑1 (MQ=14)
                                        gCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGATCGAGTCGAACCTCCCGGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAAc                                                                                                            <  1:504282/149‑1 (MQ=255)
                                          gggCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCt                                                                                                          <  1:877095/149‑1 (MQ=255)
                                           ggCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGCGTCGACCCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTGGAATTTGGTCGGTGGTAGAGGATTCGAGCCTc                                                                                                         >  1:687757/1‑149 (MQ=18)
                                             ccGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTGTCATCCGCATGAGCTGCCGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTAGGTGATAGAGGAGTCGAACCTCCg                                                                                                       >  1:664653/1‑149 (MQ=18)
                                              cGTTGCGGGCGCATTGTTATCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGATCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGa                                                                                                      <  2:110574/149‑1 (MQ=17)
                                                ttGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGAcc                                                                                                    <  1:597529/149‑1 (MQ=255)
                                                  gCTGGCGCAACGTTCTCTCGCTTCGGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGAcccc                                                                                                  <  1:422211/149‑1 (MQ=255)
                                                       cgcAACGTGCTCTCGCTTCGCTCGCGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGt                                                                                             <  2:212361/149‑1 (MQ=255)
                                                        gcAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTc                                                                                            >  2:416253/1‑149 (MQ=255)
                                                         gaaCGTTCTCTCGCTGCTCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTcc                                                                                           <  1:233813/148‑1 (MQ=255)
                                                          atcggtatctcGCGTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAACAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTccc                                                                                          <  2:573641/142‑1 (MQ=255)
                                                              ttCTCTCGCTTCGCTCGAGTCGAACCTCCGTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCcgaac                                                                                      <  2:339345/149‑1 (MQ=255)
                                                                 tctcGCTTAGGTCGAGGCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGGTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCcgaacgaa                                                                                   <  2:475326/149‑1 (MQ=255)
                                                                   tcgcttcgctCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGGAGATGCAAAAAAGCCTGCGCGATGAGCAGGCGTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTGCGTCCCGAACGAAGt                                                                                 >  1:776856/1‑149 (MQ=18)
                                                                       ttcgctCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGTTAGAGGATTCGAACCTCCGACCCCTTCGGCCCTAACGAAGTgcgg                                                                             >  2:688339/1‑148 (MQ=255)
                                                                       ttcgctCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTGGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTgcgc                                                                             >  1:128787/1‑149 (MQ=255)
                                                                         ggatctaGTCTAACCACCCTGCCCCAGTTGGATATCATCTTCATGATATGAAGATGGATAAAACGCTGATCGTTGATCAGGCTTTTCTAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTa                                                                           <  1:941650/143‑1 (MQ=16)
                                                                          gctCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTAc                                                                          <  2:820994/149‑1 (MQ=255)
                                                                            tCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTAccc                                                                        >  1:1009229/1‑148 (MQ=255)
                                                                             ggtgTCGAACATCCCTGCCCCCGGAGTTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCTTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACGAg                                                                       <  2:178289/146‑1 (MQ=18)
                                                                              gAGTCGAACGTCCCTTCCCCCGTAGGTTCTCATCCTAATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAgg                                                                      <  2:583086/149‑1 (MQ=255)
                                                                              gAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAgg                                                                      >  1:305899/1‑149 (MQ=255)
                                                                               aGTCGAACCTCCCTGCCCCCGGAGGTTCGCATCCTCATGAGCGGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGc                                                                     <  1:822763/149‑1 (MQ=255)
                                                                                gTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTAGTCCCGAACGAAGTGCGCTACCAGGCt                                                                    >  2:4269/1‑149 (MQ=255)
                                                                                 tCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAAAGAAGTGCGCTACCAGGCTg                                                                   >  2:329440/1‑149 (MQ=255)
                                                                                   gcaCCTCCCTGCCCCCGGAGGTTCTCATCAGCATGAGCTGCCGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGGGCGCTAACAGGCTgcg                                                                 >  1:961337/3‑149 (MQ=18)
                                                                                   gAACCTCCCTGCCCCCGGAGGTTCTCAGCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGCACGAAGTGCGCTACCAGGCTggg                                                                 >  2:328995/1‑147 (MQ=255)
                                                                                     aCCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGGTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGcc                                                               >  1:498754/1‑149 (MQ=255)
                                                                                       cTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCaa                                                             <  2:723260/149‑1 (MQ=255)
                                                                                         cccTGCCCCCGGCGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGCACCTCCGACCCCTTCGTCCCGCACGAAGTGCGCTACCAGGCGGCGCCACTc                                                           >  1:175009/1‑149 (MQ=18)
                                                                                         cccTGCCCCCGGCGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACATCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCCGGCTGCGCCAATc                                                           >  1:839121/1‑149 (MQ=255)
                                                                                         cccTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATc                                                           >  2:10494/1‑149 (MQ=255)
                                                                                           cttcccccGGATGTTCTCATTCTCATGATATGCAGATTCTAAATAGCCTGCTAGTTGAGTAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGTGTGCGCCAATCAc                                                         <  2:938512/146‑1 (MQ=21)
                                                                                           cTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCAc                                                         <  1:164365/149‑1 (MQ=255)
                                                                                            tGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCAcc                                                        <  2:247436/149‑1 (MQ=255)
                                                                                             gCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCg                                                       <  1:644305/149‑1 (MQ=255)
                                                                                               ccccGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCTGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGaa                                                     >  1:992551/1‑149 (MQ=255)
                                                                                               ccccGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGATCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCATCGTCACGAACCAAGTGCGCTACCACGCTGCGCCAATCACCGaa                                                     >  1:796807/1‑149 (MQ=21)
                                                                                                 ccGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTc                                                                                                                <  1:569352/88‑1 (MQ=31)
                                                                                                 ccGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATg                                                   >  1:1022410/1‑149 (MQ=255)
                                                                                                    gAGGTTCTCATCCTCATGAGCTGCAGATTCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCgg                                                <  1:985072/149‑1 (MQ=255)
                                                                                                          ctcATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCCCCGAATGCGGGGCGCa                                          >  2:765732/1‑149 (MQ=255)
                                                                                                          ctcATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCa                                          >  2:875956/1‑149 (MQ=255)
                                                                                                             aTCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGCCCCCTTCGTCCCCAACGAAGAGCGCTACCAGCCTGCGCCAATCACCGAATGCGGGGCGCAACt                                       >  1:1010942/1‑149 (MQ=25)
                                                                                                             aTCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCt                                       <  1:352462/149‑1 (MQ=255)
                                                                                                                 ccctGAGCTGCCGATGCAAAAAAGCCTGCTCGTAGAGCAGGCTTTTCGAATTTGGCCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACt                                   <  2:51801/146‑1 (MQ=255)
                                                                                                                  cATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTg                                  >  1:878339/1‑149 (MQ=255)
                                                                                                                   aTGAGCTGCAGATGCAAAAAATCCTGCTCGTTGAGCAGGCTTATCGAATTTTGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCACGAACGAAGTGCGCTACCAGGCGGCGCCAATCACCGAATGCGAGGCGCATCTTACTgc                                 >  2:190545/1‑149 (MQ=37)
                                                                                                                   aTGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTgc                                 >  2:380363/1‑149 (MQ=255)
                                                                                                                   aTGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTgc                                 <  1:206903/149‑1 (MQ=255)
                                                                                                                          gCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATa                          >  1:631236/1‑149 (MQ=255)
                                                                                                                          gCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCAACTTACTGCGCAGATa                          >  2:747826/1‑149 (MQ=255)
                                                                                                                           cAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATAc                         >  1:27463/1‑149 (MQ=255)
                                                                                                                                gCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGcccg                    >  1:270519/1‑148 (MQ=255)
                                                                                                                                 cAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTc                   <  2:435666/149‑1 (MQ=255)
                                                                                                                                  aaaaaaGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCg                  >  1:39626/1‑149 (MQ=255)
                                                                                                                                    aaacGCCTGCGCGTGGAGCAGGCTTTTCGAATATGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTc                <  2:945465/149‑1 (MQ=255)
                                                                                                                                    aaaaGCCTGCTCAACGAGCAGACTGTTCGAAGTTGGTCGGTGATAGAGGATTCGAGCCTCCGACCCCTTCGACCCGCACGCAGTGCGCTACCAGGCTGCGCCAATCACCGCATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCggc                >  2:215197/1‑147 (MQ=255)
                                                                                                                                       aGCCTGCTCAACGAGCAGGCTGTTCGAATTTGGTCGGTGATAGAGGATTCGCACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCTAATGCGGGGCGCAACTTACTGCGCAGATACGCCCTCGTCAAt             >  2:391129/1‑149 (MQ=38)
                                                                                                                                        gCCTGCTCAACGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCt                                                                                                 <  1:221262/64‑1 (MQ=18)
                                                                                                                                        gCCTGCTCAACGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCt                                                                                                 <  1:603939/64‑1 (MQ=18)
                                                                                                                                         ccTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATcc           >  1:923848/1‑149 (MQ=255)
                                                                                                                                         ccTGCTCAACGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATcc           >  1:684589/1‑149 (MQ=255)
                                                                                                                                          cTGCTCGTTGAGCAGGCTTTTCTAATTTGGTCGGTGATTGAGGATTCGAACCTCCGACCCCATCGTCCCGTACGAAGTGCGCTACCAGGCTGCGCCAAGCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATccc          <  2:966234/149‑1 (MQ=39)
                                                                                                                                           tGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCt         <  2:305899/149‑1 (MQ=255)
                                                                                                                                             cgcGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTa       <  2:600825/147‑1 (MQ=255)
                                                                                                                                              tCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTaa      <  1:71155/149‑1 (MQ=255)
                                                                                                                                              tCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTaa      <  1:74394/149‑1 (MQ=255)
                                                                                                                                               cGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAAt     <  1:417261/149‑1 (MQ=255)
                                                                                                                                                gTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATa    <  2:536236/149‑1 (MQ=255)
                                                                                                                                                 ttGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAg   <  2:618090/149‑1 (MQ=255)
                                                                                                                                                  tGAGCAGGCTTTTCGAATTTGGTCGGTGATCGAGGATTCGAACCTACGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGc  <  2:561455/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGAATGCGGGGCGCATCTTACTGCGCAGATACGCCCTCGTCAATCCCTTAATAGC  >  NC_000913/3708451‑3708745

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: