Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,542,055 A→G 6.1% intergenic (+35/+68) nadE → / ← PP_4870 NH(3)‑dependent NAD(+) synthetase/azurin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,542,0550AG 6.1% 584.0 / 18.6 212intergenic (+35/+68)nadE/PP_4870NH(3)‑dependent NAD(+) synthetase/azurin
Reads supporting (aligned to +/- strand):  ref base A (92/92);  new base G (2/10);  total (99/113)
Fisher's exact test for biased strand distribution p-value = 3.47e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.21e-02

Only 100 of 213 total aligned reads displayed.

CGCAGATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTT  >  NC_002947/5541908‑5542202
                                                                                                                                                   |                                                                                                                                                    
 gcAGATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGa                                                                                                                                                   <  1:1687446/150‑1 (MQ=255)
 gcAGATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGa                                                                                                                                                   >  2:4497533/1‑150 (MQ=255)
  cAGATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGGTGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAg                                                                                                                                                  <  2:4408044/150‑1 (MQ=255)
   aGATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCCAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGAAGGAGCGGGTTCACCCGCGAATACGAgg                                                                                                                                                 <  1:862888/150‑1 (MQ=255)
   aGATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAgg                                                                                                                                                 <  1:1045910/150‑1 (MQ=255)
   aGATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAgg                                                                                                                                                 <  2:2383400/150‑1 (MQ=255)
   aGATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAgg                                                                                                                                                 <  2:5783640/150‑1 (MQ=255)
   aGATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAgg                                                                                                                                                 >  2:3720992/1‑150 (MQ=255)
     aTCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGt                                                                                                                                               <  1:941086/150‑1 (MQ=255)
     aTCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGt                                                                                                                                               <  1:218877/150‑1 (MQ=255)
     aTCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGt                                                                                                                                               <  2:5200507/150‑1 (MQ=255)
         aTGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTGTCCCGCGAATAAGAGGGTGTAc                                                                                                                                           >  2:5760022/1‑150 (MQ=255)
           gCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAAccc                                                                                                                                         <  1:35694/150‑1 (MQ=255)
            ccTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCa                                                                                                                                        <  1:3103096/150‑1 (MQ=255)
             cTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTGACCCGCGCATTCGAGGGGGAACCCAc                                                                                                                                      >  2:5150367/1‑150 (MQ=255)
             cTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCAc                                                                                                                                      >  2:1298772/1‑150 (MQ=255)
             cTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCAc                                                                                                                                      >  2:6857877/1‑150 (MQ=255)
             cTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTCACCCCCGAATACGAGGGTGGACCCAc                                                                                                                                      >  2:42603/1‑150 (MQ=255)
              ttCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTCACCCGCGAAGACGAGGGGGGACCCAcc                                                                                                                                     >  1:6272664/1‑150 (MQ=255)
               tCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCg                                                                                                                                    <  2:6212207/150‑1 (MQ=255)
                   gCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGc                                                                                                                                <  2:5938384/150‑1 (MQ=255)
                   gCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTAACCCCCCAAAAAGATGGTGAACccacccaccc                                                                                                                                >  1:1058858/1‑148 (MQ=255)
                      tGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCCCcgccgccat                                                                                                                             >  2:6082852/1‑147 (MQ=255)
                         gCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTcgc                                                                                                                          <  1:4341278/150‑1 (MQ=255)
                           aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACCAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                        <  2:5912657/150‑1 (MQ=255)
                           aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACACGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                        <  1:4790609/150‑1 (MQ=255)
                           aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                        <  1:6649315/150‑1 (MQ=25)
                           aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                        <  1:3152845/150‑1 (MQ=25)
                           aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                        <  2:862933/150‑1 (MQ=255)
                           aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                        <  2:235132/150‑1 (MQ=255)
                           aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                        <  1:6061204/150‑1 (MQ=255)
                           aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                        <  2:4730583/150‑1 (MQ=255)
                           aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                        <  2:2984150/150‑1 (MQ=255)
                           aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                        <  1:3396585/150‑1 (MQ=255)
                           aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                        <  1:3396633/150‑1 (MQ=255)
                             ccGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGt                                                                                                                      <  2:2491303/150‑1 (MQ=25)
                             ccGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGt                                                                                                                      <  1:5777961/150‑1 (MQ=255)
                             ccGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGt                                                                                                                      <  1:2311310/150‑1 (MQ=255)
                                 ttGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGAAc                                                                                                                  <  2:1191625/150‑1 (MQ=25)
                                 ttGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGAAc                                                                                                                  <  2:3291215/150‑1 (MQ=25)
                                    aCCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAAAACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCg                                                                                                               >  2:4151870/1‑150 (MQ=255)
                                            gCGTTCGACAACATCGTCGCCACCTCCCGCAAGACCCCGCACCCGCGCGAACTGCCGTTCGCTCCTTGATACGTGGCCTGCAGGAGCGGGTTCCCCCGCGAATACGCGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCa                                                                                                       <  2:5516160/150‑1 (MQ=37)
                                                  gACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGGGACCCGCGCAAACGGGGGGGGACCCACCGCCGCCATCGGGGGTGTACCCGCGCCCGGAggggtg                                                                                                 >  1:5543942/1‑147 (MQ=255)
                                                   cccccagcgtcgCCACCCACGGCACGACCCTGCACAGGCGCGACCTGGGGTCGGCCCCTTTATATGTGGCGTGCAGGTGCGGGTACCCCGGCGACTACGCGGTTGACCCCACCTACGCTCTCGCGGGTGACCCCGCTCCTGCAGGGTTCg                                                                                                <  2:6506315/143‑1 (MQ=255)
                                                   aCATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCAACGACGCATTCGCGGGTGAACCCGCTCCAGCAGGGTTCg                                                                                                >  2:5677271/1‑150 (MQ=255)
                                                     atcatcGTCGCCACCTACCGCACGACCCAGCACACGCGCGAACTGCCTTTCGCTCCTTTATACGTGGCCTGCAGGAGGGGGTTCCCCCGCGAATACGCGGGTGACCCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGcc                                                                                              <  1:4455859/150‑1 (MQ=37)
                                                         tcgtcgCCACCTACCGCAAGACCCAGCACACGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCa                                                                                          <  2:6220512/150‑1 (MQ=255)
                                                           gtcgCCACCTACCGCAAGACCCAGCACACGCGCGAACTGCCGTTCGCTCCTTGATACGGGGCCTGCTGGAGCGGGTTCACCCGCGAATATGAGGGTGAACCCACCGACGCATACGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGc                                                                                        <  1:3950672/150‑1 (MQ=39)
                                                               ccACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGTACCCAACGAAGCATTAGCGGGTGAAACCGCTCCTGCAGGGTTGGCCGTCAGCCatt                                                                                    >  2:2800058/1‑150 (MQ=255)
                                                               caaCCTACCCCAACACCAACCACACGCGCTAACTGCCGTTCGCTCCTTTATATGTGGCCTGCAGGTGGGGGTTCACCCGCTATTCCGCGGGTCAACCCACCTACGAATTCGCGGGTTAACCCGCTCCTGCAGGGTTCGCCGTCAGCCttt                                                                                    <  2:1865558/148‑1 (MQ=255)
                                                                cACCTACCGCAAGACCCCGCACACGCGCGAACTGCCGTTCGCTCCTTGATACGTGGCCTGCAGGAGCGGGTTCACCCGCGAATAGGCGGGTGAACCCACCGACGCATTCGGGGGTGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTg                                                                                   <  1:6285112/150‑1 (MQ=39)
                                                                 cccTACCTCAACACCCAGCACACGCGCGAACTGCCGTTCGCTCCTTTATACGTGCCCTGCTGGAGGGGGTCCACCCGCTAATACGAGGTTGACCCAACCGCCGCTTTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGa                                                                                  <  2:5152511/149‑1 (MQ=255)
                                                                   cTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAAAACGAGGGTGTAACCCCCGACGCAGTCGCGGGTGAACCCCCTCCGGCAGGGTTCGCCGTCAGccgatgttt                                                                                >  2:2597329/1‑143 (MQ=255)
                                                                     aCCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACTAGGGTGAACCCACCCACGCCTTCGCGGGTGAACCCGCTCCTGCAGGGGTCGCCGTCAGCCGTTGATTAc                                                                              >  2:3806496/1‑150 (MQ=39)
                                                                       cgaaaGACCCAGCACACGCGCTAACGGCGGTTCGCTCCTTGATATGTGGCCTGCAGGAGGGGGTTCACCCGCTAATACGCGGGTGAACCCACCGACGCTTTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACtt                                                                            <  1:5057749/147‑1 (MQ=38)
                                                                       cGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTGGATCAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACtt                                                                            <  2:5076838/150‑1 (MQ=255)
                                                                        gCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATTCGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGGTCGCCGTCAGCCTTTGATTACTTg                                                                           >  2:4445621/1‑150 (MQ=255)
                                                                               cccgCACAAGCGCTACCTGCGGTTCGCCCCTTTATATGGGGCCTGCAGGCGCGGGTCCACCCGCTAATACGTGGTTTACCCCACCGACGATTTCGGGGGTTAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGa                                                                    <  1:4932555/147‑1 (MQ=255)
                                                                               ccAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCCCCCAAAAGAGGGTGAACCCACCCACGCATTCGCGGGTGGACCCGCGCCAGCAGGGTTCGCCCTCAGCCGTTGGTTACATGACGATGa                                                                    >  2:3558693/1‑150 (MQ=255)
                                                                               ccAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTCAACCCCGAAAACCAGGGTGAACCCACCCGACCCTTCGCCGGTGAAACCGCGCCTGGCGGGTTCGCCGTCAGCCATTTGTTAATTGTCGAAGa                                                                    >  2:3131087/1‑150 (MQ=255)
                                                                                  gCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCCCGAATACGAGGGGGCACCCCCCCACGCAATACCGGGTGAACCCCCGCCCGCAGGGTTCGCCCTCCCCCGTTGATTACGTGACGATGaaca                                                                 >  2:5226328/1‑147 (MQ=255)
                                                                                           cgAACGGCCGTTCGCTCATTGATATGTGGCCTGCTGGAGCGGTTTAACCCGCTAATCCGAGGGTGAACCCACCGACGCTTTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTtca                                                        <  1:2820015/150‑1 (MQ=255)
                                                                                                 gCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGTTTCACCCGCTAATACTAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCg                                                  <  1:5675110/150‑1 (MQ=255)
                                                                                                  ccGTTCGCTCATTGATACGTGGCCTGCTGGAGCGTGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCga                                                 <  2:361726/150‑1 (MQ=255)
                                                                                                  ccGTTCGCCCATTGATACGGGGCCGGCTGGTGCGGGTACCCCCGCTATTCCGCGGGTTAACCCACCGACGCTTTCGGGGTTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCga                                                 <  1:387192/150‑1 (MQ=255)
                                                                                                    gTTCGCTCATTGATCTGTGGCCTGCTGGTGGGGTTACCCCCGGGAATTCTTGGGTGAACCAACCTCCGAATTGGCGTGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCgaga                                               <  1:1114295/150‑1 (MQ=255)
                                                                                                          tCCTTGATAAGTGGCCTGCAGGAGCGGGGTCACCCCCGCAAACGTGGGGGCACCCACCCACGCAGTCGCGGGGGGCCCCCCACCAGCACGGGTGGCCGGCCGCCGTTGATTTCGTGAAGATGGCCCTGGCGTTCATCAGCGAGAAGTggg                                         >  2:3648665/1‑149 (MQ=255)
                                                                                                            cTTGATAAGTGGCCTGCAGGAGCGGGTTCAACCCCCAATAAGAGGGTGAACCCAACCACGCAATCGCCGGTGAAACCCCAACTTCAGGGTTCGCCCTCAGCCTATGATTTATTTACCATTAACGTTCCTTTCATCAAAGAGATGTGGccc                                       >  2:2629478/1‑150 (MQ=255)
                                                                                                                   cgtgccgtGCTGGCGCGGGTACACCCGCCAATCCGCGGTTCACCCCACCGACGCCTTCGGGGGTCACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAgg                                <  1:659720/143‑1 (MQ=255)
                                                                                                                    gTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGCGGGTGAACCCACCGCCGCATTCGCGGGTGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGAGGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGt                               <  2:4993684/150‑1 (MQ=255)
                                                                                                                       gCCTGCTGGCGCGGGTTCACCCGCGAATACGCGGGTGACCCCCCCGACGAATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCa                            <  2:1246817/150‑1 (MQ=39)
                                                                                                                          tGCAGGAGCGGGTTCACCCCCGCATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAAACCCCTCCTGCAGGGTTCGCCGTCAGCCGTTTATTTCTTGTCGATGACCATGGCTTTCATCATACAGATGTGGCCCGGGCAGGTGCagaa                         >  2:6701849/1‑150 (MQ=255)
                                                                                                                           gCAGGAGCGGGTTCACCCGCGAATACGAGGGTGACCCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCagaag                        <  2:5665909/150‑1 (MQ=255)
                                                                                                                           gCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCCCTCCTGCAGGGTTCGCCGTCCGCCGTTGATTACTTGACGGTGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCagaag                        >  2:6734059/1‑150 (MQ=255)
                                                                                                                              ggAGCGGGGTCACCCGCGAATACGCGGGTGAACCCACCGACGCATACGCTGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGaagaag                     <  2:4830107/150‑1 (MQ=255)
                                                                                                                                  gtgtttCCCCGGCTAATAAGTGGGTTAACCAACCTACTCATTCGCGGGTTAACCCGCTCCTGCAGTGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGt                 <  1:4561407/146‑1 (MQ=255)
                                                                                                                                   gggTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCTAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTa                <  2:1058858/150‑1 (MQ=255)
                                                                                                                                     gTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGt              >  1:5186005/1‑150 (MQ=255)
                                                                                                                                     gTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGt              >  1:5926712/1‑150 (MQ=255)
                                                                                                                                     gTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGt              >  1:1282370/1‑150 (MQ=255)
                                                                                                                                      ttCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTc             >  2:3987337/1‑150 (MQ=255)
                                                                                                                                       tCACCCGCGAATACGAGGGTGAACCCACCGACTCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCCTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTca            <  1:4518811/150‑1 (MQ=255)
                                                                                                                                       tCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTca            >  2:2895598/1‑150 (MQ=255)
                                                                                                                                        cACCCTCGTATTCGCGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcac           <  2:667575/150‑1 (MQ=255)
                                                                                                                                         aCCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcacc          >  1:679237/1‑150 (MQ=255)
                                                                                                                                          cccGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCGGAAGAAGCTGTAGTcacca         >  2:834048/1‑150 (MQ=25)
                                                                                                                                           ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTg                                                                                                         <  1:2601349/53‑1 (MQ=255)
                                                                                                                                           ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac        >  1:1928182/1‑150 (MQ=255)
                                                                                                                                           ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac        >  2:3891478/1‑150 (MQ=255)
                                                                                                                                           ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac        >  1:3434075/1‑150 (MQ=255)
                                                                                                                                           ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac        >  2:3890322/1‑150 (MQ=255)
                                                                                                                                           ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac        >  1:4107092/1‑150 (MQ=255)
                                                                                                                                           ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac        >  2:3748000/1‑150 (MQ=255)
                                                                                                                                           ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac        >  2:2211957/1‑150 (MQ=255)
                                                                                                                                           ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac        >  2:4562518/1‑150 (MQ=255)
                                                                                                                                           ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCGGTAGTcaccac        >  1:6193623/1‑150 (MQ=255)
                                                                                                                                           ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGAAGGGTACGCCTGCACCGTTTGTTTAGTTGACGTTTACCTTGCCTTAAAACTACGAGATGTGGCCCGGGAAGTTGCAGAAGAAGATTTAGTcaccac        >  1:3931460/1‑150 (MQ=21)
                                                                                                                                            cgcgAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCCTCGAGATGTGGCCCGGGACGGTGCAGAAGAAGCTGTAGTCaccacc       >  2:6013363/1‑150 (MQ=255)
                                                                                                                                             gcgAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCg      >  2:1400614/1‑150 (MQ=255)
                                                                                                                                              cgAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCgg     >  1:813535/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGCAGATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTT  >  NC_002947/5541908‑5542202

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: