Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,988,018 C→G 6.9% 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,988,0180CG 6.9% 697.1 / 31.9 204P97R (CCG→CGG) flgAflagella basal body P‑ring formation protein
Reads supporting (aligned to +/- strand):  ref base C (118/72);  new base G (2/12);  total (120/84)
Fisher's exact test for biased strand distribution p-value = 9.89e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 206 total aligned reads displayed.

GGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGTGACGTGGCC  >  NC_002947/4987871‑4988166
                                                                                                                                                   |                                                                                                                                                    
 gTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCt                                                                                                                                                   <  2:3008102/150‑1 (MQ=255)
 gTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCACGGCCCAGCGGCt                                                                                                                                                   <  2:140/150‑1 (MQ=255)
  tttCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTg                                                                                                                                                  <  1:3100169/150‑1 (MQ=255)
  tttCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTg                                                                                                                                                  <  2:3065722/150‑1 (MQ=255)
    tCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTggg                                                                                                                                                >  1:23779/1‑150 (MQ=255)
    tCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTggg                                                                                                                                                >  1:1618068/1‑150 (MQ=255)
        gAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCGGCTGGGCCgg                                                                                                                                            <  2:3814984/150‑2 (MQ=255)
        gAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCgt                                                                                                                                            >  2:1120747/1‑150 (MQ=255)
         aaTTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCGGCTGGGCCggg                                                                                                                                           <  1:4495938/150‑3 (MQ=255)
          aTTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCGGCTGGGCCgggc                                                                                                                                          <  1:2940315/150‑4 (MQ=255)
          aTTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCACGGCCCAGCGGCTGGGCCgggc                                                                                                                                          <  1:2709582/150‑4 (MQ=255)
           ttCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGAGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGtgtg                                                                                                                                         >  2:4417319/1‑150 (MQ=255)
            tCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGa                                                                                                                                        >  2:2060013/1‑150 (MQ=255)
              aCCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACg                                                                                                                                      >  1:3479359/1‑150 (MQ=255)
              aCCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCACGGCCCAGCGGCTGGGCCgggctttc                                                                                                                                      <  1:1825544/150‑8 (MQ=255)
                 gTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTa                                                                                                                                   >  2:1465458/1‑150 (MQ=255)
                  ttGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTAc                                                                                                                                  >  1:1415819/1‑150 (MQ=255)
                   tGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACg                                                                                                                                 >  2:3254218/1‑150 (MQ=255)
                     aaGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGa                                                                                                                               <  1:4405906/150‑1 (MQ=255)
                       gATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGTCGTGTGACGGTACGGATa                                                                                                                             >  1:3725689/1‑150 (MQ=255)
                       gATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATa                                                                                                                             >  1:2309727/1‑150 (MQ=255)
                       gATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATa                                                                                                                             >  1:3745301/1‑150 (MQ=255)
                       gATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATa                                                                                                                             <  2:1226998/150‑1 (MQ=255)
                             cTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGc                                                                                                                       >  1:3289544/1‑150 (MQ=255)
                               gTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCaa                                                                                                                     <  1:732512/150‑1 (MQ=255)
                               gTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCaa                                                                                                                     >  2:3328435/1‑150 (MQ=255)
                               gTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCaa                                                                                                                     >  1:1963035/1‑150 (MQ=255)
                                  caccacCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACgg                                                                                                                  >  1:936159/1‑150 (MQ=255)
                                      accCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACTCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCgca                                                                                                              >  2:3842816/1‑150 (MQ=255)
                                       cccAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCgcag                                                                                                             <  1:1064289/150‑1 (MQ=255)
                                        ccAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCgcagc                                                                                                            >  2:3920077/1‑150 (MQ=255)
                                        ccAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCgcagc                                                                                                            >  1:3114200/1‑150 (MQ=255)
                                         cAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCgcagca                                                                                                           >  1:203382/1‑150 (MQ=255)
                                             cTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGt                                                                                                       <  1:4468691/150‑1 (MQ=255)
                                             cTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGt                                                                                                       <  2:4507867/150‑1 (MQ=255)
                                                  ggccgCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACg                                                                                                  >  2:1732092/1‑150 (MQ=255)
                                                  ggccgCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACg                                                                                                  <  1:1111867/150‑1 (MQ=255)
                                                    ccgCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGCCGGt                                                                                                >  1:245485/1‑150 (MQ=255)
                                                    ccgCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGt                                                                                                >  1:2130800/1‑150 (MQ=255)
                                                      gCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTAt                                                                                              >  1:3022061/1‑150 (MQ=255)
                                                        tATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTc                                                                                            >  1:2523208/1‑150 (MQ=255)
                                                        tATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTc                                                                                            <  1:2414022/150‑1 (MQ=255)
                                                            aaaTCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGc                                                                                        <  1:3412662/150‑1 (MQ=255)
                                                              aTCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCg                                                                                      >  1:4337107/1‑150 (MQ=255)
                                                               tCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCgg                                                                                     >  2:48415/1‑150 (MQ=255)
                                                               tCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCgg                                                                                     >  2:2456415/1‑150 (MQ=255)
                                                                ccAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGc                                                                                    >  1:2408389/1‑150 (MQ=255)
                                                                ccAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGc                                                                                    <  2:2197329/150‑1 (MQ=255)
                                                                 cAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGTAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGcc                                                                                   >  1:3334557/1‑150 (MQ=255)
                                                                       aaCCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTg                                                                             <  1:2956692/150‑1 (MQ=255)
                                                                       aaCCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTg                                                                             >  1:1758612/1‑150 (MQ=255)
                                                                            gCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACt                                                                        >  1:3344445/1‑150 (MQ=255)
                                                                              tGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACttt                                                                      >  1:453603/1‑150 (MQ=255)
                                                                                 aCCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTcc                                                                   <  1:3685340/150‑1 (MQ=255)
                                                                                   ccACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCcgc                                                                 <  1:1517378/150‑1 (MQ=255)
                                                                                    cACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCcgcg                                                                >  2:87145/1‑150 (MQ=255)
                                                                                          tGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACgtgg                                                          <  2:1063360/150‑1 (MQ=255)
                                                                                          tGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACgtgg                                                          >  1:173457/1‑150 (MQ=255)
                                                                                             gcATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTTCGCGACgtggtgg                                                       >  1:638562/1‑150 (MQ=255)
                                                                                              cATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACgtggtggt                                                      <  1:3914346/150‑1 (MQ=255)
                                                                                               aTGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACgtggtggtg                                                     >  2:4136206/1‑150 (MQ=255)
                                                                                                 gCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACgtgggggtggt                                                   >  1:1085251/1‑150 (MQ=255)
                                                                                                   cGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGa                                                 <  1:3071237/150‑1 (MQ=255)
                                                                                                         gCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGAcccgcc                                           <  2:3619262/150‑1 (MQ=255)
                                                                                                          cagccagcAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGAcccgccc                                          <  1:724593/150‑1 (MQ=255)
                                                                                                          cagccagcAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGAcccgccc                                          <  1:3181015/150‑1 (MQ=255)
                                                                                                          cagccagcAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGAcccgccc                                          <  2:2074759/150‑1 (MQ=255)
                                                                                                             ccagcAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCt                                       <  1:3842816/150‑1 (MQ=255)
                                                                                                             ccagcAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCt                                       >  2:1670114/1‑150 (MQ=255)
                                                                                                                gcagcTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCaa                                    >  2:2287444/1‑150 (MQ=255)
                                                                                                                    cTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAgcgc                                <  1:2238396/150‑1 (MQ=255)
                                                                                                                     tGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGcgcg                               >  2:2002837/1‑150 (MQ=255)
                                                                                                                          gCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGacaac                          >  1:2650081/1‑150 (MQ=255)
                                                                                                                            cTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAacac                        >  1:1910766/1‑150 (MQ=255)
                                                                                                                             tCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACg                       >  1:2583759/1‑150 (MQ=255)
                                                                                                                             tCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACg                       >  2:3525477/1‑150 (MQ=255)
                                                                                                                             tCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACg                       >  1:2984836/1‑150 (MQ=255)
                                                                                                                             tCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACg                       >  2:2728927/1‑150 (MQ=255)
                                                                                                                             tCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACg                       <  1:945566/150‑1 (MQ=255)
                                                                                                                             tCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACg                       >  2:2204637/1‑150 (MQ=255)
                                                                                                                                cTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTc                    <  1:1887538/150‑1 (MQ=255)
                                                                                                                                cTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTc                    <  2:2130081/150‑1 (MQ=255)
                                                                                                                                cTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTc                    >  2:3108833/1‑150 (MQ=255)
                                                                                                                                cTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTc                    <  1:1259837/150‑1 (MQ=255)
                                                                                                                                   gAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGc                 <  2:3808478/150‑1 (MQ=255)
                                                                                                                                       gCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAg             >  2:2558148/1‑150 (MQ=255)
                                                                                                                                        cccGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAgg            <  1:48415/150‑1 (MQ=255)
                                                                                                                                        cccGGCCCAGCCGCCGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAgg            >  1:764075/1‑150 (MQ=255)
                                                                                                                                        cacGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAgg            >  1:1656251/3‑150 (MQ=255)
                                                                                                                                         ccGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGt           <  1:1048129/150‑1 (MQ=255)
                                                                                                                                         ccGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGt           >  1:2806301/1‑150 (MQ=255)
                                                                                                                                              ccAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGTGACGt      >  2:4425424/1‑150 (MQ=255)
                                                                                                                                              ccAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGTGACGt      >  1:1317888/1‑150 (MQ=255)
                                                                                                                                              ccAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGTGACGt      >  2:2009383/1‑150 (MQ=255)
                                                                                                                                               cAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGTGACAACACGGTCGGCGAAGGTGACGTg     >  2:4282936/1‑150 (MQ=255)
                                                                                                                                               cAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGTGACGTg     >  1:585902/1‑150 (MQ=255)
                                                                                                                                               cAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGTGACGTg     >  2:1675283/1‑150 (MQ=255)
                                                                                                                                                 gCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGTGACGTGGc   >  2:4088166/1‑150 (MQ=255)
                                                                                                                                                  ccGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGTGACGTGGcc  <  1:1520974/150‑1 (MQ=255)
                                                                                                                                                  ccGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGTGACGTGGcc  <  1:2275894/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GGTTTCTTGAATTCACCGTTGAAGATTATCTGTCCACCACCCAGACTGCCGGCCGCTATGAAATCCAGGTCAACCCGCTGGACCCACGCCTGCGCATGCCGTTGTGCAGCCAGCAGCTGGACGCCTCACTGGAAAGCCCGGCCCAGCCGCTGGGCCGTGTGACGGTACGGATACGCTGCAACGGCGCAGCACCGTGGACGGTATTCGTGCCGGCCACCGTGCGACTTTTCCGCGACGTGGTGGTGGTGACCCGCCCGCTCAAGCGCGACAACACGGTCGGCGAAGGTGACGTGGCC  >  NC_002947/4987871‑4988166

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: