Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,005,053 A→T 10.8% intergenic (+35/‑141) PP_5259 → / → ydcJ LysR family transcriptional regulator/metalloprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,005,0530AT10.8% 285.3 / 19.0 100intergenic (+35/‑141)PP_5259/ydcJLysR family transcriptional regulator/metalloprotein
Reads supporting (aligned to +/- strand):  ref base A (41/48);  new base T (3/8);  total (44/56)
Fisher's exact test for biased strand distribution p-value = 3.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 101 total aligned reads displayed.

GAACCATGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGATGCCCGC  >  NC_002947/6004904‑6005201
                                                                                                                                                     |                                                                                                                                                     
gAACCATGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGa                                                                                                                                                       >  2:1540190/1‑150 (MQ=255)
 aaCCATGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGAc                                                                                                                                                      >  2:1781551/1‑150 (MQ=255)
    cATGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACcgc                                                                                                                                                   >  1:1606845/1‑150 (MQ=255)
        ccaTGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGCGGTCCCCGTAg                                                                                                                                               <  1:1043492/148‑1 (MQ=255)
         aaTGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAgg                                                                                                                                              >  2:1856544/1‑150 (MQ=255)
         aaTGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAgg                                                                                                                                              >  2:1976141/1‑150 (MQ=255)
          aTGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGgc                                                                                                                                             >  2:1869932/1‑150 (MQ=255)
            gcgcAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCCGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGgcgc                                                                                                                                           <  1:732961/150‑1 (MQ=255)
            gcgcAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGgcgt                                                                                                                                           <  2:1606845/150‑2 (MQ=255)
               cAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGc                                                                                                                                        >  2:1769000/1‑150 (MQ=255)
               cAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGc                                                                                                                                        <  1:772122/150‑1 (MQ=255)
               cAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCGCCTACGCGGTCCCCGTAGGAGCGGc                                                                                                                                        <  1:2256220/150‑1 (MQ=38)
               cAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTCGACGGCCCTTTCCCGGCTAAAGCCGCGCCTACGCGGTCCCCGTAGGAGCGGc                                                                                                                                        <  2:156648/150‑1 (MQ=25)
                  tGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCttt                                                                                                                                     >  2:1379847/1‑150 (MQ=255)
                     cGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGc                                                                                                                                  >  1:1821431/1‑150 (MQ=255)
                      ggAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGCGCACTTGCAGGTTTGACGGCCTCTTCGCGCCTAAAGCCTCGCCTACGCGGTCCCCGTAGGAGCGGCTTTAGcc                                                                                                                                 <  1:1331708/150‑1 (MQ=37)
                       gAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCg                                                                                                                                >  2:1531756/1‑150 (MQ=255)
                       gAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCg                                                                                                                                <  1:1869932/150‑1 (MQ=255)
                        aGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCgg                                                                                                                               <  1:1672328/150‑1 (MQ=255)
                        aGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCgg                                                                                                                               >  1:444204/1‑150 (MQ=255)
                            cATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGaaa                                                                                                                           >  1:1462638/1‑150 (MQ=255)
                             aTTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAg                                                                                                                          <  2:2213314/150‑1 (MQ=255)
                                 ccTGGGGTATGCCGAGCATGCGGCAGAGGGGCCCACAAAGTCGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGcc                                                                                                                      <  1:712080/150‑1 (MQ=255)
                                 ccTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGcc                                                                                                                      <  2:1821431/150‑1 (MQ=255)
                                  cTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTATAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCg                                                                                                                     >  1:1591401/1‑150 (MQ=255)
                                   tGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCgg                                                                                                                                         <  2:1330320/129‑1 (MQ=255)
                                   tGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCgg                                                                                                                                         >  1:1330320/1‑129 (MQ=255)
                                   tGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCgg                                                                                                                    >  2:458506/1‑149 (MQ=255)
                                   tGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGt                                                                                                                    >  1:2174863/1‑150 (MQ=255)
                                   tGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCAGGGAAAGGGCCGt                                                                                                                    <  1:1976141/150‑1 (MQ=255)
                                   tGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCGCCTACGCGGTCCCCGTAGGAGCGGc                                                                                                                                        >  1:310335/1‑130 (MQ=38)
                                   tGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCGCCTACGCGGTCCCCGTAGGAGCGGc                                                                                                                                        <  2:310335/130‑1 (MQ=38)
                                   gggCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGTGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAGGCCTGTGCTGCGGGGTGCCCGTAGGAGCGGCTTTAGCCGCGAAGAGGCCGt                                                                                                                    <  1:990074/149‑1 (MQ=25)
                                     gCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCg                                                                                                                  <  2:1538545/150‑1 (MQ=255)
                                       gTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAg                                                                                                                <  1:1769000/150‑1 (MQ=255)
                                          tGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTc                                                                                                             >  1:372383/1‑150 (MQ=255)
                                          tGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCGCCTACGCGGTCCCCGTAGGAGCGGCtt                                                                                                                                      <  1:217932/125‑1 (MQ=37)
                                          tGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCGCCTACGCGGTCCCCGTAGGAGCGGCtt                                                                                                                                      >  2:217932/1‑125 (MQ=37)
                                           gCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCCGTTTGTTCTGGGAAAGGGCCGTCGAGGTCg                                                                                                            <  1:1484010/150‑1 (MQ=255)
                                            ccGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGa                                                                                                           >  1:2002353/1‑150 (MQ=255)
                                             cGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTTGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGAt                                                                                                          >  1:1775762/1‑150 (MQ=255)
                                              gAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATc                                                                                                         <  1:1377743/150‑1 (MQ=255)
                                              gAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATc                                                                                                         <  2:1265929/150‑1 (MQ=255)
                                                    gCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGc                                                                                                                       >  2:1127967/1‑130 (MQ=255)
                                                    gCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGc                                                                                                                       <  1:1127967/130‑1 (MQ=255)
                                                      ggCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAgg                                                                                                >  1:2267214/1‑150 (MQ=255)
                                                       gCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGa                                                                                               >  1:828150/1‑150 (MQ=255)
                                                        cAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCGGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGaa                                                                                              <  2:265282/150‑1 (MQ=255)
                                                         agagGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAAt                                                                                             <  1:2172310/150‑1 (MQ=255)
                                                         agagGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAAt                                                                                             <  2:461087/150‑1 (MQ=255)
                                                         agagGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAAt                                                                                             <  1:445018/150‑1 (MQ=255)
                                                          gagGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATg                                                                                            >  1:1017381/1‑150 (MQ=255)
                                                          gagGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATg                                                                                            >  1:2076565/1‑150 (MQ=255)
                                                          gagGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCGCCTACGCGGTCCCCGTAGGAGCGGc                                                                                                                                        >  2:866855/1‑107 (MQ=25)
                                                          gagGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCGCCTACGCGGTCCCCGTAGGAGCGGc                                                                                                                                        <  1:866855/107‑1 (MQ=37)
                                                             gTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAgg                                                                                                               <  2:2021871/129‑1 (MQ=255)
                                                             gTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAgg                                                                                                               >  1:2021871/1‑129 (MQ=255)
                                                             gTGCCAAAGGTGCGGACGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGaca                                                                                         <  2:1626427/150‑1 (MQ=255)
                                                                   aaGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACAcactca                                                                                   <  2:2174863/150‑1 (MQ=255)
                                                                            gcgcTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTc                                                                          <  1:1531756/150‑1 (MQ=255)
                                                                            gatcTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTc                                                                          <  1:1093353/147‑1 (MQ=255)
                                                                                  gtggAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTtgt                                                                    >  1:641850/1‑150 (MQ=255)
                                                                                     gAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAg                                                                 >  2:852019/1‑150 (MQ=255)
                                                                                       aTGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCa                                                               <  2:2267214/150‑1 (MQ=255)
                                                                                        tGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGACCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCaa                                                              <  1:999653/150‑1 (MQ=255)
                                                                                         ggATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGt                                                                                                              <  1:2087040/102‑1 (MQ=255)
                                                                                         ggATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGt                                                                                                              >  2:2087040/1‑102 (MQ=255)
                                                                                          gATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTc                                                                                                                   >  1:781064/1‑96 (MQ=255)
                                                                                          gATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTc                                                                                                                   <  2:781064/96‑1 (MQ=255)
                                                                                           aTTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGt                                                           <  2:1462638/150‑1 (MQ=255)
                                                                                            ttGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGGTGCTGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTg                                                          <  1:613872/150‑1 (MQ=255)
                                                                                                  aGCAGTTTCAGGTTTGACGGCCTCTTCGCGGCTAAATTATCCCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTcc                                                    <  1:122000/150‑1 (MQ=255)
                                                                                                  aGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAAAAAAAAAAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTcc                                                    <  2:1591401/150‑1 (MQ=255)
                                                                                                   gCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCg                                                                                                                     >  2:1275029/1‑85 (MQ=255)
                                                                                                   gCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCg                                                                                                                     <  1:1275029/85‑1 (MQ=255)
                                                                                                     aGTTGCAGGTTTGACGGACTTTTCGAGGCGAGAGGCGCTCCTTCGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACCCTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAg                                                 <  2:1775762/150‑1 (MQ=255)
                                                                                                      gTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTGCTACGGGGACCGGGGAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGc                                                <  1:2233481/150‑1 (MQ=255)
                                                                                                             gTTTGACTGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGAGCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCGAGAAATGAATGACACTCTCACTTTCTTTTGCCTTTGTTTGAAGTGCATTCCCAGCTGGCATg                                         >  2:713179/1‑150 (MQ=255)
                                                                                                             gTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGac                                                                                          <  2:1635190/101‑1 (MQ=255)
                                                                                                             gTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGac                                                                                          >  1:1635190/1‑101 (MQ=255)
                                                                                                             gTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAATTGCATTCCCAGCTTTCATg                                         >  2:1254347/1‑150 (MQ=255)
                                                                                                              tttGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCTCTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGt                                        <  2:1693335/150‑1 (MQ=255)
                                                                                                                tGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTaa                                      <  2:828150/150‑1 (MQ=255)
                                                                                                                tGACGGCCTCTTCGAGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTaa                                      <  1:648150/150‑1 (MQ=255)
                                                                                                                 gACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAg                                     >  1:1913483/1‑150 (MQ=255)
                                                                                                                   cGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTc                                                                          >  2:2304686/1‑111 (MQ=255)
                                                                                                                   cGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTc                                                                          <  1:2304686/111‑1 (MQ=255)
                                                                                                                          ttCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAAc                            <  2:2076565/150‑1 (MQ=255)
                                                                                                                               ggCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCTAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAAc                       <  1:1781551/150‑1 (MQ=255)
                                                                                                                               ggCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAAc                       >  1:653389/1‑150 (MQ=255)
                                                                                                                                gCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAAc                       <  2:2002353/150‑1 (MQ=255)
                                                                                                                                gCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCAAAcc                      >  1:680391/1‑150 (MQ=255)
                                                                                                                                 cTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCa                     <  1:1856544/150‑1 (MQ=255)
                                                                                                                                     aGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGaa                 <  1:1540190/150‑1 (MQ=255)
                                                                                                                                         gCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTa             >  1:1493670/1‑150 (MQ=255)
                                                                                                                                          cTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAg            >  2:492224/1‑150 (MQ=255)
                                                                                                                                              tACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGAt        >  2:2304774/1‑150 (MQ=255)
                                                                                                                                               aCGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGATg       >  1:883819/1‑150 (MQ=255)
                                                                                                                                                 ggggACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGATGcc     <  2:1017381/150‑1 (MQ=255)
                                                                                                                                                    gACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAATTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGATGCCCGc  <  1:1254347/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAACCATGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGATGCCCGC  >  NC_002947/6004904‑6005201

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