Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,988,365 G→C 7.0% intergenic (+35/+59) ycgR → / ← PP_4398 flagellar brake protein YcgR/MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,988,3650GC 7.0% 228.0 / 6.5 71intergenic (+35/+59)ycgR/PP_4398flagellar brake protein YcgR/MFS transporter
Reads supporting (aligned to +/- strand):  ref base G (32/34);  new base C (3/2);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.80e-01

TTTGCCGGCGTGCGTTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAGCCGCCACCGGCGCCTGC  >  NC_002947/4988220‑4988511
                                                                                                                                                 |                                                                                                                                                   
tttGCCGGCGTGCGTTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAgg                                                                                                                                                  >  1:281276/1‑149 (MQ=255)
       gcgtgcgtTTCCACAACCTGAGCGGGCAGGCTCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTCAGCAGGTGGCGtt                                                                                                                                           <  2:279684/149‑1 (MQ=255)
       gcgtgcgtTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGtt                                                                                                                                           >  3:323280/1‑149 (MQ=255)
          tgcgtTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTc                                                                                                                                       <  3:161051/149‑1 (MQ=255)
           gcgtTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTc                                                                                                                                       <  4:178739/148‑1 (MQ=255)
             gtTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCAtt                                                                                                                                    >  8:297727/1‑149 (MQ=255)
             gtTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTCAGCAGGTGGCGTTTTctgt                                                                                                                                    <  7:138576/149‑4 (MQ=255)
             ggttCCACAACCTGAGCGGGCAGGCGCAGCGCAAGAGCGAGAGGTTCGTGTACCAGCTGCAGCGTGTGGCGCGGCGGTTCGACAAGGACGAGTATTGATCGGCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCAtt                                                                                                                                    <  8:267351/147‑1 (MQ=255)
              tttCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCg                                                                                                                                             <  4:143857/140‑1 (MQ=255)
              tttCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCg                                                                                                                                             >  3:143857/1‑140 (MQ=255)
              tttCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCAttt                                                                                                                                   <  4:96584/149‑1 (MQ=255)
               ttCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTtg                                                                                                                                  >  6:131099/1‑149 (MQ=255)
                  cacaACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTgtgt                                                                                                                               >  3:198295/1‑149 (MQ=255)
                   acaACCTGAGCGGGCTGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCAttt                                                                                                                                   <  7:77364/144‑1 (MQ=255)
                   acaACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCAttt                                                                                                                                   >  8:77364/1‑144 (MQ=255)
                   acaACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGtt                                                                                                                              >  7:43694/1‑149 (MQ=255)
                    caACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCAttt                                                                                                                                   >  2:18744/1‑143 (MQ=255)
                    caACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCAttt                                                                                                                                   <  1:18744/143‑1 (MQ=255)
                       ccTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTctg                                                                                                                          >  8:137657/1‑149 (MQ=255)
                       ccTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTctg                                                                                                                          <  4:198295/149‑1 (MQ=255)
                       ccTGAGCGGGCAGGCGCAGCGCAACATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTctg                                                                                                                          >  7:143745/1‑149 (MQ=255)
                                  aGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGtttt                                                                                                                                         <  6:11176/124‑1 (MQ=255)
                                  aGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGtttt                                                                                                                                         >  5:11176/1‑124 (MQ=255)
                                  aGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGcc                                                                                                               >  5:217320/1‑149 (MQ=255)
                                  aGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGcc                                                                                                               >  4:161922/1‑149 (MQ=255)
                                    gcgcagcgcaAGATCGAGAGTTTCGTGTACCAGCTGCAGAGTGAGGCGCGGCGGTTCGACAAGGAAGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGc                                                                                                             >  4:95267/1‑149 (MQ=255)
                                        agcgcaAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTcg                                                                                                         >  5:65233/1‑149 (MQ=255)
                                           gcaAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCgg                                                                                                      <  1:152595/149‑1 (MQ=255)
                                           gcaAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCgg                                                                                                      <  3:161922/149‑1 (MQ=255)
                                           ctaaGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCgg                                                                                                      <  5:103385/147‑1 (MQ=255)
                                              aGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCgg                                                                                                      >  6:103385/1‑146 (MQ=255)
                                                 tCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAAc                                                                                                <  6:65233/149‑1 (MQ=255)
                                                       gTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTc                                                                                          <  7:137657/149‑1 (MQ=255)
                                                               tACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACAcc                                                                                   >  1:261876/1‑148 (MQ=255)
                                                               tACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACAcc                                                                                   >  5:71466/1‑148 (MQ=255)
                                                                 ccAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGa                                                                                 <  7:297727/148‑1 (MQ=255)
                                                                 ccAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGAc                                                                                <  7:263318/149‑1 (MQ=255)
                                                                   aGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCATCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACtt                                                                              <  3:113344/149‑1 (MQ=255)
                                                                       gCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTc                                                                                          >  5:13999/1‑133 (MQ=255)
                                                                       gCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTc                                                                                          <  6:13999/133‑1 (MQ=255)
                                                                        cAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCt                                                                           <  2:281276/147‑1 (MQ=255)
                                                                          gCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGAc                                                                        >  6:256782/1‑148 (MQ=255)
                                                                                ggCGCGGCGGTTCGACAAGGACGATTTTTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGTGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCa                                                                 <  8:280888/149‑1 (MQ=255)
                                                                                   gcggcgGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGgc                                                              >  4:19686/1‑149 (MQ=255)
                                                                                   gcggcgGTTCGACAAGGAAGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGgc                                                              <  3:95267/149‑1 (MQ=255)
                                                                                    cggcggTTCGAAAAGGAAGATTATAGATAGTAGGTCGGTAAAAAAGAAAACGCCACCTGCTGAGGTGTCGTTTTCATTTGTGTTTCTGCTGGGCTCGCCGCTTCGCGGGTAAACCCGCTACTACACATACTTGCTGACCTCTTCAGccg                                                             >  6:324717/1‑149 (MQ=255)
                                                                                            cGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTtg                                                      <  3:178177/148‑1 (MQ=255)
                                                                                            cGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTtg                                                      <  6:71466/148‑1 (MQ=255)
                                                                                              aCAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGTTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGa                                                                         <  1:262158/127‑1 (MQ=255)
                                                                                              aCAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTTCACCGACTTGCTGACCTCTTCAGGCGCCTGTtgctg                                                   >  4:179562/1‑149 (MQ=255)
                                                                                              aCAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGa                                                                         >  2:262158/1‑127 (MQ=255)
                                                                                                accgACGATCAATAATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGtt                                                       <  5:98018/140‑1 (MQ=255)
                                                                                                accgACGATCAATAATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGtt                                                       >  6:98018/4‑143 (MQ=255)
                                                                                                  ggACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGAcc                                                                       >  6:354459/1‑125 (MQ=255)
                                                                                                  cgACGATCAATAATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTgcgcg                                               >  8:276657/2‑149 (MQ=255)
                                                                                                   gACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTgcgcgc                                              <  6:217320/149‑1 (MQ=255)
                                                                                                     cGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGccc                                            <  3:19686/149‑1 (MQ=255)
                                                                                                      gATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCt                                                                                           <  8:251039/101‑1 (MQ=255)
                                                                                                              aTCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTc                                    <  4:323280/148‑1 (MQ=255)
                                                                                                              aTCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCa                                   <  8:143745/149‑1 (MQ=255)
                                                                                                               tcgtcgGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTtctc                                       >  1:249987/1‑144 (MQ=255)
                                                                                                               tcgtcgGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTtctc                                       <  2:249987/144‑1 (MQ=255)
                                                                                                                 gtcggtcgGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTGCTACCCCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATcc                                <  7:135121/149‑1 (MQ=255)
                                                                                                                        ggtAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGgtgt                         >  6:167365/1‑149 (MQ=255)
                                                                                                                          tAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCt                       <  5:167365/149‑1 (MQ=255)
                                                                                                                              aCAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCTCCCTTCTCTTCATCCTCGGTGTCTTCAg                   <  5:131099/149‑1 (MQ=255)
                                                                                                                                aGAAAACGCCACCTGCTGAGGTGGCGTTTTCCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAgc                  >  1:34620/1‑149 (MQ=255)
                                                                                                                                aGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTATGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAgcc                 >  1:55887/1‑149 (MQ=255)
                                                                                                                                 gAAAACGCCACCTCAGCAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAgccg                >  7:21224/1‑149 (MQ=255)
                                                                                                                                 gAAAACGCCACCTCAGCAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGAGTCTTCAgccg                >  4:263411/1‑149 (MQ=255)
                                                                                                                                   aaaCGCCACCTCAGCAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAgccgc               >  7:200890/1‑148 (MQ=255)
                                                                                                                                            ctgctgAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAGCCGCCACCGGCGcc     <  7:276657/149‑1 (MQ=255)
                                                                                                                                              gctgAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAGCCGCCACCGGCGCCTg   >  4:7874/1‑149 (MQ=255)
                                                                                                                                               ctgAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCg                             >  2:135798/1‑122 (MQ=255)
                                                                                                                                               ctgAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCg                             <  1:135798/122‑1 (MQ=255)
                                                                                                                                                tgAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTTGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTAAGGCGCCTGTTGCTGCGAGCCGTTCTCTTAATCCTCGGTGTATTCAGCCGCCACCGGCGCCTGc  >  3:109683/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TTTGCCGGCGTGCGTTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAGCCGCCACCGGCGCCTGC  >  NC_002947/4988220‑4988511

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