Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,722,986 G→A 5.0% P217S (CCA→TCA)  PP_2383 ← AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,722,9860GA 5.0% 240.3 / 2.1 79P217S (CCA→TCA) PP_2383AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (33/42);  new base A (2/2);  total (35/44)
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.17e-01

CCCCATAGGCCTGCTTGAACTGCCGGCTCAAATGGCTCTGGTCGGCAAAGCCGAGGTCGAACGCCACTTGCGAAGCGGCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGACAGGTTGACTGCTGCCGCCAGCGCTTCCAGTGAGGGCGGGTCGGCCAGTTGGCTTTCCAGCAGTTCGCGGGCGCGGGCTACGGCCAGGGGTTCGTGG  >  NC_002947/2722842‑2723128
                                                                                                                                                |                                                                                                                                              
ccccATAGCTCTGCTTTATCTGCCGGCTCAAATGGCTCTGGTCGGCAAAGCCGAGGTCGCAGACGTGTGGTGATTCGGCCAGGCCGTGCGTCAGCCTTACGCGAGCGAGGCACCAGCCGTGTTGCTTCAGCCAGGCATGGGGTGGCAGTc                                                                                                                                           <  1:1391275/150‑1 (MQ=255)
 cccATAGGCCTGCTTGAACTGCCGGCTCAAATGGCTCTGGTCGGCAAAGCCGAGGTCGAACGCCACTTGCGAAGCGGCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTgg                                                                                                                                                <  1:147759/144‑1 (MQ=255)
 cccATAGGCCTGCTTGAACTGCCGGCTCAAATGGCTCTGGTCGGCAAAGCCGAGGTCGAACGCCACTTGCGAAGCGGCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTgg                                                                                                                                                >  2:147759/1‑144 (MQ=255)
  ccATAGGCCTGCTTGAACTGCCGGCTCAAATCACTCTGGACGACAAACCCGAAGTCGAACGCCACTTGCGAAGCGGCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCg                                                                                                                                         >  1:249336/1‑150 (MQ=255)
       ggCCTGCTTGAACTGCCGGCTCAAATGGCTCTGGTCGGCAAAGCCGAGGTCGAACGCCACTTGCGAAGCGGCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGc                                                                                                                                    <  1:1323159/150‑1 (MQ=255)
         ccTTCTTGAACTTCCGGCTCAAATGGCTCTGGTCGGCAAAGCCGAGGTCGAACGCCACTTGCGAAGCGGCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCt                                                                                                                                  <  2:1012096/150‑1 (MQ=255)
         ccTGCTTGAACTGCCGGCTCAAATGGCTCTGGTCGGCAAAGCCGAGGTCGAACGCCACTTGCGAAGCGGCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCt                                                                                                                                  >  2:1284294/1‑150 (MQ=255)
         ccTGCTTGAACTGCCGGCTCAAATGGCTCTGGTCGGCAAAGCCGAGGTCGAACGCCACTTGCGAAGCGGCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCt                                                                                                                                  <  1:907663/150‑1 (MQ=255)
           tGCTTGAACTGCCGGCTCAAATGGCTCTGGTCGGCAAAGCCGAGGTCGAACGCCACTTGCGAAGCGGCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTgg                                                                                                                                >  1:829264/1‑150 (MQ=255)
            gCTTGAACTGCCGGCTCAAATGGCTCTGGTCGGCAAAGCCGAGGTCGAACGCCACTTGCGAAGCGGCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTgga                                                                                                                               <  1:105942/150‑2 (MQ=255)
            gCTTGAACTGCCGGCTCAAATGGCTCTGGTCGGCAAAGCCGAGGTCGAACGCCACTTGCGAAGCGGCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGc                                                                                                                               <  1:1187907/150‑1 (MQ=255)
             cTTGAACTGCCGGCTCAAATGGCTCTGGTCGGCAAAGCCGAGGTCGAACGCCACTTGCGAAGCGGCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCg                                                                                                                              >  2:1245556/1‑150 (MQ=255)
              ttGAACTGCCGGCTCAAATGGCTCTGGTCGGCAAAGCCGAGGTCGAACGCCACTTGCGAAGCGGCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCGGACCGGGGGCCTGGCgg                                                                                                                             >  2:396939/1‑150 (MQ=255)
                aaaCTGACGGCTCAAATGGGTCTGGTCGGCAAAGCCGAGGTCGAACGCCACTTGCGAAGCGGCAATGCCGTGCTTCAGCATTTCGCGTGCGCGGGACAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGACTGGCGGaa                                                                                                                           <  2:713773/149‑1 (MQ=255)
                       cGGCTCAAATGGCTCTGGTCGGCAAAGCCGAGGTCGAACGCCACTTGCGAAGCGGCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCgg                                                                                                                    <  1:287398/150‑1 (MQ=255)
                           tCAAATGGCTCTGTTCGGCAAAGCCGAGGTCGAACGCCACTTGCGAAGCGGCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGCACACCCGGGCGa                                                                                                                >  2:169955/1‑150 (MQ=255)
                            cAAATGGCTCTGGTCGGCAAAGCCGAGGTCGAACGCCACTTGCGAAGCGGCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGaa                                                                                                               <  1:943558/150‑1 (MQ=255)
                              aaTGGCTCTGGTCGGCAAAGCCGAGGTCGAACGCCACTTGCGAAGCGGCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGt                                                                                                             <  2:358284/150‑1 (MQ=255)
                              aaTGGCTCTGGTCGGCAAAGCCGAGGTCGAACGCCACTTGCGAAGCGGCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGt                                                                                                             <  1:479448/150‑1 (MQ=255)
                                 ggCTCTGGTCGGCAAAGCCGAGGTCGAACGCCACTTGCGAAGCGGCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGa                                                                                                          <  1:1284294/150‑1 (MQ=255)
                                    tctGGTCGGCAAAGCCGAGGTCGAACGCCACTTGCGAAGCGGCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACg                                                                                                       >  2:778205/1‑150 (MQ=255)
                                    tctGGTCGGCAAAGCCGAGGTCGAACGCCACTTGCGAAGCGGCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACg                                                                                                       >  1:1280061/1‑150 (MQ=255)
                                     ctGGTCGGCAAAGCCGAGGTCGAACGCCACTTGCGAAGCGGCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGt                                                                                                             >  2:945915/1‑143 (MQ=255)
                                     ctGGTCGGCAAAGCCGAGGTCGAACGCCACTTGCGAAGCGGCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGt                                                                                                             <  1:945915/143‑1 (MQ=255)
                                     ctGGTCGGCAAAGCCGAGGTCGAACGCCACTTGCGAAGCGGCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACgg                                                                                                      <  1:1184606/150‑1 (MQ=255)
                                      tGGTCGGCAAAGCCGAGGTCGAACGCCACTTGCGAAGCGGCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGc                                                                                                     >  2:621369/1‑150 (MQ=255)
                                      tGGTCGGCAAAGCCGAGGTCGAACGCCACTTGCGAAGCGGCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGc                                                                                                     <  2:487873/150‑1 (MQ=255)
                                      tGGTCGGCAAAGCCGAGGTCGAACGCCACTTGCGAAGCGGCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGc                                                                                                     <  1:1275568/150‑1 (MQ=255)
                                        gTCGGCAAAGCCGAGGTCGAACGCCACTTGCGAAGCGGCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGa                                                                                                   >  1:1389044/1‑150 (MQ=255)
                                        gTCGGCAAAGCCGAGGTCGAACGCCACTTGCGAAGCGGCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGACAGTc                                                                                                                                           >  2:303839/1‑110 (MQ=255)
                                        gTCGGCAAAGCCGAGGTCGAACGCCACTTGCGAAGCGGCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGACAGTc                                                                                                                                           <  1:303839/110‑1 (MQ=255)
                                                    gAGGTCGAACGCCACTTGCGAAGCGGCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGACAGGTTGActgc                                                                                       <  1:1044983/150‑1 (MQ=255)
                                                    gAGGTCGAACGCCACTTGCGAAGCGGCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCAgg                                                                                                                                <  1:997038/109‑1 (MQ=255)
                                                    gAGGTCGAACGCCACTTGCGAAGCGGCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCAgg                                                                                                                                >  2:997038/1‑109 (MQ=255)
                                                          gAACGCCACTTGCGAAGCGGCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGACAGGTTGACTGCTgccgc                                                                                 >  1:427756/1‑150 (MQ=255)
                                                              gCCACTTGCGAAGCGGCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGACAGGtt                                                                                             <  1:845555/134‑1 (MQ=255)
                                                              gCCACTTGCGAAGCGGCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGACAGGtt                                                                                             >  2:845555/1‑134 (MQ=255)
                                                               ccACTTGCGAAGCGGCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGACAgaaca                                                                                                                                         >  1:597853/1‑85 (MQ=255)
                                                               ccACTTGCGAAGCGGCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGACAgaaca                                                                                                                                         <  2:597853/89‑5 (MQ=255)
                                                                cACTTGCGAAGCGGCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGACAGGTTGACTGCTGCCGCCAgcgc                                                                           >  1:511719/1‑150 (MQ=255)
                                                                 aCTTGCGAAGCGGCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGACAGGTTGACTGCTGCCGCCAGCGCt                                                                          <  1:778205/150‑1 (MQ=255)
                                                                    tGCGAAGCGGCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGACAGGTTGACTGCTGCCGCCAGCGCTTcc                                                                       <  2:67209/150‑1 (MQ=255)
                                                                       gAAGCGGCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGACAGGTTGACTGCTGCCGCCAGCGCTTCCAGt                                                                    >  1:401287/1‑150 (MQ=255)
                                                                         aGCGGCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGACAGGTTGACTGCTGCCGCCAGCGCTTCCAGTGa                                                                  <  1:1413169/150‑1 (MQ=255)
                                                                             gCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCGAGGCGGCGTTGTTTCCGCCAGGCATGGGGTGGCAGTTCGGGGGCCTGGTGGAACACGCGGGCGAAGTGGAAAGGCGACGGGTTGACTCCTGCCGCCAGCGCTTCCAGTGAGGGc                                                              <  2:375892/150‑1 (MQ=255)
                                                                             gCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGACAGGTTGACTGCTGCCGCCAGCGCTTCCAGTGAGGGc                                                              >  1:906969/1‑150 (MQ=255)
                                                                              ccAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGaa                                                                                                               <  1:670734/100‑1 (MQ=255)
                                                                              ccAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGaa                                                                                                               >  2:670734/1‑100 (MQ=255)
                                                                               cAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGACAGGTTGACTGCTGCCGCCAGCGCTTCCAGTGAGGGCgg                                                            <  1:261595/150‑1 (MQ=255)
                                                                                 ggCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGACAGGTTGACTGCTGCCGCCAGCGCTTCCAGTGAGGGCGGGt                                                          <  2:957974/150‑1 (MQ=255)
                                                                                  gCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGACAGGTTGACTGCTGCCGCCAGCGCTTCCAGTGAGGGCGGGTc                                                         <  2:401287/150‑1 (MQ=255)
                                                                                  gCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGACAGGTTGACTGCTGCCGCCAGCGCTTCCAGTGAGGGAGGGTc                                                         >  2:530141/1‑150 (MQ=255)
                                                                                         ttCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGACAGGTTGACTGCTGCCGCCAGCGCTTCCAGTGAGGGCGGGTCGGCCAGt                                                  >  1:316162/1‑150 (MQ=255)
                                                                                          tCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGACAGGTTGACTGCTGCCGCCAGCGCTTCCAGTGAGGGCGGGTCGGCCAGtt                                                 <  1:1294210/150‑1 (MQ=255)
                                                                                          tCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGACAGGTTGACTGCTGCCGCCAGCGCTTCCAGTGAGGGCGGGTCGGCCAGtt                                                 <  1:1245556/150‑1 (MQ=255)
                                                                                           cAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGACAGGTTGACTGCTGCCGCCAGCGCTTCCAGTGCGGGCGGGTCGGCCAGTTg                                                >  2:795191/1‑150 (MQ=255)
                                                                                           aaGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGACAGGTTGACTGCTGCCGCCAGCGCTTCCAGTGAGGGCGGGTCGGCCAGTTg                                                <  2:1160218/149‑1 (MQ=255)
                                                                                               aTTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGACAGGTTGACTGCTGCCGCCa                                                                               <  1:711240/115‑1 (MQ=255)
                                                                                               aTTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGACAGGTTGACTGCTGCCGCCa                                                                               >  2:711240/1‑115 (MQ=255)
                                                                                               aTTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGACAGGTTGACTGCTGCCGCCAGCGCTTCCAGTGAGGGCGGGTCGGCCAGTTGGCtt                                            >  1:1299790/1‑150 (MQ=255)
                                                                                               aTTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGACAGGTTGACTGCTGCCGCCAGCGCTTCCAGTGAGGGCGGGTCGGCCAGTTGGCtt                                            <  2:1389044/150‑1 (MQ=255)
                                                                                               aTTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGACAGGTTGACTGCTGCCGCCAGCGCTTCCAGTGAGGGCGGGTCGGCCAGTTGGCtt                                            >  2:1222212/1‑150 (MQ=255)
                                                                                                      gggcgcgggcCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGACAGGTTGACTGCTGCCGCCAGCGCTTCCAGTGAGGGCGGGTCGGCCAGTTGGCTTTCcagca                                     <  1:396939/150‑1 (MQ=255)
                                                                                                           cgggcCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGACAGGTTGACTGCTGCCGCCAGCGCTTCCAGTGAGGGCGGGTCGGCCAGTTGGCTTTCCAGCAGTTcg                                <  2:511719/150‑1 (MQ=255)
                                                                                                              gcctggcggcGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGACAGGTTGACTGCTGCCGCCAGCGCTTCCAGTGAGGGCGGGTCGGCCAGTTGGCTTTCCAGCAGTTcgcgg                             <  1:329082/146‑1 (MQ=255)
                                                                                                              gcCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGACAGGTTGACTGCTGCCGCCAGCGCTTCCAGTGAGGGCGGGTCGGCCAGTTGGCTTTCCAGCAGTTcgcgg                             >  1:498394/1‑150 (MQ=255)
                                                                                                                 aGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGACGACAGGTTGACTGCTGCCGCCAGCGCTTCCAGTGAGGGCGGGTCGGCCAGTTGGCTTTCCAGCAGTTcgcgggcg                          >  2:1183874/1‑150 (MQ=255)
                                                                                                                    cggcgTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGACAGGTTGACTGCTGCCGCCAGCGCTTCCAGTGAGGGCGGGTCGGCCAGTTGGCTTTCCAGCAGTTcgcgggcgcgg                       <  2:1141921/150‑1 (MQ=255)
                                                                                                                           gCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGACAGGTTGACTGCTGCCGCCAGCGCTTCCAGTGAGGGCGGGTCGGCCAGTTGGCTTTCCAGCAGTTCGCGGGCGCGGGCTACgg                >  2:950802/1‑150 (MQ=255)
                                                                                                                             ttCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGACAGGTTGACTGCTGCCGCCAGCGCTTCCAGTGAGGGCGGGTCGGCCAGTTGGCTTTCCAGCAGTTCGCGGGCGCGGGCTACGGcc              >  2:1041753/1‑150 (MQ=255)
                                                                                                                             ttCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGACAGGTTGACTGCTGCCGCCAGCGCTTCCAGTGAGGGCGGGTCGGCCAGTTGGCTTTCCAGCAGTTCGCGGGCGCGGGCTACGGcc              >  1:367949/1‑150 (MQ=255)
                                                                                                                             ttCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGACAGGTTGACTGCTGCCGCCAGCGCTTCCAGTGAGGGCGGGTCGGCCAGTTGGCTTTCCAGCAGTTCGCGGGCGCGGGCTACGGcc              >  2:1399115/1‑150 (MQ=255)
                                                                                                                              tCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGACAGGTTGACTGCTGCCGCCAGCGCTTCCAGTGAGGGCGGGTCGGCCAGTTGGCTTTCCAGCAGTTCGCGGGCGCGGGCTACGGCCa             >  2:975821/1‑150 (MQ=255)
                                                                                                                               cAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGACAGGTTGACTGCAGCCGCCAGCGCGTCCAGTGAGGGCGGGTCGGCCAGTTGGCTTTCCAGCAGTTCGCGGGCGCGGGCTACGGCCaa            <  2:1280061/150‑2 (MQ=255)
                                                                                                                                aGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGACAGGTTGACTGCTGCCGCCAGCGCTTCCAGTGAGGGCGGGTCGGCCAGTTGGCTTTCCAGCAGTTCGCGGGCGCGGGCTACGGCCAgg           <  2:367949/150‑1 (MQ=255)
                                                                                                                                  ccAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGACAGGTTGACTGCTGCCGCCAGCGCTTCCAGTGAGGGCgg                                                            >  1:691616/1‑99 (MQ=255)
                                                                                                                                  ccAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGACAGGTTGACTGCTGCCGCCAGCGCTTCCAGTGAGGGCgg                                                            <  2:691616/99‑1 (MQ=255)
                                                                                                                                     ggCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGACAGGTTGACTGCTGCCGCCAGCGCTTCCAGTGAGGGCGGGTCGGCCAGTTGGCTTTCCAGCAGTTCGCGGGCGCGGGCTACGGCCAGGGGTTc      <  1:542409/150‑1 (MQ=255)
                                                                                                                                         tGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGACAGGTTGACTGCTGCCGCCAGCGCTTCCAGTGAGGGCGGGTCGGCCAGTTGGCTTTCCAGCAGTTCGCGGGCGCGGGCTACGGCCAGGGGTTCGTgg  <  1:621369/150‑1 (MQ=255)
                                                                                                                                             gTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGACAGGTTGACTGCTGCCGCCAGCGCTTCCAGTGAGGGCGGGTCGGCCAGTTGGCTTTCCAGCAGTTCGCGGGCGCGGGCTACGGCCAGGGGTTCGTg   >  1:342226/1‑145 (MQ=255)
                                                                                                                                             gTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGACAGGTTGACTGCTGCCGCCAGCGCTTCCAGTGAGGGCGGGTCGGCCAGTTGGCTTTCCAGCAGTTCGCGGGCGCGGGCTACGGCCAGGGGTTCGTg   <  2:342226/145‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                              
CCCCATAGGCCTGCTTGAACTGCCGGCTCAAATGGCTCTGGTCGGCAAAGCCGAGGTCGAACGCCACTTGCGAAGCGGCCAGGCCGTGCTTCAGCATTTCGCGGGCGCGGGCCAGGCGGCGTTGCTTCAGCCAGGCATGGGGTGGCAGTCCGGTGGCCTGGCGGAACACCCGGGCGAAGTGGAACGGCGACAGGTTGACTGCTGCCGCCAGCGCTTCCAGTGAGGGCGGGTCGGCCAGTTGGCTTTCCAGCAGTTCGCGGGCGCGGGCTACGGCCAGGGGTTCGTGG  >  NC_002947/2722842‑2723128

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