Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,050,054 C→A 9.1% L15M (CTG→ATG) ‡ oprB‑III → carbohydrate‑selective porin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,050,0540CA 9.1% 243.2 / 9.8 77L15M (CTG→ATG) ‡oprB‑IIIcarbohydrate‑selective porin
Reads supporting (aligned to +/- strand):  ref base C (41/29);  new base A (3/4);  total (44/33)
Fisher's exact test for biased strand distribution p-value = 4.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

GCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTACAACGCCCACGGT  >  NC_002947/4049905‑4050203
                                                                                                                                                     |                                                                                                                                                     
gCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGctggccc                                                                                                                                                       <  2:2100914/150‑1 (MQ=255)
     aaaaCCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGctggccctggcc                                                                                                                                                  <  2:158954/150‑1 (MQ=255)
              tCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCTGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCAGGGCCAGcagcagg                                                                                                                                         <  2:1967436/150‑6 (MQ=39)
                              cAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCa                                                                                                                                                 <  2:1538249/126‑1 (MQ=255)
                              cAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCa                                                                                                                                                 >  1:1538249/1‑126 (MQ=255)
                               aaGCCCGCTCCTACAGGGGCCTGGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCGATCCGCAATCCGCTTCACCCCCACCCTTCTGCTGGCCCTGGCCAGCACCACTGCCCTTCCCAGcgacga                                                                                                                        >  2:931991/1‑150 (MQ=38)
                               aaGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGcgacga                                                                                                                        >  2:1690860/1‑150 (MQ=255)
                                  cccGCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtct                                                                                                                     >  1:664198/1‑150 (MQ=255)
                                  cccGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCAGGGCCAGcagcagg                                                                                                                                         >  1:711978/1‑125 (MQ=255)
                                  cccGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCAGGGCCAGcagcagg                                                                                                                                         <  2:711978/130‑6 (MQ=255)
                                   ccGCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtctc                                                                                                                    >  2:2150513/1‑150 (MQ=255)
                                    cGCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCa                                                                                                                   >  2:1408924/1‑150 (MQ=255)
                                    cGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCa                                                                                                                   <  2:2290804/150‑1 (MQ=255)
                                     gCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCa                                                                                                                   <  2:785167/149‑1 (MQ=255)
                                     gCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCa                                                                                                                   >  1:785167/1‑149 (MQ=255)
                                     gCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCAt                                                                                                                  >  2:246565/1‑150 (MQ=255)
                                     gCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCAt                                                                                                                  >  1:2225999/1‑150 (MQ=255)
                                     gCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCAt                                                                                                                  >  2:174972/1‑150 (MQ=255)
                                       tCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGcga                                                                                                                           <  2:1845029/139‑1 (MQ=255)
                                       tCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGcga                                                                                                                           >  1:1845029/1‑139 (MQ=255)
                                                   ccTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGctggccctgg                                                                                                                                                    <  2:1958121/102‑1 (MQ=255)
                                                   ccTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGctggccctgg                                                                                                                                                    >  1:1958121/1‑102 (MQ=255)
                                                    cTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCTCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCt                                                                                                   >  2:197079/1‑150 (MQ=255)
                                                       tAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTg                                                                                                                                        >  2:739228/1‑110 (MQ=255)
                                                       tAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTg                                                                                                                                        <  1:739228/110‑1 (MQ=255)
                                                       tAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCCCAGATAGGCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGAc                                                                                                >  2:1289212/1‑150 (MQ=255)
                                                           aGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGt                                                                                            <  1:246565/150‑1 (MQ=255)
                                                             cGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGa                                                                                          <  2:1432408/150‑1 (MQ=255)
                                                                    gCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGggc                                                                                   <  1:174972/150‑1 (MQ=255)
                                                                     cccGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGggcg                                                                                  <  2:2225999/150‑1 (MQ=255)
                                                                     cccGCAAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCAGGGCCAgcagca                                                                                                                                           <  2:445049/93‑1 (MQ=255)
                                                                     cccGCAAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCAGGGCCAGcagca                                                                                                                                           >  1:445049/1‑90 (MQ=255)
                                                                      ccGCGAATGCCCCTCAAATACAACCCGAGATTCTGTAATGCCATCCGCAATCCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCt                                                                                                                                    <  1:654761/99‑1 (MQ=255)
                                                                      ccGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGggcgg                                                                                 >  1:1816687/1‑150 (MQ=255)
                                                                      ccGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGggcgg                                                                                 >  1:77591/1‑150 (MQ=255)
                                                                      ccGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATCCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCt                                                                                                                                    >  2:654761/1‑99 (MQ=255)
                                                                              gCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGCGACTGGGGCGGCCTGCGcc                                                                         >  1:617789/1‑150 (MQ=255)
                                                                                 ccTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGAc                                                                                                >  1:1265524/1‑124 (MQ=255)
                                                                                 ccTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGAc                                                                                                <  2:1265524/124‑1 (MQ=255)
                                                                                 ccTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGACTGCGCCATc                                                                      >  1:1412051/1‑150 (MQ=255)
                                                                                      aaGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGTGACGATCTCa                                                                                                                   >  1:489992/1‑100 (MQ=255)
                                                                                      aaGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGTGACGATCTCa                                                                                                                   <  2:489992/100‑1 (MQ=255)
                                                                                         aCAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGggcggc                                                                                <  2:355076/132‑1 (MQ=255)
                                                                                         aCAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGggcggc                                                                                >  1:355076/1‑132 (MQ=255)
                                                                                           aaCCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGa                                                                                                                >  1:2146977/1‑98 (MQ=255)
                                                                                           aaCCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGa                                                                                                                <  2:2146977/98‑1 (MQ=255)
                                                                                             ccggagaTTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGaagaag                                                          >  1:1895079/4‑150 (MQ=255)
                                                                                              ccGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGTGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAagaaga                                                         <  1:1319833/150‑1 (MQ=255)
                                                                                              ccGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAagaaga                                                         >  2:1607584/1‑150 (MQ=255)
                                                                                              ccGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAagaaga                                                         >  2:420086/1‑150 (MQ=255)
                                                                                              ccGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAagaaga                                                         >  2:878506/1‑150 (MQ=255)
                                                                                                       tGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCaa                                                <  2:664198/150‑1 (MQ=255)
                                                                                                       tGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCaa                                                <  1:1690860/150‑1 (MQ=255)
                                                                                                         tAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCa                                                                                                                   >  2:734087/1‑81 (MQ=255)
                                                                                                         tAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCa                                                                                                                   <  1:734087/81‑1 (MQ=255)
                                                                                                         tAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCgg                                                                                             >  2:881188/1‑103 (MQ=255)
                                                                                                         tAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCgg                                                                                             <  1:881188/103‑1 (MQ=255)
                                                                                                         tAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGa                                              >  1:1470826/1‑150 (MQ=255)
                                                                                                           aTGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGggcg                                                                                  <  2:1294056/112‑1 (MQ=255)
                                                                                                           aTGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGggcg                                                                                  >  1:1294056/1‑112 (MQ=255)
                                                                                                           aTGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATc                                            >  2:145080/1‑150 (MQ=255)
                                                                                                              ccATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCCGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCAcc                                         >  2:542408/1‑150 (MQ=255)
                                                                                                                 tCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGc                                      >  1:2043991/1‑150 (MQ=255)
                                                                                                                   cGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAAATCGAAGAAGACGGCGTCCAGATCACCGGCGa                                    >  2:305564/1‑150 (MQ=255)
                                                                                                                        ttCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGCGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACa                               <  2:617789/150‑1 (MQ=255)
                                                                                                                             aTCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGt                          <  1:878506/150‑1 (MQ=255)
                                                                                                                             aTCACCCCCACCCTCCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGt                          <  1:1120503/150‑1 (MQ=255)
                                                                                                                                 cccccACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGaaa                      >  1:2266365/1‑150 (MQ=255)
                                                                                                                                  ccccACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGggc                                                                                   <  2:329035/88‑1 (MQ=255)
                                                                                                                                  ccccACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGggc                                                                                   >  1:329035/1‑88 (MQ=255)
                                                                                                                                  ccccACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAAc                     >  2:2019497/1‑150 (MQ=255)
                                                                                                                                     cACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGc                  >  1:2324789/1‑150 (MQ=255)
                                                                                                                                     cACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGc                  <  1:1408924/150‑1 (MQ=255)
                                                                                                                                     cACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGc                  <  2:2324789/150‑1 (MQ=255)
                                                                                                                                       cccTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTCCGGCGGTGAAACCGCCt                <  1:1289212/150‑1 (MQ=255)
                                                                                                                                        ccTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTa               >  2:1612836/1‑150 (MQ=255)
                                                                                                                                          tgctgctgGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAg                              >  2:444207/1‑133 (MQ=255)
                                                                                                                                          tgctgctgGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAg                              <  1:444207/133‑1 (MQ=255)
                                                                                                                                           ggtgctgGCCAGGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGa                                    <  2:82067/124‑1 (MQ=255)
                                                                                                                                           ggtgctgGCCAGGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGa                                    >  1:82067/3‑126 (MQ=255)
                                                                                                                                             tgctgGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTACAACg          <  1:1607584/150‑1 (MQ=255)
                                                                                                                                                   ccctggccAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTACAACGCCCACg    >  2:989201/1‑150 (MQ=255)
                                                                                                                                                    cctggccAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTTCAACGCCCACgg   >  2:1309068/1‑150 (MQ=255)
                                                                                                                                                     ctggccAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTACAACGCCCACGGt  <  1:989201/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTACAACGCCCACGGT  >  NC_002947/4049905‑4050203

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: