Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,051,976 C→A 5.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,051,9760CA 5.1% 761.5 / 18.8 225L15M (CTG→ATG) ‡oprB‑IIIcarbohydrate‑selective porin
Reads supporting (aligned to +/- strand):  ref base C (154/59);  new base A (3/9);  total (157/68)
Fisher's exact test for biased strand distribution p-value = 1.34e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 252 total aligned reads displayed.

GCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTACAACGCCCACGG  >  NC_002947/4051827‑4052124
                                                                                                                                                     |                                                                                                                                                    
 ccGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGctggccct                                                                                                                                                     <  2:6371428/150‑1 (MQ=255)
 ccGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCGCCTACAAGGCCCCGGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATTCCATCCGCAAATCGCATCACCCCCACCCTGCTGcgggcccg                                                                                                                                                     >  2:422114/1‑149 (MQ=39)
      aaaCCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCa                                                                                                                                                <  2:3610779/150‑1 (MQ=255)
       aaCCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAg                                                                                                                                               <  2:6348965/150‑1 (MQ=255)
         ccGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTTGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCAGGGCCAGca                                                                                                                                             <  1:7778849/150‑1 (MQ=255)
         ccGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGca                                                                                                                                             <  1:6860414/150‑1 (MQ=255)
         ccGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGca                                                                                                                                             <  2:2284463/150‑1 (MQ=255)
            cgTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACACGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCAGGGCCAgcagca                                                                                                                                          <  1:4511395/150‑1 (MQ=255)
            cgTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGcacca                                                                                                                                          <  1:2805418/150‑1 (MQ=255)
            cgTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCAGGGCCAgcagca                                                                                                                                          <  2:268454/150‑1 (MQ=255)
              tCGCCTGCTTCGCGGGCATGCCCGCTCCTACAAGGCCCGTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGGTGCTGGCCAGGGCCAgcagcagg                                                                                                                                        <  1:7147451/150‑3 (MQ=38)
                gCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGc                                                                                                                                      <  2:8433891/150‑1 (MQ=255)
                gCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCAGGGCCAgcagcagggt                                                                                                                                      <  2:7907364/150‑5 (MQ=255)
                 ccTGCTTCGCGGGCAAGCCCGCTCCTACACCGCCCCTGGAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTTTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGcc                                                                                                                                     <  1:8843169/150‑1 (MQ=255)
                 ccTGCTTCGCGGGCAAGCCCCCTCCTACAAGGCCCCTGTAGGAGCGGTCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGAAATGCCATCCGCAATTCGCATCACCCCCCCCCTGCTGCTGGCCAGGGCCAgcagcagggtg                                                                                                                                     <  1:2799623/150‑6 (MQ=255)
                    gCTTCGCGGGCAAGCCCGCTCCTACACGGCCCCTGTTGGTGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCtt                                                                                                                                  <  1:1056442/150‑1 (MQ=255)
                      ttCGCGGGCAAGCCCGCTCCTACAGGGGCCTTGTAGGATCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATTCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGc                                                                                                                                >  2:7695229/1‑150 (MQ=38)
                      ttCGCGGGCAAGCCCGCTCCTACACGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGc                                                                                                                                <  2:6585918/150‑1 (MQ=255)
                               aaGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGcgacga                                                                                                                       >  2:6259322/1‑150 (MQ=255)
                                  cccGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTAGCCAGCACCACTGCCCTTGCCAGCGACGAtct                                                                                                                    >  1:9031352/1‑150 (MQ=255)
                                   ccGCTCCTCCAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCGGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGAtctc                                                                                                                   <  2:6934214/150‑1 (MQ=255)
                                    cGCTCCTACAGGGGCCTTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCa                                                                                                                  >  2:6132906/1‑150 (MQ=255)
                                    cGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAACGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCa                                                                                                                  <  2:7619455/150‑1 (MQ=255)
                                    cGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCa                                                                                                                  >  1:3917053/1‑150 (MQ=255)
                                     gCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCAt                                                                                                                 <  2:7817782/150‑1 (MQ=255)
                                     gCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCAt                                                                                                                 >  2:4781959/1‑150 (MQ=255)
                                     gCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCAt                                                                                                                 >  1:3132752/1‑150 (MQ=255)
                                     gCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCAt                                                                                                                 >  1:7030600/1‑150 (MQ=255)
                                       tCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGa                                                                                                               >  1:4576975/1‑150 (MQ=255)
                                         cTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGAcc                                                                                                             >  1:3942321/1‑150 (MQ=255)
                                          tACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGAccc                                                                                                            >  1:4707941/1‑150 (MQ=255)
                                          tACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGAccc                                                                                                            >  1:1863884/1‑150 (MQ=255)
                                          tACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGAccc                                                                                                            >  1:5207862/1‑150 (MQ=255)
                                          tACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGAccc                                                                                                            >  1:7505150/1‑150 (MQ=255)
                                              aGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCgcag                                                                                                        >  2:4055772/1‑150 (MQ=255)
                                                  cccTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCAcc                                                                                                    >  1:4886058/1‑150 (MQ=255)
                                                  cccTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCAcc                                                                                                    >  1:5562637/1‑150 (MQ=255)
                                                  cccTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCAcc                                                                                                    >  2:7475751/1‑150 (MQ=255)
                                                  cccTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCAcc                                                                                                    >  2:7408780/1‑150 (MQ=255)
                                                     ggcattcgCGGGCAAGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCCTCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTg                                                                                                 >  2:3419049/8‑150 (MQ=255)
                                                      gTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGa                                                                                                <  1:9098417/150‑1 (MQ=255)
                                                       tAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGAc                                                                                               >  1:2506469/1‑150 (MQ=255)
                                                           aGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGt                                                                                           <  2:5789702/150‑1 (MQ=255)
                                                            gCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTg                                                                                          >  2:8341479/1‑150 (MQ=255)
                                                                gCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTg                                                                                      >  1:3642690/1‑150 (MQ=255)
                                                                 cTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTgg                                                                                     <  1:5836120/150‑1 (MQ=255)
                                                                     cccGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGggcg                                                                                 >  1:5045650/1‑150 (MQ=255)
                                                                      ccGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGggcgg                                                                                >  2:119326/1‑150 (MQ=255)
                                                                        gcgAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGcc                                                                              >  2:3105676/1‑150 (MQ=255)
                                                                        gcgAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGAAGCACCCTGACCGGTGACTGGGGCGGcc                                                                              >  1:8668308/1‑150 (MQ=255)
                                                                               ccccTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCa                                                                       <  1:6378087/150‑1 (MQ=255)
                                                                                cccTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCAt                                                                      >  1:6576008/1‑150 (MQ=255)
                                                                                cccTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCAt                                                                      >  1:7068065/1‑150 (MQ=255)
                                                                                  cTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCa                                                                    >  1:5998532/1‑150 (MQ=255)
                                                                                  cTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCa                                                                    >  2:6642881/1‑150 (MQ=255)
                                                                                    cAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAAc                                                                  <  1:3420612/150‑1 (MQ=255)
                                                                                     aaaGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACt                                                                 >  1:2146083/1‑150 (MQ=255)
                                                                                      aaGACACCCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCACCTc                                                                <  1:530624/150‑1 (MQ=255)
                                                                                         aCAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGc                                                                                                                                <  1:8536428/83‑1 (MQ=255)
                                                                                         aCAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGc                                                                                                                                >  2:8536428/1‑83 (MQ=255)
                                                                                         aCAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCgaa                                                             >  1:2911698/1‑150 (MQ=255)
                                                                                         aCAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCgaa                                                             <  2:3784516/150‑1 (MQ=255)
                                                                                           aaCCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCgaaga                                                           >  2:4676808/1‑150 (MQ=255)
                                                                                           aaCCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCgaaga                                                           >  1:232860/1‑150 (MQ=255)
                                                                                            aCCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCgaagag                                                          >  2:328593/1‑149 (MQ=255)
                                                                                             cccGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGaagaag                                                         >  2:1551708/1‑150 (MQ=255)
                                                                                             cccGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGaagaag                                                         >  1:5266472/1‑150 (MQ=255)
                                                                                              ccGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAagaaga                                                        <  2:8680447/150‑1 (MQ=255)
                                                                                               cGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGAc                                                       >  2:7785883/1‑150 (MQ=255)
                                                                                               cGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTACTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCACGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGAc                                                       >  1:6784974/1‑150 (MQ=255)
                                                                                                       tGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGTCCTGCGCCATCAACTCGAAGAAGACGGCGTCaa                                               <  2:8809492/150‑1 (MQ=255)
                                                                                                       tGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCaa                                               >  2:1816374/1‑150 (MQ=255)
                                                                                                            tGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGCCTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAgggca                                          >  1:6204107/1‑146 (MQ=255)
                                                                                                                aTCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCgg                                      >  1:5318725/1‑150 (MQ=255)
                                                                                                                aTCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCgg                                      >  2:4949123/1‑150 (MQ=255)
                                                                                                                aTCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCgg                                      >  2:5185490/1‑150 (MQ=255)
                                                                                                                 tCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGc                                     >  2:3828846/1‑150 (MQ=255)
                                                                                                                   cGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGa                                   >  2:2831827/1‑150 (MQ=255)
                                                                                                                    gCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGAc                                  >  1:1345562/1‑150 (MQ=255)
                                                                                                                    gCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGAc                                  <  1:5876925/150‑1 (MQ=255)
                                                                                                                       aTTCGCATCACCCCCACCCTGCTGCTGGCGCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTAc                               <  2:5333153/150‑1 (MQ=255)
                                                                                                                       aTTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTAc                               >  1:481789/1‑150 (MQ=255)
                                                                                                                        ttCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACa                              >  1:2247154/1‑150 (MQ=255)
                                                                                                                             aTCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGt                         >  1:4626648/1‑150 (MQ=255)
                                                                                                                                aCCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGaa                      >  1:1601219/1‑150 (MQ=255)
                                                                                                                                 cccccACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGaaa                     <  2:4548090/150‑1 (MQ=255)
                                                                                                                                  ccccACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAAc                    >  1:3310518/1‑150 (MQ=255)
                                                                                                                                   cccACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAAcc                   <  2:4460336/150‑1 (MQ=255)
                                                                                                                                     cACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGc                 <  2:249795/150‑1 (MQ=255)
                                                                                                                                       cccTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCt               <  1:7924392/150‑1 (MQ=255)
                                                                                                                                          tgctgctgGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTACa            <  1:6132906/150‑1 (MQ=255)
                                                                                                                                          tgctgctgGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTACa            >  1:7209841/1‑150 (MQ=255)
                                                                                                                                          tgctgctgGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTACa            >  1:3755983/1‑150 (MQ=255)
                                                                                                                                             tgctgGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTACAACg         >  1:9129374/1‑150 (MQ=255)
                                                                                                                                              gctgGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTACAACGc        >  1:4963321/1‑150 (MQ=255)
                                                                                                                                               ctggccctggccAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTACAACGcc       >  1:5578768/1‑150 (MQ=255)
                                                                                                                                                tggccctggccAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTACAACGccc      >  1:4371891/1‑150 (MQ=255)
                                                                                                                                                   ccctggccAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTACAACGccgaac   >  1:8595742/1‑146 (MQ=255)
                                                                                                                                                   ccctggccAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTACAACGCCCACg   >  1:5626766/1‑150 (MQ=255)
                                                                                                                                                    cctggccAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTACAACGCCCACgg  >  2:5555999/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCCGTAAAACCGCGTCGCCTGCTTCGCGGGCAAGCCCGCTCCTACAAGGCCCCTGTAGGAGCGGGCTTGCCCGCGAATGCCCCTCAAAGACAACCCGAGATTCTGTAATGCCATCCGCAATTCGCATCACCCCCACCCTGCTGCTGGCCCTGGCCAGCACCACTGCCCTTGCCAGCGACGATCTCATGACCCGCAGCACCCTGACCGGTGACTGGGGCGGCCTGCGCCATCAACTCGAAGAAGACGGCGTCAAGATCACCGGCGACTACAGCGGTGAAACCGCCTACAACGCCCACGG  >  NC_002947/4051827‑4052124

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