Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,450,173 C→G 9.7% A68G (GCC→GGC)  nicA → nicotinate dehydrogenase subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,450,1730CG 9.7% 438.7 / 32.8 147A68G (GCC→GGC) nicAnicotinate dehydrogenase subunit A
Reads supporting (aligned to +/- strand):  ref base C (59/73);  new base G (6/8);  total (66/81)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.89e-01

Only 100 of 155 total aligned reads displayed.

CATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATCGCA  >  NC_002947/4450024‑4450322
                                                                                                                                                     |                                                                                                                                                      
    ccTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCgccgcc                                                                                                                                                    <  2:2915285/150‑1 (MQ=255)
      tGACACCCGGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCAAGCTGCATCCCCTTGGCCGGCGCCgccgg                                                                                                                                                  >  1:2008745/1‑150 (MQ=255)
      tGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCgccgg                                                                                                                                                  >  2:1937435/1‑150 (MQ=255)
        acacCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCgccggcc                                                                                                                                                <  1:394711/150‑1 (MQ=255)
        acacCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCgccggcc                                                                                                                                                >  1:1868634/1‑150 (MQ=255)
        acacCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCgccggcc                                                                                                                                                <  1:1429354/150‑1 (MQ=255)
          acCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGGCgccggccaa                                                                                                                                              <  2:142572/150‑3 (MQ=255)
            cccGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGGCgccggccaagg                                                                                                                                            <  2:2571384/150‑5 (MQ=255)
            cccGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGCCATCCCCTTGGCCGGCGCCGCCCGCCGCaa                                                                                                                                            >  1:283998/1‑150 (MQ=255)
             cggctggtgctGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAat                                                                                                                                           >  2:2295230/3‑150 (MQ=255)
             ccGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGGCgccggccaaggg                                                                                                                                           <  2:1457891/150‑6 (MQ=255)
               gctgctgctGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAatat                                                                                                                                         >  2:620349/1‑150 (MQ=255)
                 tgctgctgATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGGCgccggccaag                                                                                                                                             >  2:1977052/1‑141 (MQ=255)
                   ctgctgATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATcacc                                                                                                                                     <  1:2526319/150‑1 (MQ=255)
                    tgctgATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATcacca                                                                                                                                    <  2:466136/150‑1 (MQ=255)
                          tCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCTGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCg                                                                                                                              >  2:2664250/1‑150 (MQ=255)
                          tCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCg                                                                                                                              >  2:795969/1‑150 (MQ=255)
                            cTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGaa                                                                                                                            <  1:1199492/150‑1 (MQ=255)
                            cTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGACCGCCAATATTTCGTCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACACTCGaa                                                                                                                            >  1:188498/1‑150 (MQ=255)
                             tGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGGCgccggccaag                                                                                                                                             <  2:2144786/132‑4 (MQ=255)
                             tGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGGCgccggccaag                                                                                                                                             >  1:2144786/1‑129 (MQ=255)
                             tGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAg                                                                                                                           <  2:1212613/150‑1 (MQ=255)
                                  aaCGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGNCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTg                                                                                                                      <  1:2295230/150‑1 (MQ=255)
                                   aCGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTgg                                                                                                                     <  1:1996158/150‑1 (MQ=255)
                                     gACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGgc                                                                                                                   <  1:2426689/150‑1 (MQ=255)
                                      aCCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGgca                                                                                                                  <  2:2618706/150‑1 (MQ=255)
                                      aCCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGgca                                                                                                                  <  2:1984067/150‑1 (MQ=255)
                                      aCCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGgca                                                                                                                  >  2:1523033/1‑150 (MQ=255)
                                       ccTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGgcag                                                                                                                 >  1:2601459/1‑150 (MQ=255)
                                        cTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGgcagc                                                                                                                >  2:2694597/1‑150 (MQ=255)
                                        cTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGgcagc                                                                                                                >  2:1787992/1‑150 (MQ=255)
                                         tgtgCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGgcagca                                                                                                               >  2:635502/1‑150 (MQ=255)
                                            gCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAgg                                                                                                            <  2:2601459/150‑1 (MQ=255)
                                             ccTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCaa                                                                                                              >  1:41761/1‑147 (MQ=255)
                                             ccTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCaa                                                                                                              <  2:41761/147‑1 (MQ=255)
                                              cTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGcc                                                                                                          <  1:2640865/150‑1 (MQ=255)
                                               tGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCaccacc                                                                                                                                  >  1:987263/1‑125 (MQ=255)
                                               tGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCaccacc                                                                                                                                  <  2:987263/125‑1 (MQ=255)
                                               tGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGccg                                                                                                         >  2:2921702/1‑150 (MQ=255)
                                                gAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATcac                                                                                                                                      <  1:2792897/120‑1 (MQ=255)
                                                gAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATcac                                                                                                                                      >  2:2792897/1‑120 (MQ=255)
                                                  aCGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGccgcac                                                                                                     >  1:2565124/1‑150 (MQ=255)
                                                      ccccAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGccgcaccgca                                                                                                 <  2:2530801/150‑1 (MQ=255)
                                                        ccAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGgca                                                                                                                  >  2:1879725/1‑132 (MQ=255)
                                                        ccAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTCGGCCGGCGGCgccggccaag                                                                                                                                             <  2:1037103/105‑4 (MQ=255)
                                                        ccAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTCGGCCGGCGGCgccggccaag                                                                                                                                             >  1:1037103/1‑102 (MQ=255)
                                                            gTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGt                                                                                           <  1:1937435/150‑1 (MQ=255)
                                                             tACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCCCCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTg                                                                                          <  1:544287/150‑1 (MQ=255)
                                                                ggCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCaa                                                                                       <  1:795969/150‑1 (MQ=255)
                                                                     cGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGc                                                                                  <  1:1787992/150‑1 (MQ=255)
                                                                     cGGCCGGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGACCCTGGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGc                                                                                  <  2:106323/150‑1 (MQ=255)
                                                                         cTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTc                                                                              <  2:1076246/150‑1 (MQ=255)
                                                                          tGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCa                                                                             >  2:1181017/1‑150 (MQ=255)
                                                                          tGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCa                                                                             >  1:1040623/1‑150 (MQ=255)
                                                                           gggCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGcc                                                                                                          <  1:971622/121‑1 (MQ=255)
                                                                           gggCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCAt                                                                            <  1:501606/150‑1 (MQ=255)
                                                                            ggCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGTCGGCTCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATc                                                                           <  1:2359650/150‑1 (MQ=255)
                                                                            ggCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATc                                                                           >  2:2142004/1‑150 (MQ=255)
                                                                              cGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATcga                                                                         <  1:3054318/150‑1 (MQ=255)
                                                                                 aTGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATcgacga                                                                      >  1:1687493/1‑150 (MQ=255)
                                                                                     ggggCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATCTCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAg                                                                  <  1:1523033/150‑1 (MQ=255)
                                                                                      gggCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAgg                                                                 >  1:1497341/1‑150 (MQ=255)
                                                                                       ggCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCATCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAggc                                                                <  2:742735/150‑1 (MQ=255)
                                                                                        gCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTGATCCCCTTCGCCGGCGCCGCCGGCCGCCATATCACCACCAGCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAggcg                                                               <  1:581589/150‑1 (MQ=255)
                                                                                          gtgtACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGccgca                                                                                                      <  2:521103/109‑1 (MQ=255)
                                                                                               cGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAgt                                                        <  2:1254436/150‑1 (MQ=255)
                                                                                               cGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAgt                                                        <  1:3034901/150‑1 (MQ=255)
                                                                                               cGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAgt                                                        >  2:3034901/1‑150 (MQ=255)
                                                                                               cGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCCTCATCCCGTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAgt                                                        <  2:1111926/150‑1 (MQ=255)
                                                                                                ggTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGgcag                                                                                                                 <  1:155451/93‑1 (MQ=255)
                                                                                                  tGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAg                                                                  >  2:144042/1‑137 (MQ=255)
                                                                                                  tGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAg                                                                  <  1:144042/137‑1 (MQ=255)
                                                                                                   tctcatcGACGGCGTGGCCACCCGCTCCTGCGTCATCCCCTTGTCCGGCGCCGCTTGCCGCTATATCACCACCCTCGAAGGGCTTGGCAGCAAGGCCGCACCGCCCCCCGTGCAAGAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTgg                                                    <  2:1195587/148‑1 (MQ=255)
                                                                                                    atcatcGACGGCGTGGCCGCCCGCTCCTGCGTCAAAAGCGTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGt                                                   <  1:2464059/150‑1 (MQ=255)
                                                                                                     tcatcGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGtt                                                  >  1:1296062/1‑150 (MQ=255)
                                                                                                      catcGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCACACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTa                                                 <  2:1541440/150‑1 (MQ=255)
                                                                                                       atcGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTAc                                                >  2:1864319/1‑150 (MQ=255)
                                                                                                              gCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATcgac                                                                        <  1:197659/119‑1 (MQ=255)
                                                                                                                 tGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACg                                      <  1:1220687/150‑1 (MQ=255)
                                                                                                                   gccgccCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCCAGGCCGACGAAGGCCTGTGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGc                                    <  1:2664250/150‑1 (MQ=255)
                                                                                                                     cgcccgcTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCAt                                  >  1:2536118/1‑150 (MQ=255)
                                                                                                                       cccgcTCCTGCGTCATCCCCTGGGCCGGCGCCGCCGGCCGCAATATCCCCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGa                                <  1:1195442/150‑1 (MQ=255)
                                                                                                                         cgcTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCAt                                  <  2:1478903/146‑1 (MQ=255)
                                                                                                                              cTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATg                           <  1:1821436/148‑1 (MQ=255)
                                                                                                                              cTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATg                           >  2:1821436/1‑148 (MQ=255)
                                                                                                                                     aTCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAg                  <  1:635502/150‑1 (MQ=255)
                                                                                                                                     aTCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAg                  >  1:499499/1‑150 (MQ=255)
                                                                                                                                     aTCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAg                  <  2:499499/150‑1 (MQ=255)
                                                                                                                                       ccccTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGccgcacc                                                                                                    <  1:137902/66‑1 (MQ=255)
                                                                                                                                       ccccTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGccgcacc                                                                                                    >  2:137902/1‑66 (MQ=255)
                                                                                                                                         ccTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGccc              >  1:895081/1‑150 (MQ=255)
                                                                                                                                            cggCCGGCGGCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATcg                                                                          <  1:1805088/86‑1 (MQ=255)
                                                                                                                                             ggCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCt          <  1:2189507/150‑1 (MQ=255)
                                                                                                                                                cGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGa       <  2:1040623/150‑1 (MQ=255)
                                                                                                                                                 ggCGCCGCCGGCCGCAATATCACCACCCTCGGAGGGCTGGGCAGCAAGGCCGCACCGCACCCCCTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGTCCCTGCTCGAt      <  2:283998/150‑1 (MQ=255)
                                                                                                                                                  gcgcCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATc     <  1:2694597/150‑1 (MQ=255)
                                                                                                                                                   ccccgccgGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGAAAAGCTGTCACCAAAAATACTCATCCACAGTTCAGCGACGAGCATTCGCCCCAGTGTGGTTACTGCATGAACGTCAGGATCATGACCGCCAAGGCCCTGCTCGATCg    <  1:1465617/148‑1 (MQ=255)
                                                                                                                                                     ccgccgGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAAc                                       >  2:1213459/1‑113 (MQ=255)
                                                                                                                                                     ccgccgGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAAc                                       <  1:1213459/113‑1 (MQ=255)
                                                                                                                                                     ccgccgGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATCGCa  <  1:2962016/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CATGCCTGACACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATCGCA  >  NC_002947/4450024‑4450322

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