Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,973,885:1 +T 6.0% coding (1719/2043 nt) flgK ← flagellar hook‑associated protein FlgK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,973,8851.T 6.0% 273.5 / 10.2 67coding (1719/2043 nt)flgKflagellar hook‑associated protein FlgK
Reads supporting (aligned to +/- strand):  ref base . (32/31);  new base T (2/2);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.67e-01

AATGGTGGCCTGGCTGGCCTTGGCCCCTACCCGCTCGACCAGCGAGGCGTAGGCCGAGGTGAAGCTGGTACCGGTGCCGGAATTGGTGCCCACGGTCGACTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATG  >  NC_002947/4973741‑4974029
                                                                                                                                                  |                                                                                                                                                
aaTGGTGGCCTGGCTGGCCTTGGCCCCTACCCGCTCGACCAGCGAGGCGTAGGCCGAGGTGAAGCTGGTACCGGTGCCGGAATTGGTGCCCACGGTCGACTTGGTCTGCAGGTCCAAGCAGCGCGTTGGCATTGCGGTTGTCCGCCtt                                                                                                                                                <  1:42341/148‑1 (MQ=255)
aaTGGTGGCCTGGCTGGCCTTGGCCCCTACCCGCTCGACCAGCGAGGCGTAGGCCGAGGTGAAGCTGGTACCGGTGCCGGAATTGGTGCCCACGGTCGACTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGc                                                                                                                                              >  3:403405/1‑149 (MQ=255)
  tggtggCCTGGCTGGCCTTGGCCCCTACCCGCTCGACCAGCGAGGCGTAGGCCGAGGTGAAGCTGGTACCGGTGCCGGAATTGGTGCCCACGGTCGACTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCg                                                                                                                                            >  3:104269/1‑149 (MQ=255)
       gCCTGGCTGGCCTTGGCCCCTACCCGCTCGACCAGCGAGGCGTAGGCCGAGGTGAAGCTGGTACCGGTGCCGGAATTGGTGCCCACGGTCGACTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCgg                                                                                                                                        >  5:28456/1‑148 (MQ=255)
                  cTTGGCCCCTACCCTCTCGACCAGCGAGGCGTAGGCCGAGGTGAAGCTGGTACCGGTGCCGGAATTGGTGCCCACGGTCGACTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACa                                                                                                                            <  4:310430/149‑1 (MQ=255)
                  cTTGGCCCCTACCCGCTCGACCAGCGAGGCGTAGGCCTAGGTGAAGCTGGTACCGGTGCCGGAATTGGTGCCCACGGTCGACTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACa                                                                                                                            <  7:248842/149‑1 (MQ=255)
                     ggCCCCTACCCGCTCGCCCCGCGAGGCGTAGGCCGAGGTGAAGCTGGTACCGGTGCCGGAATTGGTGCCCACGGTCGACTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAg                                                                                                                           <  4:316616/147‑1 (MQ=255)
                     ggCCCCTACCCGCTCGACCAGCGAGGCGTAGGCCGAGGTGAAGCTGGTACCGGTGCCGGAATTGGTGCCCACGGTCGACTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAgg                                                                                                                          <  8:196509/148‑1 (MQ=255)
                      gCCCCTACCCGCTCGACCAGCGAGGCGTAGGCCGAGGTGAAGCTGGTACCGGTGCCGGAATTGGTGCCCACGGTCGACTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGt                                                                                                                        <  6:28456/149‑1 (MQ=255)
                        cccTACCCGCTCGACCAGCGAGGCGTAGGCCGAGGTGAAGCTGGTACCGGTGCCGGAATTGGTGCCCACGGTCGACTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTa                                                                                                                       >  4:53209/1‑148 (MQ=255)
                        cccTACCCGCTCGACCAGCGAGGCGTAGGCCGAGGTGAAGCTGGTACCGGTGCCGGAATTGGTGCCCACGGTCGACTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTa                                                                                                                       >  3:56090/1‑148 (MQ=255)
                        cccTACCCGCTCGACCAGCGAGGCGTAGGCCGAGGTGAAGCTGGTACCGGTGCCGGAATTGGTGCCCACGGTCGACTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGACCGCCTTGCCGTCGGCATTGAACGACAGGGTa                                                                                                                       >  8:108298/1‑148 (MQ=255)
                         ccTACCCGCTCGACCAGCGAGGCGTAGGCCGAGGTGAAGCTGGTACCGGTGCCGGAATTGGTGCCCACGGTCGACTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTa                                                                                                                       >  5:406970/1‑147 (MQ=255)
                          cTACCGGCTCGACCCGCGAGGAGTAGGCCGAGGTGAAGCTGGTACCGGTGCCGGAATTGGTGCCCACGGTCGACTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAg                                                                                                                    <  8:310764/149‑1 (MQ=255)
                           tACCCGCTCGACCAGCGAGGCGTAGGCCGAGGTGAAGCTGGTACCGGTGCCGGAATTGGTGCCCACGGTCGACTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGc                                                                                                                   >  3:298273/1‑149 (MQ=255)
                              ccGCTCGACCAGCGAGGCGTAGGCCGAGGTGAAGCTGGTACCGGTGCCGGAATTGGTGCCCATGGTCGACTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGt                                                                                                                >  6:348621/1‑149 (MQ=255)
                              ccGCTCGACCAGCGAGGCGTAGGCCGAGGTGAAGCTGGTACCGGTGCCGGAATTGGTGCCCACGGTCGACTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGt                                                                                                                >  7:319451/1‑149 (MQ=255)
                                    gACCAGCGAGGCGTAGGCCGAGGTGAAGCTGGTACCGGTGCCGGAATTGGTGCCCACGGTCGACTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTgg                                                                                                          <  6:406970/149‑1 (MQ=255)
                                         gCGAGGCGTAGGCCGAGGTGAAGCTGGTACCGGTGCCGGAATTGGTGCCCACGGTCGACTTGGTCTGCAGGTCCCGCCGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGCCAGGGTAAAGCTGTCATTGGCCg                                                                                                       >  6:263701/1‑147 (MQ=255)
                                         gCGAGGCGTAGGCCGAGGTGAAGCTGGTACCGGTGCCGGAATTGGTGCCCACGGTCGACTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCa                                                                                                              >  7:20725/1‑140 (MQ=255)
                                         gCGAGGCGTAGGCCGAGGTGAAGCTGGTACCGGTGCCGGAATTGGTGCCCACGGTCGACTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCa                                                                                                              <  8:20725/140‑1 (MQ=255)
                                         gCGAGGCGTAGGCCGAGGTGAAGCTGGTACCGGTGCCGGAATTGGTGCCCACGGTCGACTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCa                                                                                                     <  3:53209/149‑1 (MQ=255)
                                            aGGCGTAGGCCGAGGTGAAGCTGGTACCGGTGCCGGAATTGGTGCCCACGGTCGACTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAg                                                                                                                    <  1:123394/131‑1 (MQ=255)
                                            aGGCGTAGGCCGAGGTGAAGCTGGTACCGGTGCCGGAATTGGTGCCCACGGTCGACTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAg                                                                                                                    >  2:123394/1‑131 (MQ=255)
                                                   ggCCGAGGTGAAGCTGGTACCGGTGCCGGAATTGGTGCCCATGGTCGACTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCa                                                                                            <  5:348621/148‑1 (MQ=255)
                                                       gAGGTGAAGCTGGTACCGGTGCCGGAATTGGTGCCCACGGTCGACTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGcc                                                                                                        >  4:189922/1‑132 (MQ=255)
                                                       gAGGTGAAGCTGGTACCGGTGCCGGAATTGGTGCCCACGGTCGACTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGcc                                                                                                        <  3:189922/132‑1 (MQ=255)
                                                       gAGGTGAAGCTGGTACCGGTGCCGGAATTGGTGCCCACGGTCGACTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGAt                                                                                       <  8:319451/149‑1 (MQ=255)
                                                          gTGAAGCTGGTACCGGTGCCGGAATTGGTGCCCACGGTCGACTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATtg                                                                                     <  4:403405/148‑1 (MQ=255)
                                                           tGAAGCTGGTACCGGTGCCGGAATTGGTGCCCACGGTCGACTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATtgtg                                                                                   <  5:260816/149‑1 (MQ=255)
                                                           tGAAGCTGGTACCGGTGCCGGAATTGGTGCCCACGGTCGACTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGACCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGTCACCGATtgt                                                                                    <  7:108298/148‑1 (MQ=255)
                                                                    tACCGGTGCCGGAATTGGTGCCCACGGTCGACTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAggt                                                                                                                  >  5:356243/1‑107 (MQ=255)
                                                                    tACCGGTGCCGGAATTGGTGCCCACGGTCGACTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCt                                                                                                                  <  6:356243/109‑1 (MQ=255)
                                                                         gTGCCGGAATTGGTGCCCACGGTCGACTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATtgt                                                                                    <  4:73864/134‑1 (MQ=255)
                                                                         gTGCCGGAATTGGTGCCCACGGTCGACTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATtgt                                                                                    >  3:73864/1‑134 (MQ=255)
                                                                                  ttGGTGCCCACGGTCGACTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTGCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCgg                                                            <  7:127684/149‑1 (MQ=255)
                                                                                     gTGCCCACGGTCGACTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGt                                                          >  2:268727/1‑149 (MQ=255)
                                                                                     gTGCCCACGGTCGACTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGt                                                          >  4:193672/1‑149 (MQ=255)
                                                                                      tGCCCACGGTCGACTTGGTCTGCAGGTCCATCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCgg                                                            <  2:57666/145‑1 (MQ=255)
                                                                                      tGCCCACGGTCGACTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCgg                                                            >  1:57666/1‑145 (MQ=255)
                                                                                         ccACGGTCGACTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCg                                                     >  5:379325/1‑149 (MQ=255)
                                                                                          cACGGTCGACTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGc                                                    >  3:236643/1‑149 (MQ=255)
                                                                                          cACGGTCGACTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGc                                                    >  6:251703/1‑149 (MQ=255)
                                                                                               tCGACTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCa                                               >  6:368008/1‑149 (MQ=255)
                                                                                                  aCTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATc                                             <  1:268727/149‑1 (MQ=255)
                                                                                                  aCTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATc                                             <  3:193672/149‑1 (MQ=255)
                                                                                                   cTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGGAAAGGTACGCGGGTTACCGCTAcc                                                >  5:371311/1‑144 (MQ=255)
                                                                                                   cTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTAcc                                                <  6:371311/144‑1 (MQ=255)
                                                                                                    ttGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATc                                             >  7:175526/1‑146 (MQ=255)
                                                                                                            cAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGt                                  >  5:36172/1‑149 (MQ=255)
                                                                                                               gTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAg                                                                    <  3:124401/112‑1 (MQ=255)
                                                                                                               gTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAg                                                                    >  4:124401/1‑112 (MQ=255)
                                                                                                                    agcagcGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCgccggc                           <  5:251703/149‑1 (MQ=255)
                                                                                                                          gcgTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATc                     >  7:218778/1‑149 (MQ=255)
                                                                                                                           cgTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATcc                    <  6:36172/149‑1 (MQ=255)
                                                                                                                            gTTGGTATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATc                     >  2:88781/1‑147 (MQ=255)
                                                                                                                            gTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATc                     >  8:134661/1‑147 (MQ=255)
                                                                                                                                gCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTATAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCgg                                    >  1:288555/1‑128 (MQ=255)
                                                                                                                                gCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCgg                                    <  2:288555/128‑1 (MQ=255)
                                                                                                                                 cATTGCTGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATc              <  4:298273/149‑1 (MQ=255)
                                                                                                                                  aTTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCg             >  4:196198/1‑149 (MQ=255)
                                                                                                                                   ttGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCgg            <  8:175526/149‑1 (MQ=255)
                                                                                                                                      cGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCAc         <  4:236643/149‑1 (MQ=255)
                                                                                                                                       ggTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACg        <  1:88781/149‑1 (MQ=255)
                                                                                                                                       ggTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTATCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACg        >  1:79969/1‑149 (MQ=255)
                                                                                                                                        gTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACgg       <  8:218778/149‑1 (MQ=255)
                                                                                                                                            tCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGAt   <  1:57815/149‑1 (MQ=255)
                                                                                                                                             ccGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATg  >  7:77236/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
AATGGTGGCCTGGCTGGCCTTGGCCCCTACCCGCTCGACCAGCGAGGCGTAGGCCGAGGTGAAGCTGGTACCGGTGCCGGAATTGGTGCCCACGGTCGACTTGGTCTGCAGGTCCAGCAGCGCGTTGGCATTGCGGTTGTCCGCCTTGCCGTCGGCATTGAACGACAGGGTAAAGCTGTCATTGGCCGCAGGGCTGCCACCGATTGTGGTCTCGACGGCAAAGGTACGCGGGTTACCGCTACCATCCAGAATCGGGTCGCCGGCCGCATCCAGCATCGGCACGGCGATG  >  NC_002947/4973741‑4974029

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: