Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 668,094 A→C 8.0% intergenic (+55/+97) yfhM → / ← PP_0571 lipoprotein/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947668,0940AC 8.0% 562.8 / 5.1 215intergenic (+55/+97)yfhM/PP_0571lipoprotein/membrane protein
Reads supporting (aligned to +/- strand):  ref base A (64/116);  new base C (19/3);  total (89/126)
Fisher's exact test for biased strand distribution p-value = 9.00e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 218 total aligned reads displayed.

ACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTC  >  NC_002947/667945‑668243
                                                                                                                                                      |                                                                                                                                                     
          gccgccgcAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGcccc                                                                                                                                             <  2:1795241/150‑1 (MQ=255)
            cgccgcAGGTCGAATCGTTGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCATGGGGAGATGACGGTGAAAGCCCGCTATAGACCCGCTCTACGCGTAAAGTGTGGGTGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGa                                                                                                                                           <  1:6893730/150‑1 (MQ=39)
            cgccgcAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAACGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATGGGGCCGCTTTGCGGccctat                                                                                                                                           <  2:1281174/150‑4 (MQ=255)
             gccgcAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGAGAGGGCCGCTTTGCGGccctatc                                                                                                                                          <  2:1206498/150‑5 (MQ=255)
                  aGGTCGAACCGATGTACCGGCCGAACCTGCAGGCCGTGGGCTATGGGCATGGGGAGTTGACGGTGAAAGCCCGATTGAGACCCGCTGTACGCGGAATGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTc                                                                                                                                     <  1:7083518/150‑1 (MQ=38)
                      cGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAgc                                                                                                                                 <  2:5978779/150‑1 (MQ=255)
                       gactCGATGTACCGGCCGACCCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGCGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAgcg                                                                                                                                <  2:7104983/147‑1 (MQ=255)
                       gAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGGCGCGCGAGGGCCGCACAGCGGCCCCCACCCATCAAgcg                                                                                                                                >  2:3100447/1‑150 (MQ=255)
                        aaTCGATGTACCGGCCGAACCTGCTGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAGGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGAAAAGCGGCCCCGACAATTCAAgcgg                                                                                                                               <  2:2244850/150‑1 (MQ=255)
                        aaTCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTTGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAgcgg                                                                                                                               <  1:8702533/150‑1 (MQ=255)
                        aaTCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCGGGGCGGCGCGGTGGGGCGGCCGAGCGGCGCCGACCAttacaccgg                                                                                                                               >  1:2503843/1‑143 (MQ=255)
                         aTCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGGCGTGCGGGGCGCTGGGGCCGCCAACCGGCCCCCCCAAAGCAAacggc                                                                                                                              >  1:3864077/1‑150 (MQ=255)
                         aTCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAACGCGGCCCCGACAATTCAAgcggc                                                                                                                              <  1:6523997/150‑1 (MQ=255)
                         aTCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGGCGCGGGAGGGCCGCAAAGCGGCCCCGACAAATAAAgcggc                                                                                                                              >  2:848416/1‑150 (MQ=255)
                             aTGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGGCGCGATAGGGCCGCCAAGCGGCCCCGACAAGTCCAGCGGCGAAg                                                                                                                          >  2:8142396/1‑150 (MQ=255)
                               gTACCGGCCGAACCTGCAGGCGGGGGGCGAGGGGCCTGGGGATATGCGGGTGACCGCCCGCTAGCGACCCGCTCTCCGCGGAAAGTGGGGGCGGGGCCTTGCGTCGCGATCGGGCCGAAAAGCGGCCCCGACAATTCAAGCGGCGAAGCt                                                                                                                        <  1:6611502/150‑1 (MQ=255)
                               gTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCt                                                                                                                        <  2:7892158/150‑1 (MQ=255)
                               gTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCt                                                                                                                        <  2:7318515/150‑1 (MQ=255)
                               gTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCt                                                                                                                        <  2:654187/150‑1 (MQ=255)
                                tACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGCAGCTg                                                                                                                       >  2:3343192/1‑150 (MQ=255)
                                tACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTg                                                                                                                       <  2:5262072/150‑1 (MQ=255)
                                 aCCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGGCGCGATAGGGCCCCAAAGCCGCCCCGACAATACAATCGgcggagcgga                                                                                                                      >  2:83580/1‑150 (MQ=255)
                                  ccGGCCGAACCTGCAGGCCGTGGGCGTTGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGATACCCGCTCTACGCGGAAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGa                                                                                                                      <  2:2815495/150‑1 (MQ=255)
                                   cGGCCGAACCTGCTGGCCGTGGGCGATGGGCAAGGGGAGTTGACGGTTAAAGCCCGCTAGAGACCCGCTCTACGCGTAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGaaa                                                                                                                    <  1:5537959/150‑1 (MQ=255)
                                   cGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGaaa                                                                                                                    <  2:7815068/150‑1 (MQ=255)
                                    ggCCGAACCTGCTGGCCGTGGGCGTTGGGCATGGGGAGATGACGGTGAACGCCCGCTAGAGACCCGCTCTACGCGTAATGTGTGGGTGCGGCCTTGCGTCGCTTTAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAg                                                                                                                   <  1:5183488/150‑1 (MQ=255)
                                     gCCGACCCTGCAGGCGGGGGGGGTTGGGCATGGGGAGTTGACGGTGAAAGCCCGCTCGCGACCCGCTCTAGGGGGAATGTGGGGGTGGGGCCTTGCGGCGGTATCGGGCCGCAACGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGc                                                                                                                  <  1:2989842/150‑1 (MQ=255)
                                        gAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGTCGCGGTAGGGCCGCCAAGCGGCCCCGACAATTGAAGCGGCGTCGCGGAAAGccgg                                                                                                               >  1:1370219/1‑148 (MQ=38)
                                           ccggccggCCGTGGGCGTTGGGCATGGGGAGAGGACGGTGAACGCCCGCTAGAGACCCGCTCTACGCGGAATGTGTGGGCGCGGCCTTGCGTCGCGATAGGGCCGCAACGGGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTgggg                                                                                                            <  1:1497149/144‑1 (MQ=37)
                                           ccTGCTGGCGGTGGGCGAGGGGCATGGGGAGTTGACGGTTAAAGCCCGCTAGAGACCCGCTCTACGGGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATCGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTgggg                                                                                                            <  1:4480142/150‑1 (MQ=37)
                                                aGGCGGTGGGCGAGGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGgctgct                                                                                                       <  2:2185762/150‑1 (MQ=255)
                                                aGGCCGTGGGCGATGGGCAAGGGGAGATGACTGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCCGCGCGATAGGGCCGCAAAGCCGCACCGACAAATCAAGCGCCCAACCAGACAGACTGGGgcagca                                                                                                       >  2:7452607/1‑149 (MQ=255)
                                                aGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGGGCGGCGCGCTTGGGCCGCCAAGCGGCCCCGCCCCGTGACGCGGCGAAGCGGAAAGCCTGGGgcggct                                                                                                       >  1:3617728/1‑150 (MQ=255)
                                                  gCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGTCGCGGTATGGGCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCGGGGGCGGCggt                                                                                                     >  2:8224934/1‑150 (MQ=255)
                                                   ccGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGGCGGGCGGGGCGGTAGGGCCGCCCCGCGGCCCCGCCCAAAAAAGCCGCGCAGCGGCGAGCCGGGggcggcgggg                                                                                                    >  1:2715387/1‑146 (MQ=255)
                                                   ccGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGGCGCGGTAGGGCCGCAAAGCGGCCCCGCCAAATCAAGCGGCGCAGCGGTAAGCCGGGGGCTGCggtg                                                                                                    >  2:2363085/1‑150 (MQ=255)
                                                   ccGTGGGCGATGGCAAAGGGTAGATGACGGTTAAACCCCGCTATACACCCCCTCTACGCGGAAAGTGTGGGAGCGGCATTGCGTCGCGATAGGCCCCAAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCtgtg                                                                                                    <  1:901405/150‑1 (MQ=255)
                                                    cGTGGGCGATGGGCATGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGc                                                                                                   <  2:573335/150‑1 (MQ=255)
                                                     gTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCa                                                                                                  <  1:8618762/150‑1 (MQ=255)
                                                     gTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCCACAAATCAAGCCGCGAAGCTGAAAGCCTGGGGCTGCTGTGCa                                                                                                  >  2:8515696/1‑150 (MQ=255)
                                                     gTGGGCGATGGGCAAGGGGAGATGAAGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGAGTGGGAGCGGCCATGCCGCGCCATAGCGCCCCCCAGCCGCCCCCACAAATCAACCGGCGCACCGGCAACCCTGGGGCGGCtgtgaa                                                                                                  >  2:3805193/1‑148 (MQ=255)
                                                       gggggttGGGCATGGGGAGATGCCGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGTAATGTGTGGGTGCGGCCTTGCGTCGCTATTGGGCCCCAACGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGc                                                                                                <  1:3690017/144‑1 (MQ=255)
                                                       gggCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCAGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGc                                                                                                <  1:3557094/150‑1 (MQ=255)
                                                        ggCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGcc                                                                                               <  2:7058323/150‑1 (MQ=255)
                                                        ggCGAGGGGCACGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGcc                                                                                               <  2:4800484/150‑1 (MQ=255)
                                                          cgtggggCTAGGGGCGAGTAGGGTGAACGCCCGCTCCCGCCCCGCTCTACGCGTAACGTGTGGGTGCGGCCTTGCGTCGCGATTGGGCCGCAACGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCt                                                                                             <  1:275128/146‑1 (MQ=255)
                                                          cgcggggCACGGGGAGATGCCGGTGAACGCCCGCTCGAGCCCCGCTCTACGCGGAACGGGGGGGAGCGGCCTGGCGTCGCGATCGGGCCGCAACGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCt                                                                                             <  1:2540664/146‑1 (MQ=255)
                                                          cGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGTAAAGTGTGGGAGCGGCCTTGCGTCGCTATAGGCCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCt                                                                                             <  1:116934/150‑1 (MQ=255)
                                                          cGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGTCGCGCAAGGGCCGCAAAGCGGCCCCGACAAAAAAACCGGCGAAGCCGAAACCCAGGCGCTGCTGAGCAGacca                                                                                             >  1:7275898/1‑149 (MQ=37)
                                                                 cAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGGCAGGCGTGGCGATAGGGCCGCCAAGCGGCCCCGACCATACAAGCGGCCAAGCGGTACGCCAGGGGCGGCGGTGCCGCCCAGGCGCGa                                                                                      >  1:1027646/1‑150 (MQ=255)
                                                                 cAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATCGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcga                                                                                      >  2:5360253‑M2/1‑144 (MQ=255)
                                                                 cAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcga                                                                                      >  2:8052571‑M2/1‑144 (MQ=255)
                                                                  aaGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCCTCGCGATAGGGCCGCCAAGCGGCCCCCACAAATCAAGCGGCGAACCGGAAAGCCTGGGGCTGCAGTGCAGCCCGGtcgcgac                                                                                     <  2:4440916‑M2/150‑8 (MQ=255)
                                                                        agaTGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGACAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCACAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaagg                                                                               >  2:3507025‑M2/1‑137 (MQ=255)
                                                                         gttgaGGGTCAACGCCCGCTATACACCGGCTCTACGGGTAACGTGTGGGAGCGGCCTTGCGTCGCGATCGGGCCGCAACGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaaggc                                                                              >  1:1751967‑M2/3‑136 (MQ=255)
                                                                            gACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCCAAGCGGCCCCGACAAATAAAGCGGCGAAGCTGAAAGCCAGGGGCTGCGGTGCAGCCCAGtcgcgacgccaggccgc                                                                           <  2:2588190‑M2/150‑18 (MQ=255)
                                                                               ggTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGCGCGGCCGTGCGGCGCGCTAGGGCCGCACAGCGGCCCCGACACGGCCCGCGGCGCAGCGGCAAGCCGGGGGCGGCGGGGCCGCCCTGGCGCGGCGCCCGGCCGCGcc                                                                        >  2:5098460/1‑150 (MQ=255)
                                                                               ggTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCCAAGCGGGCCCGACAAATATAGCGGCGCAGCTGTAAGCCTGGGGCTGCTGTGCAGCCCGGtggcggcgcaaggccgctcc                                                                        <  2:2567195‑M2/150‑21 (MQ=255)
                                                                               ggTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCCAGCGGCGAAGCTGTAAGCCGGGGGCGGCTTTGCAGCCCTGtcgcgacgccaggccgcgcc                                                                        <  2:5209047‑M2/150‑21 (MQ=255)
                                                                                     aGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTAAAGCGGCGACGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCCAGGCCGCTCCCACaaa                                                                  >  2:2039655/1‑150 (MQ=255)
                                                                                       cccGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGGCGCGATAGGGCCGCAAAGCGGCCCCGACAAATCAAGCGGCGAAGCGGAAAGCCAGGGGCTGCTGGGCAGCCCGGTCGCGCCGCAAGGCCGCGCCCACCAAGa                                                                >  2:8624150/1‑150 (MQ=37)
                                                                                           cTAGAGACCCGCTCTCCGGGTAATGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAACGCGGCCCCGAAAATTCACGCGGCGAAGCTGAACGCCTGGGGCTGCTGTGCAGCCCTGTGGCGCCGCAAGGCCGCTCCCACAAAGACcgc                                                            <  2:5831604/150‑1 (MQ=25)
                                                                                              gagaCCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCAAAGCGGCCCCCACAAAAAAAGCGGCGAAGCGGGAAGCCAGGGGCTGCGGTGCAGCCCAGTCGCGACGCAAGGCCGCGCCCACAAACACcgcccc                                                         >  2:2923815/1‑147 (MQ=25)
                                                                                                caCCCGCTCTCCGCGGACTGTGGGGGCGCGGCCTTGCGTCGCGATGGGGCCGAAACGCGGCCCCGACACTTCACGCGGCGACGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGTAGa                                                       <  2:8242546/149‑1 (MQ=25)
                                                                                                 aCCCGCTCTACGCGGAAAGTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCCAAGCGGCCCCGACAAATCAAGCGGCGGAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCAGTCGCGACGCAAGGCCGCTCCCCCAAAGACCgcgccgct                                                      >  2:3398330/1‑146 (MQ=37)
                                                                                                         tACGCGGATTGTGGGGGCGCGGCGTGGGGTCGCTATCGGGCCCCAAGGCGCCCCCGACAATTCACGCGGCGACGCTGAACGCGTGGGGCTGGGGTGCCGCCGTGCCGCGCGGCACGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGAt                                              <  2:2096383/150‑1 (MQ=255)
                                                                                                           cgcgGAAAGTGTGGGAGCGGCCTTGCCTCGCGATAGGGCCGCAAAGCGGCCCCGACAAATCAACCGGCCAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCCAGGCCCCTCCCACCAAGACCGCGCCGATCAGTGGATCa                                            >  2:8727765/1‑150 (MQ=37)
                                                                                                           cgcgGAAAGTGTGGGAGCGGCCGTGCGTCGCGATAGGGCCGCAAAGCCGCCCCGACAAATCAAGCGGCCCAGCTGAAAGCCGGGGGCTGCGGGGCAGCCCAGTCGCGACGCAAGGCCGCGCCCACCAAGACCGCGCAGATCAGTGGATCa                                            >  2:3982047/1‑150 (MQ=37)
                                                                                                             cgGAACGTGGGGGTGCGGCCTTGCGTCGCGATCGGGCCGCAACGCGGCCCCGACAATTCACGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCAcc                                          <  2:7935680/150‑1 (MQ=37)
                                                                                                                    ttgtgGGTGCGCCCTTGCGAGCCTATTGGCCCCCAACGCGCCCCCCACAATACAAGCGGCGAAGCTTAAAGCCTGGGGCTGCTGTGCCGCCCTGTCGCCACGAACGCCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCtc                                    <  2:8293367/149‑1 (MQ=255)
                                                                                                                     tgtgGGAGCGGCCTTGCGGCGCGATAGGGCCGCAAAGCGGCCCCCACAAATCAAGCGGCCGAGCTGAAAGCCTGGGGCTGCTGTGCCGCCCATTAGCCTCCCCAGGCCGCGCCCCCAAAGACCCCGCCGGTCACTGGATCACCCAGCtca                                   >  2:2491074/1‑150 (MQ=21)
                                                                                                                      gtgggtggGGGCTTGGGTCGCGATCGGGCCGCAACGGGCCCCCTACAATTCACGCGGCGACGCTTAACGGGTGGGGCTGCTGTGCCGCCCTGCCGCGCCGCAAGCCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCtcat                                  <  2:4522281/150‑1 (MQ=255)
                                                                                                                       tgGGCGGGGCCTTGGGTCGCTTTCGGGCCGCAACGGGGCCCCGACAATTCACGCGGCCAAGCTGAAAGCCTGGGGCTGCTTTGCAGCCCTTTGGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCtcatc                                 <  2:1359164/150‑1 (MQ=25)
                                                                                                                         ggcgCGGCCTGGCGTGGCGTTCGGGCCGCAACGCGGCCCCGAAAATCCACGCGGCGACGCTGAACGCGTGGGGCTGCTGGGCCGCCGTGTCGCGACGCACGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCAc                               <  2:4580374/147‑1 (MQ=255)
                                                                                                                          gcgCGGCCTTGCGCCGCTATCGGGCCGAAACGCGGCCCCGACAATTAACGCGGCGACGCTGAACGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGAAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCAcc                              <  2:4346061/148‑1 (MQ=37)
                                                                                                                          gAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGcc                                                                                               <  2:5696415/85‑1 (MQ=255)
                                                                                                                           tgCGGCCTTGGGCCGCGATTGGGCCGCAACGCGGCCCCTAAAATACACGCGGCGAAGCTCAAAGCCTGGGGCTGCTGTGCCCCCCTGTCGCGACGCACGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCAccc                             <  2:8284482/149‑1 (MQ=255)
                                                                                                                           aGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCAccc                             >  2:8412343/1‑150 (MQ=255)
                                                                                                                               gCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaaggccgctcccacactttccgcgtagagcgggtctctagcgggctttcaccgtcatc                         <  1:6001027‑M2/150‑68 (MQ=255)
                                                                                                                               gCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAg                         >  1:948292/1‑150 (MQ=255)
                                                                                                                                 cTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaaggccgctcccacactttccgcgtagagcgggtctctagcgggctttcaccgtcatctc                       <  1:8346003‑M2/150‑70 (MQ=255)
                                                                                                                                  ttGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGccc                      >  1:8481763/1‑150 (MQ=255)
                                                                                                                                    gCGTCGCTATAGGCCCGCAAAGCGCCCCCCAAAATACAACCGGCTAAGCTGAAAGCCTGGGCCTGCTGTGCACCCCTGTCCCGACGAAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCa                    <  2:4925876/150‑1 (MQ=255)
                                                                                                                                      gTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGc                  >  1:5825507/1‑150 (MQ=255)
                                                                                                                                       tCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaaggccgctcccacactttccgcgtagagcgggtctctagcgggctttcaccgtcatctccccttg                 <  1:4329158‑M2/150‑76 (MQ=255)
                                                                                                                                           gATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCAccag             >  1:5807289/1‑150 (MQ=255)
                                                                                                                                            aTAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCAccagc            >  1:3770202/1‑150 (MQ=255)
                                                                                                                                              ggggCCGCTTTGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCAccagcca          >  2:7413528/2‑150 (MQ=255)
                                                                                                                                              ggggCCGCTTTGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCAccagcca          >  2:6120325/2‑150 (MQ=255)
                                                                                                                                              ggggCCGCTTTGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCAccagcca          >  1:6206114/2‑150 (MQ=255)
                                                                                                                                              ggggCCGCTTTGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCAccagcca          >  2:1695873/2‑150 (MQ=255)
                                                                                                                                               gggCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCAccagccag         >  1:7981563/1‑150 (MQ=255)
                                                                                                                                                ggCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaaggccgctcccacactttccgcgtagagcgggtctctagcgggctttcaccgtcatctccccttgcccatcgcc        <  1:6164513‑M2/150‑85 (MQ=255)
                                                                                                                                                ggCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaaggccgctcccacactttccgcgtagagcgggtctctagcgggct                                     >  1:5246947‑M2/1‑66 (MQ=255)
                                                                                                                                                ggCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaaggccgctcacacactttccgcgtagagcgggtctctagagggctttcaccgtcatctccccttgcccatcgcc        <  2:8784541‑M2/150‑85 (MQ=255)
                                                                                                                                                 gCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGtcgcgacgcaaggccgctcccacactttccgcgtagagcgggtctctagcgggctttcaccgtcatctccccttgcccatcgccc       <  2:4264194‑M2/150‑86 (MQ=255)
                                                                                                                                                 gCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAgat       >  1:2122464/1‑150 (MQ=255)
                                                                                                                                                 gCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGAGCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAgat       >  1:3224090/1‑150 (MQ=255)
                                                                                                                                                  ccGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAgatg      >  1:7526909/1‑150 (MQ=255)
                                                                                                                                                      aaaGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTc  >  1:2964616/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
ACCGCGTGCCGCCGCCGCAGGTCGAATCGATGTACCGGCCGAACCTGCAGGCCGTGGGCGATGGGCAAGGGGAGATGACGGTGAAAGCCCGCTAGAGACCCGCTCTACGCGGAAAGTGTGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAAAGCGGCCCCGACAATTCAAGCGGCGAAGCTGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTC  >  NC_002947/667945‑668243

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