Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,161,048 A→G 9.4% intergenic (+30/‑179) PP_3661 → / → PP_3662 membrane protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,161,0480AG 9.4% 499.4 / 28.1 171intergenic (+30/‑179)PP_3661/PP_3662membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (61/91);  new base G (8/8);  total (72/99)
Fisher's exact test for biased strand distribution p-value = 5.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 172 total aligned reads displayed.

GCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAG  >  NC_002947/4160900‑4161197
                                                                                                                                                    |                                                                                                                                                         
gCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGGAA                                                                                                                                                          >  2:4024546‑M1/1‑150 (MQ=255)
  ttCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAAC                                                                                                                                                        <  1:2337863‑M1/150‑1 (MQ=255)
    cAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCC                                                                                                                                                      <  1:4098217‑M1/150‑1 (MQ=255)
       ccAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCGGCTTCGCGGGGAAAACCccg                                                                                                                                                   >  1:4459699‑M1/1‑147 (MQ=255)
        cAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTC                                                                                                                                                  <  1:1153071‑M1/150‑1 (MQ=255)
            tgacatggcgttgaggcatctaagtctttcctgtgcgggtaaaacatccttctttgacaccaccataccgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagctgcCTGCTTCGCGGGTTcacccgcgaagca                                                                                                                                              >  2:972401‑M1/124‑144 (MQ=255)
            tgacatggcgttgaggcatctaagtctttcctgtgcgggtaaaacatccttctttgacaccaccataccgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagctgcCTGCTTCGCGGGTTcacccgcgaagca                                                                                                                                              >  1:3752913‑M1/124‑144 (MQ=255)
              cGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑G                                                                                                                                         <  1:2757678‑M1/150‑1 (MQ=255)
                 cAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCGTCGCGGGTAAAACCCCTCCGACA‑‑‑ggtg                                                                                                                                      >  2:2328033/1‑150 (MQ=255)
                     cTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑ggtgggtg                                                                                                                                  <  1:115600/150‑1 (MQ=255)
                      tGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGtgc                                                                                                                                 <  1:1279164/150‑1 (MQ=255)
                       gggCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCGGCGTCGCGGGTAAACCCCCGCCTCCA‑‑‑CGgggggggg                                                                                                                                >  2:2991758‑M1/1‑142 (MQ=255)
                            aCGGGCGCTGGCCACCGGTGCGATCCTCACCCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCt                                                                                                                           <  1:2470535/150‑1 (MQ=255)
                             cGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTg                                                                                                                          <  2:1277212/150‑1 (MQ=255)
                             cGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTg                                                                                                                          <  2:72828/150‑1 (MQ=255)
                             cGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTg                                                                                                                          <  2:72833/150‑1 (MQ=255)
                              gggCGCTGGCCACCGGTGCGATCCTCCACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGa                                                                                                                         <  2:2268477/150‑1 (MQ=255)
                              gggCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGa                                                                                                                         <  1:4358675/150‑1 (MQ=255)
                               ggCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGaa                                                                                                                        <  2:1566321/150‑1 (MQ=255)
                               ggCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGaa                                                                                                                        <  2:3540629/150‑1 (MQ=255)
                                gcgcTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAg                                                                                                                       <  2:3385944/150‑1 (MQ=255)
                                    tGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGGTGGGTCGCGGGGGAAAACCCACCGACA‑‑‑GGGGGGGGGGGGCCGGAAGACgg                                                                                                                   >  2:200119/1‑150 (MQ=255)
                                       ttcctgtgcgggtaaaacatccttctttgacaccaccataccgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagctgcCTGCTTCGCGGGTGAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCaa                                                                                                                >  1:3413944‑M1/97‑150 (MQ=255)
                                       ccACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCaa                                                                                                                <  2:1803499/150‑1 (MQ=255)
                                         aCCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATGACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAgg                                                                                                              <  1:653233/150‑1 (MQ=255)
                                          ccGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGGGGCGTGGCCGGGGAACCCCCCCCAACA‑‑‑GCTGGGGGGGGCCCCGGAGACCGccaagc                                                                                                             >  1:1393519/1‑145 (MQ=255)
                                          ccGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGc                                                                                                             <  1:1178921/150‑1 (MQ=255)
                                            ggTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCc                                                                                                           >  2:1926708‑M2/1‑149 (MQ=255)
                                              tGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCcct                                                                                                         >  2:1087762‑M2/1‑147 (MQ=255)
                                               gCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGggccttgccgtcttcagggcagcacccac                                                                                                           <  1:2174084‑M1/150‑30 (MQ=255)
                                                cGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCcctgt                                                                                                       >  2:2072510‑M2/1‑145 (MQ=255)
                                                 ggtaaaacatccttctttgacaccaccataccgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagctgcCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTa                                                                                                      >  1:3167637‑M1/87‑150 (MQ=255)
                                                 gATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCcctgta                                                                                                      >  2:4369021‑M2/1‑144 (MQ=255)
                                                 gATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGCCTGCTTCGCGGGTGAACCCGCTCCTACA‑‑‑GGggccttgccgtcttcagggcagcacccacctgta                                                                                                      <  1:1911392‑M1/150‑35 (MQ=255)
                                                  gtaaaacatccttctttgacaccaccataccgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagctgcCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAg                                                                                                     >  2:1345283‑M1/86‑150 (MQ=255)
                                                  aTCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGggccttgccgtcttcagggcagcacccacctgtag                                                                                                     <  1:2337236‑M1/150‑36 (MQ=255)
                                                    ccTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGggccttgccgtcttcagggcagcacccacctgtagga                                                                                                   <  1:2345091‑M1/150‑38 (MQ=255)
                                                     cTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGCCTGCTTCGCGGGTGAACCCCCACCAACA‑‑‑GGggccgtgcccgcataagggcagcacccaccaggaggat                                                                                                  >  1:2606957‑M1/1‑112 (MQ=255)
                                                     aaacatccttctttgacaccaccataccgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagctgcCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAg                                                                                                  >  1:164796‑M1/83‑150 (MQ=255)
                                                      tCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTGAACCCGCTCCTACA‑‑‑GGggccttgccgtcttcagggcagcacccacctgtaggagc                                                                                                 <  1:4514243‑M1/150‑40 (MQ=255)
                                                      tCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCcctgtaggagc                                                                                                 >  2:3761945‑M2/1‑139 (MQ=255)
                                                      tCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCcctgtaggagc                                                                                                 >  2:1875537‑M2/1‑139 (MQ=255)
                                                        catccttctttgacaccaccataccgagaaatgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagctgcCTGCTTCGCGGGTGAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCgg                                                                                               >  1:1549013‑M1/80‑150 (MQ=255)
                                                        aaCCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCcctgtaggagcgg                                                                                               >  1:3216348‑M2/1‑137 (MQ=255)
                                                        aaCCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGCCTGCTTCGCGGGTGAACCCGCTCCTACA‑‑‑GGggccttgccgtcttcagggcagcacccacctgtaggagcgg                                                                                               <  1:3374378‑M1/150‑42 (MQ=255)
                                                        aaCCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGCCTGCTTCGCGGGTGAACCCGCTCCTACA‑‑‑GGggccttgccgtcttcagggcagcacccacctgtaggagcgg                                                                                               <  2:3663496‑M1/150‑42 (MQ=255)
                                                         ccctGTGGCTGGTTGGCGGGGGGCTGGCGATGCCCCTGGGGGTGCAGACGCTGTTGGGCTGACTGGCGACCGAGTCGGCTGCGTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCcctgtaggagcggg                                                                                              >  2:3056754‑M2/2‑136 (MQ=255)
                                                         aCCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCcctgtaggagcggg                                                                                              >  2:4253360‑M2/1‑136 (MQ=255)
                                                         aCCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCcctgtaggagcggg                                                                                              >  1:1471898‑M2/1‑136 (MQ=255)
                                                            tgtgGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTc                                                                                           <  1:2173318/150‑1 (MQ=255)
                                                              tgGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCAc                                                                                         <  1:2955000/150‑1 (MQ=255)
                                                                gCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTGAACCCGCTCCTACA‑‑‑GGggccttgccgtcttcagggcagcacccacctgtagga                                                                                                   >  2:3038951‑M1/1‑101 (MQ=255)
                                                                gCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCAccc                                                                                       <  2:3714637/150‑1 (MQ=255)
                                                                 cTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGggccttgccgtcttcagggcagcacccacctgtaggagcgg                                                                                               <  2:1348280‑M1/141‑42 (MQ=255)
                                                                  tGGTTGGCGGGGGGCTGGCGATTACCCGGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCcgc                                                                                     <  1:226271/150‑1 (MQ=255)
                                                                  tGGTTGGCGGGGGACTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCcgc                                                                                     <  2:2673565/150‑1 (MQ=255)
                                                                    acaccaccataccgagaactgtggacggcaacgaggggggctgtgaaaaaccttgagaccgagctgcTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGaa                                                                                  >  1:2309194‑M1/68‑150 (MQ=255)
                                                                     ttGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGaa                                                                                  <  1:780247/150‑1 (MQ=255)
                                                                     ttGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGaa                                                                                  <  2:755858/150‑1 (MQ=255)
                                                                        gCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAgca                                                                               <  1:1787647/150‑1 (MQ=255)
                                                                        caccctcccgagaaatgtggacggccacgaggggggctgtgaaaaaccttgagaccgagctgcCTGCTTCGCGGGTAAACCCGCTCCTCCA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAgca                                                                               >  2:3929533‑M1/64‑150 (MQ=255)
                                                                         cGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAgcag                                                                              <  2:408349/150‑1 (MQ=255)
                                                                          ggggggCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACC‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAgcagg                                                                             <  2:1121979/150‑1 (MQ=255)
                                                                                cTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTc                                                                       <  2:4129869/150‑1 (MQ=255)
                                                                                cTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTc                                                                       <  2:3003243/150‑1 (MQ=255)
                                                                                     gATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGtct                                                                  <  1:1011399/150‑1 (MQ=255)
                                                                                     gATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGtct                                                                  <  1:1011391/150‑1 (MQ=255)
                                                                                       ttACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCa                                                                <  1:3360226/150‑1 (MQ=255)
                                                                                         tggacggcaacgaggggggctgtgaaaaacctggagaccgagtcggCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAg                                                             >  2:264122‑M1/47‑150 (MQ=255)
                                                                                               ggggTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGACGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGttttt                                                       <  2:2104027/150‑1 (MQ=255)
                                                                                               ggggTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGttttt                                                       <  2:3712799/150‑1 (MQ=255)
                                                                                               gcaacgaggggggctgttaaaaacctggagaccgagtcggCTGCTTCGCGGGTAAACCCGCTCCTACC‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGttttt                                                       >  1:3522428‑M1/41‑150 (MQ=255)
                                                                                                   ggCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTcaca                                                   <  1:4030033/149‑1 (MQ=255)
                                                                                                      aGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACC‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGcc                                                <  2:1979396/150‑1 (MQ=255)
                                                                                                      aGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGcc                                                <  1:4081214/150‑1 (MQ=255)
                                                                                                       gAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCGGCGTACCCGGTAAAAACCCACCAAAA‑‑‑CGGGGGGGCTGCCCCGAAGAACGCCCGGCCCCCGGAGGGGCGGGGTCACCCCCCCAGCACGCCCCACGGGCTCAACGGtttttaaaaccc                                                >  2:3467921/1‑143 (MQ=255)
                                                                                                           cTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTTAAGACGGCAAGGCCCCTGTAGGAGGGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCt                                           <  1:3528579/150‑1 (MQ=255)
                                                                                                                      aCTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTcc                                <  2:3130445/150‑1 (MQ=255)
                                                                                                                       cTGGAGACCGTGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACC‑‑‑GGTGGGTGCTGCCCTGAAGACGGCACGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCca                               <  1:1633028/150‑1 (MQ=255)
                                                                                                                       cTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCca                               <  2:2898512/150‑1 (MQ=255)
                                                                                                                               cGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACT‑‑‑GGTGGGTGCGGCCCTGAAGACGGCACGGCCCCTGTAGGAGCGGGTTCACCCGCGACGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTtctc                       <  1:200119/150‑1 (MQ=255)
                                                                                                                               cGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTtctc                       <  1:285445/150‑1 (MQ=255)
                                                                                                                                gAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTCGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCg                      <  1:859100/150‑1 (MQ=255)
                                                                                                                                gAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCg                      <  1:1477546/150‑1 (MQ=255)
                                                                                                                                 aGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCATGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCgg                     <  1:1779282/150‑1 (MQ=255)
                                                                                                                                      gctgctTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGccc                                               >  1:1866133/1‑119 (MQ=255)
                                                                                                                                      gctgctTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGccc                                               <  2:1866133/119‑1 (MQ=255)
                                                                                                                                      gctgctTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTAtgg                <  2:1145692/150‑1 (MQ=255)
                                                                                                                                      gctgctTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTAtgg                >  2:666095/1‑150 (MQ=255)
                                                                                                                                        TGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTAtggtg              >  2:2600365‑M1/1‑150 (MQ=255)
                                                                                                                                          CTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATggtggt            >  1:4486220‑M1/1‑150 (MQ=255)
                                                                                                                                          CTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATggtggt            >  2:2445678‑M1/1‑150 (MQ=255)
                                                                                                                                          CTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATggtggt            >  2:2445684‑M1/1‑150 (MQ=255)
                                                                                                                                            TCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGgtgt          <  1:764104‑M1/150‑1 (MQ=255)
                                                                                                                                            TCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGgtgt          >  2:693981‑M1/1‑150 (MQ=255)
                                                                                                                                             CGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTc         <  2:2469254‑M1/150‑1 (MQ=255)
                                                                                                                                              gcgGGTAAACCCGCTCCTACA‑‑‑GGggccttgccgtcttcagggcagcacccacctgtaggagcgggtttacccgcgaagcagccgactcggtctccagtcagcccaacagcttctgcacccccagggtaatcgccagccccccgcccacca         >  2:2807635‑M1/1‑23 (MQ=255)
                                                                                                                                                GGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTGCACATTTCTCGGTATGGTGGTGTCaaa      <  2:2930690‑M1/150‑1 (MQ=255)
                                                                                                                                                   TAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAaagaa   >  2:2025070‑M1/1‑150 (MQ=255)
                                                                                                                                                    aaaCCCGCTCCTCCA‑‑‑GGTGGGTGCGGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAaagaag  <  1:2991758/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                         
GCTTCAGCCAGATTCGCCAGGCTGGGCCACGGGCGCTGGCCACCGGTGCGATCCTCAACCTGTGGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACA‑‑‑GGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAG  >  NC_002947/4160900‑4161197

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: