Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,788,719 T→A 5.0% H17L (CAC→CTC)  ynfL ← DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,788,7190TA 5.0% 248.4 / 2.1 79H17L (CAC→CTC) ynfLDNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (39/36);  new base A (2/2);  total (41/38)
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.36e-01

TCGAGGAACAATCGCCCGGCCTCGCTCAGTTCGACACGACGATTGGTACGGTCGAACAACCGGGCCCCCAACTCCTGCTCCAGAGCCTGGATCTGCTGACTCAAGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCAAGTCGTAAAACATCTCAAATAGGTCGAACAATATATTGGATCTAATCGTTAGCCAGCTATATGCTTTTTTCATTGCGCCAG  >  NC_002947/5788572‑5788864
                                                                                                                                                    |                                                                                                                                                 
tCGAGGAACAATCGCCCGGCCTCGCTCAGTTCGACACGACGATTGGTACGGTCGAACAACCGGGCCCCCAACTCCTGCTCCAGAGCCTGGATCTGCTGACTCAAGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTg                                                                                                                                                  >  7:311368/1‑149 (MQ=255)
     gAACAATCGCCCGGCCTCGCTCAGTTCGACACGACGATTGGTACGGTCGAACAACCGGGCCCCCAACTCCTGCTCCAGAGCCTGGATCTGCTGACTCAAGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGct                                                                                                                                             <  4:426828/149‑1 (MQ=255)
     gAACAATCGCCCGGCCTCGCTCAGTTCGACACGACGATTGGTACGGTCGAACAACCGGGCCCCCAACTCCTGCTCCAGAGCCTGGATCTGCTGACTCAAGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGct                                                                                                                                             <  3:13684/149‑1 (MQ=255)
     gAACAATCGCCAGGCCTAGCTCAGTTCGACACGACGATTGGTACGGTCGAACAACCGGGCCCCAAACTAATTCTCCAGAGACTGTATCTGCTGACTCAAGGGCGGTTGTTAAATACACAGCTGCTGGGCGGCGCGGCCAAAGagcagc                                                                                                                                              >  2:340639/1‑148 (MQ=255)
          aTCGCCCGGCCTCGCTCAGTTCGACACGACGATTGGTACGGTCGAACAACCGGGCCCCAAACTCCTGCTCCAGAGCGTGGATCTGCTGACTCAAGGCCGGTTGGGAAATACCCAGCTCCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCg                                                                                                                                        >  6:184193/1‑149 (MQ=255)
           tCGCCCGGCCTCGCTCAGTTCGACACGACGATTGGTACGGTCGAACAACCGGGCCCCCAACTCCTGCTCCAGAGCCTGGATCTGCTGACTCAAGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCgg                                                                                                                                       <  5:63492/149‑1 (MQ=255)
              cccGGCCTCGCTCAGTTCGACACGACGATTGGTACGGTCGAACAACCGGGCCCCCAACTCCTGCTCCAGAGCCTGGATCTGCTGACTCAAGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCca                                                                                                                                    <  5:48199/149‑1 (MQ=255)
              cccGGCCTCGCTCAGTTCGACACGACGATTGGTACGGTCGAACAACCGGGCCCCCAACTCCTGCTCCAGAGCCTGGATCTGCTGACTCAAGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAACTCCAGCTCTACGGCca                                                                                                                                    >  5:438122/1‑149 (MQ=255)
                   ccTCGCTCAGTTCGACACGACGATTGGTACGGTCGAACAACCGGGCCCCCAACTCCTGCTCCAGAGCCTGGATCTGCTGACTCAAGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGc                                                                                                                                >  6:181201/1‑148 (MQ=255)
                     tCGCTCAGTTCGACACGACGATTGGTACGGTCGAACAACCGGGCCCCCAACTCCTGCTCCAGAGCCTGGATCTGCTGACTCAAGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAAt                                                                                                                             >  3:410942/1‑149 (MQ=255)
                     tCGCTCAGTTCGACACGACGATTGGTACGGTCGAACAACCGGGCCCCCAACTCCTGCTCCAGAGCCTGGATCTGCTGACTCAAGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAAt                                                                                                                             >  3:219791/1‑149 (MQ=255)
                        cTCAGTTCGACACGACGATTGGTACGGTCGAACAACCGGGCCCCCAACTCCTGCTCCAGAGCCTGGATCTGCTGACTCAAGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAATGaa                                                                                                                          <  8:256547/149‑1 (MQ=255)
                        cTCAGTTCGACACGACGATTGGTACGGTCGAACAACCGGGCCCCCAACTCCTGCTCCAGAGCCTGGATCTGCTGACTCAAGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAATGaa                                                                                                                          <  3:159592/149‑1 (MQ=255)
                        cTCAGTTCGACACGACGATTGGTACGGTCGAACAACCGGGCCCCCAACTCCTGCTCCAGAGCCTGGATCTGCTGACTCAAGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAATGaa                                                                                                                          <  3:84439/149‑1 (MQ=255)
                           aGTTCGACACGACTCTTGCTACGGTCGAACAACCGGGCCCCCAACTCCTGCTCCAGAGCCTGGATCTGCTGACTCAAGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAATGAAg                                                                                                                         <  4:302406/147‑1 (MQ=255)
                              tCGACACGACGATTGGTACGGTCGAACAACCGGGCCCCCAACTCCTGCTCCAGAGCCTGGATCTGCTGACTCAAGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAATGAAGTAACg                                                                                                                    >  7:218717/1‑149 (MQ=255)
                                gACACGACGATTGGTACGGTCGAACAACCGGGCCCCCAACTCCTGCTCCAGAGCCTGGATCTGCTGACTCAAGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAAt                                                                                                                             <  6:341500/138‑1 (MQ=255)
                                gACACGACGATTGGTACGGTCGAACAACCGGGCCCCCAACTCCTGCTCCAGAGCCTGGATCTGCTGACTCAAGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAAt                                                                                                                             >  5:341500/1‑138 (MQ=255)
                                gACACGACGATTGGTACGGTCGAACAACCGGGCCCCCAACTCCTGCTCCAGAGCCTGGATCTGCTGACTCAAGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGc                                                                                                                   <  6:11536/148‑1 (MQ=255)
                                   acgacgATTGGTACGGTCGAACAACCGGGCCCCCAACTCCTGCTCCAGAGCCTGGATCTGCTGACTCAAGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGAt                                                                                                               <  8:218717/149‑1 (MQ=255)
                                     gacgaTTGGTACGGTCGAACAACCGGGCCCCCAACTCCTGCTCCAGAGCCTGGATCTGCTGACTCAAGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAATGaa                                                                                                                          <  7:398652/136‑1 (MQ=255)
                                     gacgaTTGGTACGGTCGAACAACCGGGCCCCCAACTCCTGCTCCAGAGCCTGGATCTGCTGACTCAAGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAATGaa                                                                                                                          >  8:398652/1‑136 (MQ=255)
                                           tGGTACGGTCGAACAACCGGGCCCCCAACTCCTGCTCCAGAGCCTGGATCTGCTGACTCAAGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGt                                                                                                       >  2:300059/1‑149 (MQ=255)
                                             gTACGGTCGAACAACCGGGCCCCCAACTCCTGCTCCAGAGCCTGGATCTGCTGACTCAAGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGtt                                                                                                      <  6:354647/148‑1 (MQ=255)
                                             gTACGGTCGAACAACCGGGCCCCCAACTCCTGCTCCAGAGCCTGGATCTGCTGACTCAAGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTc                                                                                                     >  2:240335/1‑149 (MQ=255)
                                                     gAACAACCGGGCCCCCAACTCCTGCTCCAGAGCCTGGATCTGCTGACTCAAGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGcc                                                                                             >  5:32492/1‑149 (MQ=255)
                                                           ccGGGCCCCCAACTCCTGCTCCAGAGCCTGGATCTGCTGACTCAAGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTcc                                                                                        >  8:276567/1‑148 (MQ=255)
                                                            cGGGCCCCCAACTCCTGCTCCAGAGCCTGGATCTGCTGACTCAAGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCaa                                                                                      <  1:240335/149‑1 (MQ=255)
                                                              gtcccccAACTCCTGATACAGAGCCTTTATCTGCTGACTCAATGGCGGTTGGGAAATACCCAGCTGCTTGGCGTCGCTTCCAAAGTGAATCTCTTCGGCCACAGCAATGAAGTAACGCCGATGACGCAGTTCCATGAGCCGCTCCAATc                                                                                    <  1:340639/147‑1 (MQ=255)
                                                                cccccAACTCCTGCTCCAGAGCCTGGATCTGCTGACTCAAGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGGTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCaa                                                                                  <  4:327679/149‑1 (MQ=255)
                                                                cccccAACTCCTGCTCCAGAGCCTGGATCTGCTGACTCAAGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCaa                                                                                  <  3:168890/149‑1 (MQ=255)
                                                                  cccAACTCCTGCTCCAGAGCCTGGATCTGCTGACTCAAGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCAAGt                                                                                >  4:218166/1‑149 (MQ=255)
                                                                  cccAACTCCTGCTCCAGAGACTGGATCTGCTGACTCAAGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAACTGCAGCTCTTCGGCCACAGCAAAGAAGTAACGCAGCTGACGCAGTTCCATGAGCCGCTCCAATCAAGt                                                                                >  5:325540/1‑149 (MQ=255)
                                                                   ccAACTCCTGCTCCAGAGCCTGGATCTGCTGACTCAAGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCAAGTc                                                                               <  7:276567/149‑1 (MQ=255)
                                                                             cTCCAGAGCCTGGATCTGCTGACTCAAGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCTGCCAAAGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCAAGTCGTAAAACa                                                                       <  6:322788/147‑1 (MQ=255)
                                                                                cAGAGCCTGGATCTGCTGACTCAAGGGCGGTTGGGAAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGATCTTCGCCAACAGAAATCAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCAAATCAATTCTTAAAACAtct                                                                    >  1:72584/1‑148 (MQ=255)
                                                                                cAGAGCCTGGATCTGCTGACTCAAGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTACATGAGCCGCTCCAATCAAGTCGTAAAACAtctc                                                                   >  5:364991/1‑148 (MQ=255)
                                                                                   agCCTGGATCTGCTGACTCAAGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCAAGTCGTAAAACATCTCAAAt                                                               <  5:181201/149‑1 (MQ=255)
                                                                                      cTGGATCTGCTGACTCAAGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCAAGTCGTAAAACATCTCAAATAgg                                                            >  4:37619/1‑149 (MQ=255)
                                                                                        ggATCTGCTGACTCAAGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCAAGTCGTAAAACATCTCAAATAGGTc                                                          >  1:16007/1‑149 (MQ=255)
                                                                                         gATCTGCTGACTCAAGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCAAGTCGTAAAACATCTCAAATAGGTc                                                          >  3:152057/1‑148 (MQ=255)
                                                                                              gctgACTCAAGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCAAGTCGTAAAACATCTCAAATAGGTCGAACaa                                                    >  2:119638/1‑149 (MQ=255)
                                                                                               ctgACTCAAGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCAAGTCGTAAAACATCTCAAATAGGTCGAACaa                                                    >  1:240032/1‑148 (MQ=255)
                                                                                                  aCTCAAGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCAAGTCGTAAAACATCTCAAATAGGTCGAACAatata                                                >  2:237399/1‑149 (MQ=255)
                                                                                                  aCTCAAGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGAAGTTCCATGAGCCGCGCCAATCAAATCGTAAAACATCTCAAAAACGTCGAACAatata                                                >  6:409850/1‑149 (MQ=255)
                                                                                                    gcacgggCGGTTCGGAAATACCCAGCTGCTGGGCGGCGCGGCCAACGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGCGCCGCTCCAATCAAGTCGTAAAACATCTCAAATAGGTCGAACAATATAtt                                              <  4:219791/145‑1 (MQ=255)
                                                                                                      aatggCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGTCCAAAGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCAAGTCGCAAAACATCTCAAATAGGTCGAACAATATATTgg                                            <  4:31938/146‑1 (MQ=255)
                                                                                                      aaGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCAAGTCGTAAAAc                                                                        >  4:53898/1‑121 (MQ=255)
                                                                                                      aaGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCAAGTCGTAAAAc                                                                        <  3:53898/121‑1 (MQ=255)
                                                                                                      aaGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCAAGTCGTAAAACATCTCAAATAGGTCGAACAATATATTg                                             >  4:379243/1‑148 (MQ=255)
                                                                                                        gggCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCAAGTCGTAAAACATCTCAAATAGGTCGAACAATATATTGGAt                                          <  3:218166/149‑1 (MQ=255)
                                                                                                                 ggAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCAAGTCGTAAAACATATCAAATATGTCGAACAATATATTGGATCTAATCGtt                                 >  4:230332/1‑149 (MQ=255)
                                                                                                                  gAAATACCCCTCTGCTGGGCGGCGCGTCCAAAGTGTAGCTCTTCGGCCACATCAATGAATTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATGAAGTCGTAAAACATCTCAAATAGGTCGAACAATATATTGTATCTAATCGTTa                                <  8:410790/149‑1 (MQ=255)
                                                                                                                    aaaTACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCATCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCAAGTCGTAAAACATCTCAAATAGGTCGAACAATATATTGGATCTAATCGTTAg                               <  3:320054/149‑1 (MQ=255)
                                                                                                                     aaTACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCAAGTCGTAAAACATCTCAAATAGGTCGAACAATATATTGGATCTAATCGTTAGc                              <  3:37619/149‑1 (MQ=255)
                                                                                                                          ccAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCAAGTCGTAAAACATCTCAAATAGGTCGAACaa                                                    >  8:96188/1‑122 (MQ=255)
                                                                                                                          ccAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCAAGTCGTAAAACATCTCAAATAGGTCGAACaa                                                    >  8:102084/1‑122 (MQ=255)
                                                                                                                          ccAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCAAGTCGTAAAACATCTCAAATAGGTCGAACaa                                                    <  7:96188/122‑1 (MQ=255)
                                                                                                                          ccAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCAAGTCGTAAAACATCTCAAATAGGTCGAACaa                                                    <  7:102084/122‑1 (MQ=255)
                                                                                                                          ccAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCAAGTCGTAAAACATCTCAAATAGGTCGAACAATATATTGGATCTAATCGTTAGCCAGc                          >  8:227206/1‑148 (MQ=255)
                                                                                                                          ccAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCAAGTCGTAAAACAACTAAAATAGGTCGAACAATATATTGGATCTAATCGTTAGCCAGc                          >  8:175682/1‑148 (MQ=255)
                                                                                                                                 ctgGGCGGCGCGGCCAAATTGCAGCTCTTCGGCCAAAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCAAGTCGTAAAACATCTCAAATAGGTCGAACAATATATTGGATCTAATCGTTAGCCAGCTatat                     <  5:207501/146‑1 (MQ=255)
                                                                                                                                    ggcggcGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCAAGTCGTAAAACATCTCAAATAGGTCGAACAATATATTGGATCTAATCGTTAGCCAGCTATATGCtt                 <  3:379243/147‑1 (MQ=255)
                                                                                                                                     gcggcgcggcCAAAGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCAAGTCGTAAAACATCTCAAATAGGTCGAACAATATATTGGATCTAATCGTTAGCCAGCTATATGCttttt              >  1:366233/1‑149 (MQ=255)
                                                                                                                                         tgtggcCGAAGAGCTGCACTTTGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCAAGTCGTAAAACATCTCAAATAGGTCGAACAATATATTGGATCTAATCGTTAGCCAGCTatat                     >  5:364215/4‑138 (MQ=255)
                                                                                                                                         tgtggcCGAAGAGCTGCACTTTGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCAAGTCGTAAAACATCTCAAATAGGTCGAACAATATATTGGATCTAATCGTTAGCCAGCTatat                     <  6:364215/135‑1 (MQ=255)
                                                                                                                                         tgtggcCGAAGAGCTGCACTTTGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCAAGTCGTAAAACATCTCAAATAGGTCGAACAATATATTGGATCTAATCGTTAGCCAGCTatat                     <  1:419246/135‑1 (MQ=255)
                                                                                                                                         tgtggcCGAAGAGCTGCACTTTGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCAAGTCGTAAAACATCTCAAATAGGTCGAACAATATATTGGATCTAATCGTTAGCCAGCTatat                     >  2:419246/4‑138 (MQ=255)
                                                                                                                                          gcgtcCAATGTGCAGCTCTTCGGCCCCAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTGCAATCAAGTCGTAAAACATCTCATATAGGTCGAACAATATATTGGATCTAATCGTTAGCCAGCTATATGCTTTTTTCAtt         <  6:325540/149‑1 (MQ=255)
                                                                                                                                          gcggcCAAAGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCAAGTCGTAAAACATCTCAAATAGGTCGAACAATATATTGGATCTAATCGTTAGCCAGCTATATGCTTTTTTCAtt         >  8:66423/1‑149 (MQ=255)
                                                                                                                                             gcCAAAGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCAAGTCGTAAAACATCTCAAATAgg                                                            <  3:84994/95‑1 (MQ=255)
                                                                                                                                             gcCAAAGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCAAGTCGTAAAACATCTCAAATAgg                                                            >  4:84994/1‑95 (MQ=255)
                                                                                                                                             gcCAAAGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCAAGTCGTAAAACATCTCAAATAGGTCGAACAATATATTGGATCTAATCGTTAGCCAGCTATATGCTTTTTTCATTgcg      >  3:154029/1‑149 (MQ=255)
                                                                                                                                               cAAAGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCAAGTCGTAAAACATCTCAAATAGGTCGAACAATATATTGGATCTAATCGTTAGCCAGCTATATGCTTTTTTCATTgcgc     >  7:204741/1‑148 (MQ=255)
                                                                                                                                                aaaGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCAAGTCGTAAAAGATCTCAAATAGGTCGAACAATATATTGGATCTAATCGTTAGCCAGCTATATCCttttt              <  1:75200/138‑1 (MQ=255)
                                                                                                                                                aaaGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCAAGTCGTAAAACATCTCAAATAGGTCGAACAAAATAATGGATCTAAACGTTAGCCAGCTAAATACttatt              >  2:75200/1‑138 (MQ=255)
                                                                                                                                                 aaGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCAAGTCGTAAAACATCTCAAATAGGTCGAACAATATATTGGATCTAATCGTTAGCCAGCTATATGCTTTTTTCAt          >  3:65290/1‑141 (MQ=255)
                                                                                                                                                 aaGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCAAGTCGTAAAACATCTCAAATAGGTCGAACAATATATTGGATCTAATCGTTAGCCAGCTATATGCTTTTTTCAt          <  4:65290/141‑1 (MQ=255)
                                                                                                                                                 aaGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCAAGTCGTAAAACATCTCAAATAGGTCGAACAATATATTGGATCTAATCGTTAGCCAGCTATATGCTTTTTTCATTGCGCCag  <  6:32492/149‑1 (MQ=255)
                                                                                                                                                  aGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCAAGTCGTAAAACATCTCAAATAGGTCGAACAATATATTGGATCTAATCGTTAGCCAGCTATATGCTTTTTTCATTgcgc     >  4:292245/1‑145 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TCGAGGAACAATCGCCCGGCCTCGCTCAGTTCGACACGACGATTGGTACGGTCGAACAACCGGGCCCCCAACTCCTGCTCCAGAGCCTGGATCTGCTGACTCAAGGGCGGTTGGGAAATACCCAGCTGCTGGGCGGCGCGGCCAAAGTGCAGCTCTTCGGCCACAGCAATGAAGTAACGCAGATGACGCAGTTCCATGAGCCGCTCCAATCAAGTCGTAAAACATCTCAAATAGGTCGAACAATATATTGGATCTAATCGTTAGCCAGCTATATGCTTTTTTCATTGCGCCAG  >  NC_002947/5788572‑5788864

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