Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,991,045 C→G 5.3% P97R (CCG→CGG)  flgA → flagella basal body P‑ring formation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,991,0450CG 5.3% 1033.8 / 34.0 303P97R (CCG→CGG) flgAflagella basal body P‑ring formation protein
Reads supporting (aligned to +/- strand):  ref base C (181/106);  new base G (11/5);  total (192/111)
Fisher's exact test for biased strand distribution p-value = 7.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.23e-01

Only 100 of 304 total aligned reads displayed.

TTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGTGACGTGGCC  >  NC_002947/4990900‑4991193
                                                                                                                                                 |                                                                                                                                                    
tttCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTg                                                                                                                                                  <  1:3203143/150‑1 (MQ=255)
    ttGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGcc                                                                                                                                              >  2:207292/1‑150 (MQ=255)
    ttGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGcc                                                                                                                                              >  2:515528/1‑150 (MQ=255)
     tGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCg                                                                                                                                             >  1:1202543/1‑150 (MQ=255)
      gAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCGGCTGGGCCgg                                                                                                                                            <  2:4122896/150‑2 (MQ=255)
       aaTTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCGGCTGGGCCggg                                                                                                                                           <  2:4275838/150‑3 (MQ=255)
         ttCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCTGCCAGCAGCTGGAGGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGtgtg                                                                                                                                         >  2:1724702/1‑150 (MQ=255)
                ttGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTAc                                                                                                                                  >  2:5418196/1‑150 (MQ=255)
                ttGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTAc                                                                                                                                  >  2:2388574/1‑150 (MQ=255)
                 tGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACg                                                                                                                                 <  2:71196/150‑1 (MQ=255)
                  gAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACgg                                                                                                                                >  1:165662/1‑150 (MQ=255)
                  gAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACgg                                                                                                                                >  1:165791/1‑150 (MQ=255)
                     gATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATa                                                                                                                             >  1:2053088/1‑150 (MQ=255)
                     gATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATa                                                                                                                             >  2:4912265/1‑150 (MQ=255)
                         aTCTGTCCACCACCCAGACTGCCGGCCGCGATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGcg                                                                                                                         <  1:5008122/150‑2 (MQ=255)
                           cTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGc                                                                                                                       <  2:535689/150‑1 (MQ=255)
                                caccacCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACgg                                                                                                                  >  1:1366669/1‑150 (MQ=255)
                                caccacCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACgg                                                                                                                  >  1:4035816/1‑150 (MQ=255)
                                   caccCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGgcgc                                                                                                               <  1:5463287/150‑1 (MQ=255)
                                    accCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCgca                                                                                                              >  2:3036980/1‑150 (MQ=255)
                                      ccAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCgcagc                                                                                                            >  1:2072417/1‑150 (MQ=255)
                                      ccAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGGGTGACGGTACGGATACGCTGCACCGGCgcagc                                                                                                            >  1:1035885/1‑150 (MQ=255)
                                       cAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCgcagca                                                                                                           >  1:2469080/1‑150 (MQ=255)
                                         gACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCAcc                                                                                                         >  2:19050/1‑150 (MQ=255)
                                                ggccgCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACg                                                                                                  >  2:294781/1‑150 (MQ=255)
                                                 gccgCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACgg                                                                                                 <  2:3675485/150‑1 (MQ=255)
                                                  ccgCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGt                                                                                                >  2:586043/1‑150 (MQ=255)
                                                   cgCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTa                                                                                               <  1:2937548/150‑1 (MQ=255)
                                                     cTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTAtt                                                                                             <  2:5249417/150‑1 (MQ=255)
                                                      tATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTCTTc                                                                                            >  2:2491190/1‑150 (MQ=255)
                                                      tATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTc                                                                                            >  2:4963287/1‑150 (MQ=255)
                                                      tATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTc                                                                                            >  1:491804/1‑150 (MQ=255)
                                                      tATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTc                                                                                            <  1:685386/150‑1 (MQ=255)
                                                      tATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTc                                                                                            >  2:3887145/1‑150 (MQ=255)
                                                      gaTGAAATCCCGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTc                                                                                            <  1:268376/149‑1 (MQ=255)
                                                       aTGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCg                                                                                           >  2:2287677/1‑150 (MQ=255)
                                                       aTGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCg                                                                                           <  2:1194929/150‑1 (MQ=255)
                                                        tGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGt                                                                                          >  2:3421780/1‑150 (MQ=255)
                                                          aaaTCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGc                                                                                        >  2:1044250/1‑150 (MQ=255)
                                                          aaaTCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGc                                                                                        <  1:5097235/150‑1 (MQ=255)
                                                           aaTCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGcc                                                                                       <  1:4176607/150‑1 (MQ=255)
                                                             tCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCgg                                                                                     >  2:5596052/1‑150 (MQ=255)
                                                              ccAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGc                                                                                    >  1:3619989/1‑150 (MQ=255)
                                                                aGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCa                                                                                  >  2:2931291/1‑150 (MQ=255)
                                                                  gTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCAcc                                                                                >  1:1128471/1‑150 (MQ=255)
                                                                      aCCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGc                                                                            >  1:3896227/1‑150 (MQ=255)
                                                                       cccGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCg                                                                           >  1:5299253/1‑150 (MQ=255)
                                                                       cccGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCg                                                                           >  1:4671460/1‑150 (MQ=255)
                                                                       cccGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCg                                                                           >  1:2612528/1‑150 (MQ=255)
                                                                       cccGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCg                                                                           >  2:1277447/1‑150 (MQ=255)
                                                                         cGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGAc                                                                         >  2:2665106/1‑150 (MQ=255)
                                                                               aCCCACGCCTGCGCATGCCATTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCAGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTcc                                                                   <  2:3830996/150‑1 (MQ=255)
                                                                                 ccACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCcgc                                                                 <  2:421373/150‑1 (MQ=255)
                                                                                   aCGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGa                                                               >  1:984054/1‑150 (MQ=255)
                                                                                     gCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACg                                                             >  1:3223095/1‑150 (MQ=255)
                                                                                             aTGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACgtggtggtg                                                     >  1:5621953/1‑150 (MQ=255)
                                                                                               gCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACgtggtggtggt                                                   <  2:3955674/150‑1 (MQ=255)
                                                                                                 cGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGa                                                 >  1:1167703/1‑150 (MQ=255)
                                                                                                     gtgCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGAcccg                                             <  1:4484831/150‑1 (MQ=255)
                                                                                                             agcagcTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCa                                     >  1:1540124/1‑150 (MQ=255)
                                                                                                             agcagcTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCa                                     >  1:1252136/1‑150 (MQ=255)
                                                                                                             agcagcTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCa                                     >  1:3486452/1‑150 (MQ=255)
                                                                                                             agcagcTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCa                                     >  2:3278532/1‑150 (MQ=255)
                                                                                                              gcagcTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCaa                                    >  2:434761/1‑150 (MQ=255)
                                                                                                                agcTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAgc                                  <  2:1583139/150‑1 (MQ=255)
                                                                                                                agcTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAgc                                  <  1:4156121/150‑1 (MQ=255)
                                                                                                                agcTGGACGCCTAACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTg                                                                                         >  1:2812453/1‑95 (MQ=255)
                                                                                                                     gACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGac                             >  2:1465482/1‑150 (MQ=255)
                                                                                                                      aCGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGaca                            <  1:552511/150‑1 (MQ=255)
                                                                                                                        gCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGacatc                          <  2:4799713/150‑3 (MQ=255)
                                                                                                                        gCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGacaac                          >  1:194987/1‑150 (MQ=255)
                                                                                                                          cTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAacac                        >  2:2378235/1‑150 (MQ=255)
                                                                                                                           tCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTCTTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACACCACg                       >  2:5065295/1‑150 (MQ=255)
                                                                                                                           tCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACg                       <  2:2838102/150‑1 (MQ=255)
                                                                                                                             aCTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGt                     <  1:4807993/150‑1 (MQ=255)
                                                                                                                                 gAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGc                 >  1:2123243/1‑150 (MQ=255)
                                                                                                                                 gAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGc                 <  1:3977186/150‑1 (MQ=255)
                                                                                                                                   aaGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTc                    >  1:1711857/1‑145 (MQ=255)
                                                                                                                                   aaGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGa               <  2:4834365/150‑1 (MQ=255)
                                                                                                                                     gCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAg             <  1:3417875/150‑1 (MQ=255)
                                                                                                                                     gCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAg             <  1:2885408/150‑1 (MQ=255)
                                                                                                                                     acacGGCCCAGCGGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAg             >  2:2804983/4‑150 (MQ=255)
                                                                                                                                     acacGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAg             >  2:5123217/4‑150 (MQ=255)
                                                                                                                                      cccGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAgg            >  1:1600954/1‑150 (MQ=255)
                                                                                                                                      cccGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAgg            >  2:785257/1‑150 (MQ=255)
                                                                                                                                      cccGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAgg            >  1:2176733/1‑150 (MQ=255)
                                                                                                                                      cccGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAgg            >  2:3134406/1‑150 (MQ=255)
                                                                                                                                      cccGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAgg            >  1:3557810/1‑150 (MQ=255)
                                                                                                                                       ccGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGt           >  1:3040488/1‑150 (MQ=255)
                                                                                                                                       ccGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGt           <  2:1607661/150‑1 (MQ=255)
                                                                                                                                       ccGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGt           >  1:799282/1‑150 (MQ=255)
                                                                                                                                       acGGCCCAGCGGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGt           >  1:2956377/2‑150 (MQ=255)
                                                                                                                                       acGGCCCAGCGGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGt           >  1:4597229/2‑150 (MQ=255)
                                                                                                                                        cGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGTg          >  1:527331/1‑150 (MQ=255)
                                                                                                                                          gCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGTGAc        <  2:5291475/150‑1 (MQ=255)
                                                                                                                                             cAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGTGACGTg     >  1:4430031/1‑150 (MQ=255)
                                                                                                                                             cAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGTGACGTg     >  2:1358297/1‑150 (MQ=255)
                                                                                                                                               gCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGTGACGTGGc   <  2:4685881/150‑1 (MQ=255)
                                                                                                                                               gCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGTGACGTGGc   >  1:1920612/1‑150 (MQ=255)
                                                                                                                                                ccGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGTGACGTGGcc  >  1:1821328/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGTGACGTGGCC  >  NC_002947/4990900‑4991193

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