Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,005,053 A→T 16.1% 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,0530AT16.1% 233.8 / 32.4 95intergenic (+35/‑141)PP_5259/ydcJLysR family transcriptional regulator/metalloprotein
Reads supporting (aligned to +/- strand):  ref base A (35/43);  new base T (6/9);  total (42/53)
Fisher's exact test for biased strand distribution p-value = 7.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

GAACCATGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGATGCCCGCC  >  NC_002947/6004904‑6005202
                                                                                                                                                     |                                                                                                                                                     
gAACCATGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGa                                                                                                                                                       <  2:1386162/150‑1 (MQ=255)
    cATGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACcgc                                                                                                                                                   <  2:1080679/150‑1 (MQ=255)
     aTGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTCGACGGCCCTTTCCCGGCTAAAGCCGCGCCTACGCGGTCcccg                                                                                                                                                  >  1:219267/1‑150 (MQ=37)
      tGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGt                                                                                                                                                 >  2:1063869/1‑150 (MQ=255)
       gCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGTGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTa                                                                                                                                                <  1:2264589/150‑1 (MQ=255)
       gCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTa                                                                                                                                                >  1:2397688/1‑150 (MQ=255)
        cAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGCGGTCCCCGTAg                                                                                                                                               <  1:1996496/150‑1 (MQ=255)
         aaTGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGAGGCTAAAGCCGCGCCTACGCGGTCCCCGTAgg                                                                                                                                              >  1:604593/1‑150 (MQ=255)
          aTGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGgc                                                                                                                                             >  1:1773649/1‑150 (MQ=255)
            gcgcAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGgcgc                                                                                                                                           <  1:1698805/150‑1 (MQ=255)
            gcgcAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGAGGCTAAAGCCGCGCCTACGCGGTCCCCGTAGgagc                                                                                                                                           <  2:604593/150‑1 (MQ=255)
            gcgcAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCCCGGCTAAAGCCGCGCCTACGCGGTCCCCGTAGgagc                                                                                                                                           <  1:368832/150‑1 (MQ=38)
            gcgcAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTCGACGGCCCTTTCCCGGCTAAAGCCGCGCCTACGCGGTCCCCGTAGgagc                                                                                                                                           <  2:1995017/150‑1 (MQ=25)
              gcAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCGCCTACGCGGTCCCCGTAGGAGCgg                                                                                                                                         <  2:2256209/150‑1 (MQ=39)
              gcAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTCGACGGCCCTTTCCCGGCTAAAGCCGCGCCTACGCGGTCCCCGTAGGAGCgg                                                                                                                                         <  2:219267/150‑1 (MQ=37)
               cAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGc                                                                                                                                        <  2:2397688/150‑1 (MQ=255)
                  tGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCttt                                                                                                                                     <  2:1967620/150‑1 (MQ=255)
                   gCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTa                                                                                                                                    <  1:588426/150‑1 (MQ=255)
                            cATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGaaa                                                                                                                           >  1:623369/1‑150 (MQ=255)
                             aTTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAg                                                                                                                          >  2:2247518/1‑150 (MQ=255)
                               tACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAggg                                                                                                                        <  2:317685/150‑1 (MQ=255)
                                  cTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCGCCTACGCGGTCCCCGTAGGAGCGGCTTTAGcc                                                                                                                                 >  1:733190/1‑138 (MQ=37)
                                  cTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCGCCTACGCGGTCCCCGTAGGAGCGGCTTTAGcc                                                                                                                                 <  2:733190/138‑1 (MQ=37)
                                   tGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTa                                                                                                                                                <  2:2094642/122‑1 (MQ=255)
                                   tGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTa                                                                                                                                                >  1:2094642/1‑122 (MQ=255)
                                        tATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAgg                                                                                                               <  2:1936596/150‑1 (MQ=255)
                                         aTGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGt                                                                                                              <  2:623369/150‑1 (MQ=255)
                                         aTGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGt                                                                                                              <  1:312713/150‑1 (MQ=255)
                                              gatcATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATc                                                                                                         <  1:165927/147‑1 (MQ=255)
                                              gAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGCCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATc                                                                                                         <  2:1819186/150‑1 (MQ=255)
                                                       gCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGa                                                                                                >  2:2176199/1‑150 (MQ=255)
                                                        atctggTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAATTGCATGTTTGACGGCCTCTTCGCGGCTAAAGCCCCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGaa                                                                                               <  1:759293/146‑1 (MQ=255)
                                                           agGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGa                                                                                            <  2:1247959/150‑1 (MQ=255)
                                                             gTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACAGCGTAGGCGCGGCGTTAAGAGGGAAAGTGTCGTCGTGGGTGATGAAAAAAGGAAtagat                                                                                          >  2:759293/1‑146 (MQ=255)
                                                                ccAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCAACCTCGCGGCTAAAGCGTCTCCTACGGGGGCCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGacacac                                                                                       >  1:902401/1‑150 (MQ=255)
                                                                 cAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACAcact                                                                                      >  2:1996832/1‑150 (MQ=255)
                                                                  aaaGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACAcactc                                                                                     >  2:2070669/1‑150 (MQ=255)
                                                                    aGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTa                                                                                                                                    <  1:436247/101‑1 (MQ=255)
                                                                    aGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTa                                                                                                                                    >  2:436247/1‑101 (MQ=255)
                                                                      gTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACtt                                                                                 >  2:1196545/1‑150 (MQ=255)
                                                                          gggCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACtttttt                                                                             <  1:1914492/150‑1 (MQ=255)
                                                                           ggCGCTGGGGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGTCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACttttttt                                                                            <  2:260608/150‑1 (MQ=255)
                                                                              gcTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGc                                                                         <  1:2070669/150‑1 (MQ=255)
                                                                                tggtggAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCg                                                                                                                     >  1:1734651/1‑104 (MQ=255)
                                                                                tggtggAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCg                                                                                                                     <  2:1734651/104‑1 (MQ=255)
                                                                                   tggAATGGATGGGCGAGCAGTTGCAGGTTTGACGCCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTtgtg                                                                    <  1:1971293/150‑1 (MQ=255)
                                                                                    ggAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGa                                                                   <  1:2176199/150‑1 (MQ=255)
                                                                                    ggAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGa                                                                   <  2:902401/150‑1 (MQ=255)
                                                                                     gAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGt                                                                                                              >  1:1387139/1‑106 (MQ=255)
                                                                                     gAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGt                                                                                                              <  2:1387139/106‑1 (MQ=255)
                                                                                       aTGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACAcactc                                                                                     <  1:1396042/129‑1 (MQ=255)
                                                                                       aTGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACAcactc                                                                                     >  2:1396042/1‑129 (MQ=255)
                                                                                       aTGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCa                                                                >  1:315623/1‑150 (MQ=255)
                                                                                               gCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCAt                                                        >  1:2166630/1‑150 (MQ=255)
                                                                                                cGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTAACTTTTTTTCGCTTGTGAGCAAAGTGCAtt                                                       <  1:2247518/150‑1 (MQ=255)
                                                                                                   gCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTccc                                                    >  2:465965/1‑150 (MQ=255)
                                                                                                      gTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGc                                                 <  1:1996832/150‑1 (MQ=255)
                                                                                                      gTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGc                                                 <  1:1196545/150‑1 (MQ=255)
                                                                                                          cAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTc                                             >  1:788529/1‑150 (MQ=255)
                                                                                                           aGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCa                                            <  1:1063869/150‑1 (MQ=255)
                                                                                                             gTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCAAGCTTTCATg                                          >  1:2370247/1‑150 (MQ=255)
                                                                                                               ttGACGGCCCCTTCGCGGCTAAAGCCGCTTCTACGGAGACCGCGTAGGTGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTa                                        <  1:1338731/150‑1 (MQ=255)
                                                                                                                tGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGGCGGCGTAGGAGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTaa                                       <  2:2370247/150‑1 (MQ=255)
                                                                                                                tGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTaa                                       <  1:465965/150‑1 (MQ=255)
                                                                                                                    ggCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGt                                   >  1:1359660/1‑150 (MQ=255)
                                                                                                                     gCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGa                                                                   <  2:2367136/117‑1 (MQ=255)
                                                                                                                     gCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGa                                                                   >  1:2367136/1‑117 (MQ=255)
                                                                                                                       ctctTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCg                                >  1:1680591/1‑150 (MQ=255)
                                                                                                                         ctTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGaa                              <  2:1773649/150‑1 (MQ=255)
                                                                                                                          ttCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAAc                             <  2:1680591/150‑1 (MQ=255)
                                                                                                                            cgcgGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAAccc                           <  1:637988/150‑1 (MQ=255)
                                                                                                                             gcgGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGAt                                                                                                          >  1:732931/1‑70 (MQ=255)
                                                                                                                             gcgGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGAt                                                                                                          <  2:732931/70‑1 (MQ=255)
                                                                                                                                gCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAAcc                       >  2:1549399/1‑150 (MQ=255)
                                                                                                                                 ttcaggCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAAc                             <  1:2328010/138‑1 (MQ=255)
                                                                                                                                 ttcaggCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAAc                             >  2:2328010/6‑143 (MQ=255)
                                                                                                                                  tAAAGCCGCTCCTATGAGGACCGCGTAGGCGCGTCTTTATCCAGGAAAGGCCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAg                     <  1:1316620/150‑1 (MQ=255)
                                                                                                                                  tAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCtt                                                                       >  2:895288/1‑100 (MQ=255)
                                                                                                                                  tAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCtt                                                                       <  1:895288/100‑1 (MQ=255)
                                                                                                                                      gCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGa                                                                   <  1:1579663/100‑1 (MQ=255)
                                                                                                                                      gCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGa                                                                   >  2:1579663/1‑100 (MQ=255)
                                                                                                                                      gCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAg                 <  2:2166630/150‑1 (MQ=255)
                                                                                                                                      gCCGCGCCTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTtgt                                                                     <  2:1027280/98‑1 (MQ=25)
                                                                                                                                      gCCGCGCCTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTtgt                                                                     >  1:1027280/1‑98 (MQ=21)
                                                                                                                                      gCCGCGCCTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGCGAAGAGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTaa                                       >  2:1299511/1‑128 (MQ=21)
                                                                                                                                      gCCGCGCCTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGCGAAGAGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTaa                                       <  1:1299511/128‑1 (MQ=21)
                                                                                                                                       ccGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGc                >  2:1662258/1‑150 (MQ=255)
                                                                                                                                         gcgccTACGCGGTCCCCGTAGGAGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTa              >  2:1379432/4‑150 (MQ=39)
                                                                                                                                           tCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGc            >  2:1411925/1‑150 (MQ=255)
                                                                                                                                            ccTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCg           >  2:1724044/1‑150 (MQ=255)
                                                                                                                                             cTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGa          >  1:211295/1‑150 (MQ=255)
                                                                                                                                             cTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGa          >  1:1275352/1‑150 (MQ=255)
                                                                                                                                             cTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATGCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGa          >  1:2242326/1‑150 (MQ=255)
                                                                                                                                              tACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGa                                                                   <  2:1937243/92‑1 (MQ=255)
                                                                                                                                              tACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGa                                                                   >  1:1937243/1‑92 (MQ=255)
                                                                                                                                                    gACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGATGCCCGc   <  2:917360/150‑1 (MQ=255)
                                                                                                                                                     aCCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGATGCCCGcc  >  2:606710/1‑150 (MQ=255)
                                                                                                                                                     aCCGCGTAGGAGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGATGCCCGcc  >  2:1371888/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAACCATGCAATGCGCAGTGCCGGAGCGCATTACCTGGCGTATGCCGAGCATGCGGCAGAGGTGCCAAAGGTGCGGGCGCTGGTGGAATGGATTGGCGAGCAGTTGCAGGTTTGACGGCCTCTTCGCGGCTAAAGCCGCTCCTACGGGGACCGCGTAGGCGCGGCTTTAGCCGGGAAAGGGCCGTCGAGGTCGATCAAAAAAGGAATGACACACTCACTTTTTTTCGCTTGTGAGCAAAGTGCATTCCCAGCTTTCATGTAAGCGTCCGAACCCAACCAGGAAGCTAGCGATGCCCGCC  >  NC_002947/6004904‑6005202

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