Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 383,420 A→G 6.0% intergenic (+77/‑106) PP_0319 → / → PP_0320 hypothetical protein/methyl‑accepting chemotaxis transducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947383,4200AG 6.0% 791.6 / 22.3 247intergenic (+77/‑106)PP_0319/PP_0320hypothetical protein/methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base A (97/134);  new base G (6/9);  total (104/143)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 248 total aligned reads displayed.

CCTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGCGGGGAGTAGGGGCATGCGTATCTGGCGTAAGAGCATCCAGTTGCAATTGATCACCAG  >  NC_002947/383274‑383569
                                                                                                                                                  |                                                                                                                                                     
ccTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAaac                                                                                                                                                    >  2:7213776‑M2/1‑147 (MQ=255)
ccTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAaac                                                                                                                                                    >  1:3248556‑M2/1‑147 (MQ=255)
  ttGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCACATTCGCGGGTAaaccc                                                                                                                                                  <  2:1912163‑M2/150‑6 (MQ=255)
      gcgAACTCAAGCGCAGATCGCCGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATCAGCCCCGTCCACCGTCCCCCCGTAGGAGCGGGTTCCCCCGCGAATGCGTCCGTGGCTTCAACACCGCATTCGCGGGTAaacccgctc                                                                                                                                              >  2:4311460‑M2/1‑141 (MQ=255)
      gcgAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTAACCCGCGAATGCGTCAGTGGCTTCAAAACCGCATTAGCGGGGAaacccgcgc                                                                                                                                              <  1:814306‑M2/150‑10 (MQ=255)
          aCTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAaacccgctcctac                                                                                                                                          <  2:2365770‑M2/150‑14 (MQ=255)
            tCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCCTTCGCGGGTAAACCCGCTCCTACAg                                                                                                                                        >  1:3003404/1‑150 (MQ=255)
            tCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAg                                                                                                                                        >  2:579158/1‑150 (MQ=255)
                gcgcAGATCGACGCCGCGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTg                                                                                                                                    >  2:6549476/1‑150 (MQ=255)
                      aTCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTg                                                                                                                              >  2:703851/1‑150 (MQ=255)
                      aTCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCACTGGCTTCAAAACCGCAATCGCGGGTAAACCCCCTCCGAAAGGGTGCAtagg                                                                                                                              >  1:6329248/1‑147 (MQ=255)
                       tCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGa                                                                                                                             <  1:1276146/150‑1 (MQ=255)
                         gACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGAtt                                                                                                                           <  1:7117982/150‑1 (MQ=255)
                         gACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTAACCCGCGCAAGCGGGAGTGGCTTAAAAAACGCATTAGCGGGGAAAACCCCACCAACACGGTGCAACAGAtt                                                                                                                           >  1:1849446/1‑150 (MQ=255)
                             ccACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCtt                                                                                                                       <  2:429477/150‑1 (MQ=255)
                             ccACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCtt                                                                                                                       <  2:5998415/150‑1 (MQ=255)
                              cACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTc                                                                                                                      <  2:1808811/150‑1 (MQ=255)
                              cACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTc                                                                                                                      <  1:6219576/150‑1 (MQ=255)
                               aCGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGGTCACCCGCCAAATCGTCACTTGCTTCAAAACCCCAATCGCGGGGAAACCCCCACCGAAAGGGTGCATCTGATTGCttga                                                                                                                     >  1:2670986/1‑148 (MQ=255)
                                 ggTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGa                                                                                                                   <  1:887132/150‑1 (MQ=255)
                                 ggTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGa                                                                                                                   <  1:4789372/150‑1 (MQ=255)
                                  gTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCAACTGATTGCGTCAGAt                                                                                                                  >  1:6536131/1‑150 (MQ=255)
                                      aTGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACt                                                                                                              <  2:3319041/150‑1 (MQ=255)
                                       tGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTg                                                                                                             <  2:1370626/150‑1 (MQ=255)
                                       tGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTg                                                                                                             <  1:6467424/150‑1 (MQ=255)
                                        gCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTgg                                                                                                            <  2:7328388/150‑1 (MQ=255)
                                         ccATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGa                                                                                                           <  2:1495068/150‑1 (MQ=255)
                                           aTCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGAtt                                                                                                         >  1:6699923/1‑150 (MQ=255)
                                            tCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCTCCCGCGAATGCGTCAGTGGCTTCAACAGCGCATTCGCGGGTAATGCCGCTGCTACAGGGTGCATGTGATTGCTTGTGATCAGTGGAttt                                                                                                        >  2:7343917/1‑149 (MQ=255)
                                                gCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTAGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTg                                                                                                    >  2:1318840/1‑150 (MQ=255)
                                                 cGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTgg                                                                                                   <  1:3680876/150‑1 (MQ=255)
                                                 cGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTgg                                                                                                   <  1:2504928/150‑1 (MQ=255)
                                                    aCCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCATATCACTGGATTAACTGGATa                                                                                                >  2:3433200/1‑150 (MQ=255)
                                                    aCCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATa                                                                                                <  2:6322328/150‑1 (MQ=255)
                                                      cATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCGACAGGGTGCATATGATTGCTTAAGATCACTGGATTAACTGGATACg                                                                                              >  1:2913310/1‑150 (MQ=255)
                                                       aTCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCGACAGGGTGCATCTGAATGCTTCAGATCACTGGATTAACTGGATACGc                                                                                             >  2:6066007/1‑150 (MQ=255)
                                                         tgCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCaa                                                                                           <  1:6033870/149‑1 (MQ=255)
                                                          gCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAg                                                                                          <  2:4079355/150‑1 (MQ=255)
                                                           cTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAgt                                                                                         <  1:2179204/150‑1 (MQ=255)
                                                                cAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGtata                                                                                    <  2:3958899/150‑1 (MQ=255)
                                                                cAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGtata                                                                                    <  1:1667624/150‑1 (MQ=255)
                                                                   tAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCCAATGCGTCAGTGGCGTCAACACCCCATTCGCGGGGAAACCCGCTCCGACAGGGTGCATCTGATTGCTTCAGATAACTGGATTAACTGGATACGCAAGTGTATacac                                                                                 >  1:4123401/1‑150 (MQ=39)
                                                                       ccccGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTGAACCCGCGGCTACAGGGTGCATCTGATTGGGTCAGATCAAAGGTTTAACGGGAGACGGCAGTGGATAAACTTTc                                                                             >  1:5809615/1‑150 (MQ=21)
                                                                           gTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCa                                                                         >  2:6977003/1‑150 (MQ=255)
                                                                                 cGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTCAACCCGCTCCGACAGGGTGCCTCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACATTCCGCAAAAcca                                                                   >  2:6245819/1‑148 (MQ=255)
                                                                                  gTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACaat                                                                  <  2:2966558/150‑1 (MQ=255)
                                                                                   tCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACaata                                                                 <  1:6894194/150‑1 (MQ=255)
                                                                                      ccccGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACaataata                                                              <  1:995189/150‑2 (MQ=255)
                                                                                       cccGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAAtgt                                                             >  1:3591021/1‑150 (MQ=255)
                                                                                       cccGTAGGAGCGGGTTCACCCGCGAATGCGCCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAAtgt                                                             >  2:6856607/1‑150 (MQ=255)
                                                                                         cGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGc                                                           <  2:6425400/150‑1 (MQ=255)
                                                                                           tAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGCgg                                                         <  1:3341729/150‑1 (MQ=255)
                                                                                           tAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGCgg                                                         >  2:5450864/1‑150 (MQ=255)
                                                                                            aGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGCggg                                                        <  1:2646523/150‑1 (MQ=255)
                                                                                                  gggTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAaacccgctcctacggggggacggtgtacggggcttattgctcagcgatggtccgcacgatggcatccaccgtggcgtcgatctgcgcttgagttcgcgcac                                                  <  1:1019465‑M2/150‑102 (MQ=255)
                                                                                                   ggTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGCGGGGAGTAgg                                                 >  2:1557690/1‑150 (MQ=255)
                                                                                                    gTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAaacccgctcctacggggggacggtgtacggggcttattgctcagcgatggtccgcacgatggcatccaccgtggcgtcgatctgcgcttgagttcgcgcaagg                                                <  1:5057843‑M2/150‑104 (MQ=255)
                                                                                                    gTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGCGGGGAGTAggg                                                >  2:2280654/1‑150 (MQ=255)
                                                                                                     ttCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGCGGGGAGTAgggg                                               >  2:5600386/1‑150 (MQ=255)
                                                                                                     ttCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGGGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGCGGGGAGTAgggg                                               >  2:6596181/1‑150 (MQ=255)
                                                                                                       cACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTg                                                                                                                                    >  2:3695806/1‑63 (MQ=255)
                                                                                                        aCCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAAc                                                                                                      >  1:1340257/1‑92 (MQ=255)
                                                                                                        aCCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGAGACGCAAGTGTATACACTTTCCGCAAAACAATACTGTGCGGGGAGTAGGGGCAt                                            <  2:5146255/150‑1 (MQ=255)
                                                                                                         cccGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGCGGGGAGTAGGGGCATg                                           <  2:1696982/150‑1 (MQ=255)
                                                                                                          ccGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTGaacccgctcctacggggggacggtgtacggggcttattgctcagcgatggtccgcacgatggcatccaccgtggcgtcgatctgcgcttgagttcgcgcaagggtctgc                                          <  2:6479331‑M2/150‑110 (MQ=255)
                                                                                                          ccGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAaacccgctcctacggggggacggtgtacggggcttattgctcagcgatggtccgcacgatggcatccaccgtggcgtcgatctgcgcttgagt                                                          <  1:3138743‑M2/134‑94 (MQ=255)
                                                                                                          ccGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAaacccgctcctacggggggacggtgtacggggcttattgctcagcgatggtccgcacgatggcatccaccgtggcgtcgatctgcgcttgagt                                                          >  2:3138743‑M2/1‑41 (MQ=255)
                                                                                                          ccGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGCGGGGAGTAGGGGCATGc                                          >  1:2930255/1‑150 (MQ=255)
                                                                                                          ccGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGCGGGGAGTAGGGGCATGc                                          >  1:1039507/1‑150 (MQ=255)
                                                                                                           cgcgAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTGaacccgctcctacggggggacggtgtacggggcttattgctcagcgatggtccgcacgatggcatccaccgtggcgtcgatctgcgcttgagttcgcgcaagggtctgct                                         <  2:3790108‑M2/150‑111 (MQ=255)
                                                                                                           cgcgAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGCGGGGAGTAGGGGCATGCg                                         >  1:4448767/1‑150 (MQ=255)
                                                                                                           cgcgAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGCGGGGAGTAGGGGCATGCg                                         >  1:734631/1‑150 (MQ=255)
                                                                                                            gcgAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTGaacccgctcctacggggggacggtgtacggggcttattgctcagcgatggtccgcacgatggcatccaccgtggcgtcgatctgcgcttgagttcgcgcaagggtctgctg                                        <  2:5864522‑M2/150‑112 (MQ=255)
                                                                                                            gcgAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGCGGGGAGTAGGGGCATGCGt                                        >  2:3620480/1‑150 (MQ=255)
                                                                                                             cgAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGCGGGGAGTAGGGGCATGCGTa                                       >  1:2375108/1‑150 (MQ=255)
                                                                                                               aaTGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAaacccgctcctacggggggacggtgtacggggcttattgctcagcgatggtccgcacgatggcatccaccgtggcgtcgatctgcgcttgagttcgcgcaagggtctgctgatg                                     <  1:2638836‑M2/150‑115 (MQ=255)
                                                                                                                 tGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGCGGGGAGTAGGGGCATGCGTATCTg                                   >  2:2883570/1‑150 (MQ=255)
                                                                                                                 tGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGCGGGGAGTAGGGGCATGCGTATCTg                                   <  1:4723706/150‑1 (MQ=255)
                                                                                                                    gTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTg                                                                                                             >  2:2742425/1‑73 (MQ=255)
                                                                                                                    gTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGCGGGGAGTAGGGGCATGCGTATCTGGCg                                <  1:2735217/150‑1 (MQ=255)
                                                                                                                       aGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGCGGGGAGTAGGGGCATGCGTATCTGGCGTaa                             >  2:4341855/1‑150 (MQ=255)
                                                                                                                       aGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGCGGGGAGTAGGGGCATGCGTATCTGGCGTaa                             >  1:1351289/1‑150 (MQ=255)
                                                                                                                         tGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGCGGGGAGTAGGGGCATGCGTATCTGGCGTAaga                           <  1:6009876/150‑1 (MQ=255)
                                                                                                                         tGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAAGGTGCGGGGAGTAGGGGCATGCGTATCTGGCGTAaga                           <  2:5676905/150‑1 (MQ=255)
                                                                                                                         gggCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGCGGGGAGTAGGGGCATGCGTATCTGGCGTAaga                           <  2:2767080/149‑1 (MQ=255)
                                                                                                                          ggCTTCAACACCGCATTCGCGGGTGaacccgctcctacggggggacggtgtacggggcttattgctcagcgatggtccgcacgatggcatccaccgtggcgtcgatctgcgcttgagttcgcgcaagggtctgctgatgctccttctgca                          >  1:5864139‑M2/1‑25 (MQ=255)
                                                                                                                            cTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGCGGGGAGTAGGGGCATGCGTATCTGGCGTAAGAGCa                        >  1:2256857/1‑150 (MQ=255)
                                                                                                                             ttCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGCGGGGAGTAGGGGCATGCGTATCTGGCGTAAGAGCAt                       >  2:1272224/1‑150 (MQ=255)
                                                                                                                             ttCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGCGGGGAGTAGGGGCATGCGTATCTGGCGTAAGAGCAt                       >  1:7407871/1‑150 (MQ=255)
                                                                                                                             ttCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGCGGGGAGTAGGGGCATGCGTATCTGGCGTAAGAGCAt                       >  1:5503276/1‑150 (MQ=255)
                                                                                                                                aaCACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGCGGGGAGTAGGGGCATGCGTATCTGGCGTAAGAGCATCCa                    >  1:7015059/1‑150 (MQ=255)
                                                                                                                                aaCACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGCGGGGAGTAGGGGCATGCGTATCTGGCGTAAGAGCATCCa                    >  2:3878418/1‑150 (MQ=255)
                                                                                                                                 acacCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGCGGGGAGTAGGGGCATGCGTATCTGGCGTAAGAGCATCCAg                   >  1:5106276/1‑150 (MQ=255)
                                                                                                                                      ccATTCGCCGCAAAACCCGCTCCTACAGGCAGCATCTCATTCCTTCAGATCACTGGAATAACTCCATACGCAAGTGTATACACTATCCCCAAAACAATAATCTGCGGGGAGTAGGGGCATGCGTATCTGGCGAAAGAGCATCCAGTTGCa              <  2:5013681/149‑1 (MQ=255)
                                                                                                                                          tCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAAGGTGCGGGGAGTAGGGGCATGCGTATCTGGCGTAAGAGCATCCAGTTGCAATTg          >  1:3720793/1‑150 (MQ=255)
                                                                                                                                          tCGCGGGTAAACCCGCTCCTACAGGGTGCAACTGATTGATTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAATACAATATTGTGCGGGGAGTAGGGGCATGCGTATCTGGCGGAAGAGCATCCAGTTGCAATTg          >  2:3001163/1‑150 (MQ=255)
                                                                                                                                              gggTAAACCCGCTCCTACAGTGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGCGGGGAGTAGGGGCATGCGTATCTGGCGTAAGAGCATCCAGTTGCAATTGATCa      <  1:2019743/150‑1 (MQ=255)
                                                                                                                                                gTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGCGGGGAGTAGGGGCATGCGTATCTGGCGTAAGAGCATCCAGTTGCAATTGATCAcc    >  2:4902540/1‑150 (MQ=255)
                                                                                                                                                 tAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGCGGGGAGTAGGGGCATGCGTATCTGGCGTAAGAGCATCCAGTTGCAATTGATCACCa   <  1:5450864/150‑1 (MQ=255)
                                                                                                                                                 tAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGCGGGGAGTAGGGGCATGCGTATCTGGCGTAAGAGCATCCAGTTGCAATTGATCACCa   <  1:2303554/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CCTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGCGGGGAGTAGGGGCATGCGTATCTGGCGTAAGAGCATCCAGTTGCAATTGATCACCAG  >  NC_002947/383274‑383569

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