Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,053,218 C→G 5.7% intergenic (+16/+27) oprB‑III → / ← PP_3571 carbohydrate‑selective porin/acetylornithine deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,053,2180CG 5.7% 706.9 / 24.5 209intergenic (+16/+27)oprB‑III/PP_3571carbohydrate‑selective porin/acetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base C (91/105);  new base G (4/8);  total (96/113)
Fisher's exact test for biased strand distribution p-value = 5.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 209 total aligned reads displayed.

TACGGCTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCAGGGCCGCACAC  >  NC_002947/4053071‑4053364
                                                                                                                                                   |                                                                                                                                                  
 aCGGCTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCTGGGGCCGCAGTGc                                                                                                                                                 <  2:2882142/150‑1 (MQ=255)
 aCGGCTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGc                                                                                                                                                 <  1:4474584/150‑1 (MQ=255)
  cGGCTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCg                                                                                                                                                <  1:3729391/150‑1 (MQ=255)
  cGGCTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCg                                                                                                                                                <  1:1453320/150‑1 (MQ=255)
       aCCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCAGTGCGGcccc                                                                                                                                           <  1:1835944/150‑1 (MQ=255)
       aCCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGcccc                                                                                                                                           >  1:3398300/1‑150 (MQ=255)
             ccACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATccgc                                                                                                                                     <  2:5271/150‑1 (MQ=255)
             ccACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATccgc                                                                                                                                     >  2:3635036/1‑150 (MQ=255)
             ccACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATccgc                                                                                                                                     >  1:750906/1‑150 (MQ=255)
                 gCCATTGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGgc                                                                                                                                 >  2:4143921/1‑150 (MQ=255)
                 gCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGgc                                                                                                                                 >  2:3737295/1‑150 (MQ=255)
                  ccATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGgcg                                                                                                                                <  1:630100/150‑1 (MQ=255)
                    aTGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCa                                                                                                                              >  2:1983207/1‑150 (MQ=255)
                    aTGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCa                                                                                                                              >  1:1151471/1‑150 (MQ=255)
                      ggCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAAc                                                                                                                            >  1:972795/1‑150 (MQ=255)
                      ggCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAAc                                                                                                                            <  2:1585969/150‑1 (MQ=255)
                        cTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAAccc                                                                                                                          >  1:4454070/1‑150 (MQ=255)
                         tGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCCACCCt                                                                                                                         >  1:2872331/1‑150 (MQ=255)
                         tGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCt                                                                                                                         >  2:37639/1‑150 (MQ=255)
                         tGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCt                                                                                                                         >  2:3038660/1‑150 (MQ=255)
                            aCCTGCGCACGGTCGTGCATTAGATCATGGAGCCAGGGGCGTTCTCGGGCATAGTCATCGTCAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGTACCGCGGGGGCCGCACTGCGGCCCCATCCGCCGGCGCATCCCTCaa                                                                                                                      <  2:1303773/150‑1 (MQ=255)
                                gcgcCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCt                                                                                                                  <  2:2359030/150‑1 (MQ=255)
                                  gcCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGt                                                                                                                >  2:2384169/1‑150 (MQ=255)
                                    ccGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGcc                                                                                                                                             <  2:565778/119‑1 (MQ=255)
                                    ccGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGcc                                                                                                                                             >  1:565778/1‑119 (MQ=255)
                                    ccGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGttt                                                                                                              <  2:2500347/150‑1 (MQ=255)
                                    ccGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGttt                                                                                                              >  2:333592/1‑150 (MQ=255)
                                      ggACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCTGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGa                                                                                                            >  2:2240337/1‑150 (MQ=255)
                                      ggACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGa                                                                                                            >  2:2499334/1‑150 (MQ=255)
                                       gACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAg                                                                                                           <  2:1093187/150‑1 (MQ=255)
                                       gACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAg                                                                                                           <  2:2310336/150‑1 (MQ=255)
                                       gACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAg                                                                                                           <  1:35136/150‑1 (MQ=255)
                                           tGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAg                                                                                                       >  2:2846074/1‑150 (MQ=255)
                                            gCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGt                                                                                                      <  2:2748339/150‑1 (MQ=255)
                                                  catcatCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTccacca                                                                                                <  2:344449/150‑1 (MQ=255)
                                                  catcatCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGACCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTccacca                                                                                                <  2:775248/150‑1 (MQ=255)
                                                   atcatcGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAg                                                                                               >  2:618385/1‑150 (MQ=255)
                                                    tcatcGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGc                                                                                              >  2:3845092/1‑150 (MQ=255)
                                                           gCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTcagca                                                                                       <  1:4209962/150‑1 (MQ=255)
                                                                gggCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCCAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATg                                                                                  <  2:1509769/150‑1 (MQ=255)
                                                                gggCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATg                                                                                  <  2:3013374/150‑1 (MQ=255)
                                                                gggCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATg                                                                                  <  2:1044408/150‑1 (MQ=255)
                                                                gggCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATg                                                                                  >  2:4052588/1‑150 (MQ=255)
                                                                 ggCGGTCGCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATgg                                                                                 <  1:302703/150‑1 (MQ=255)
                                                                       ctcGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATcac                                                                           <  2:2741647/150‑1 (MQ=255)
                                                                       ctcGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCAACCCCCGGCGCAACCCGCAAGCCTGGTTGAGGTAGTCCACCAGCCTCAGCAGCCTGGCCTcac                                                                           >  1:1305624/1‑150 (MQ=255)
                                                                        tcGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATcaca                                                                          <  2:2157106/150‑1 (MQ=255)
                                                                        tcGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATcaca                                                                          <  2:3364700/150‑1 (MQ=255)
                                                                        tcGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATcaca                                                                          >  2:3428329/1‑150 (MQ=255)
                                                                          gggCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGc                                                                        <  2:3603872/150‑1 (MQ=255)
                                                                          gggCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGc                                                                        <  2:2123817/150‑1 (MQ=255)
                                                                              aaaGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCcgcg                                                                    >  1:3065071/1‑150 (MQ=255)
                                                                              aaaGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCcgcg                                                                    >  2:906906/1‑150 (MQ=255)
                                                                                  aCATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCCAc                                                                >  1:4383576/1‑150 (MQ=255)
                                                                                    aTCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTg                                                              >  1:1700271/1‑150 (MQ=255)
                                                                                    aTCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCAGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTg                                                              >  2:3507329/1‑150 (MQ=255)
                                                                                       gACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTc                                                           >  2:241483/1‑150 (MQ=255)
                                                                                       gACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTc                                                           <  2:1901724/150‑1 (MQ=255)
                                                                                       gACAACGCACTCGTGGTAGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGTAACTGCTc                                                           <  2:2908214/150‑1 (MQ=255)
                                                                                         cAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGa                                                         <  2:1890126/150‑1 (MQ=255)
                                                                                         cAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGa                                                         <  1:1166038/150‑1 (MQ=255)
                                                                                          aaCGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGac                                                        <  2:4009948/150‑1 (MQ=255)
                                                                                          aaCGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGGTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGac                                                        >  1:3907133/1‑150 (MQ=255)
                                                                                           aCGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGaca                                                       <  1:844596/150‑1 (MQ=255)
                                                                                             gCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGGTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACt                                                     >  2:2443487/1‑150 (MQ=255)
                                                                                                cTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGAc                                                  <  2:1709146/150‑1 (MQ=255)
                                                                                                 tCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACa                                                 >  2:669101/1‑150 (MQ=255)
                                                                                                         ggCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATc                                         >  2:4256547/1‑150 (MQ=255)
                                                                                                         ggCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATc                                         >  2:2996465/1‑150 (MQ=255)
                                                                                                          gCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAACCCTGTTTGAGGTAGTCCACCAGCCACAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATcc                                        >  1:439717/1‑150 (MQ=255)
                                                                                                            cTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCgg                                      <  2:759179/150‑1 (MQ=255)
                                                                                                              gCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCt                                    >  2:2016157/1‑150 (MQ=255)
                                                                                                              gCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCt                                    >  2:470159/1‑150 (MQ=255)
                                                                                                               cAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCtt                                   <  1:4254366/150‑1 (MQ=255)
                                                                                                                 ggTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTtgt                                 <  2:1670780/150‑1 (MQ=255)
                                                                                                                 ggTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTtgt                                 <  2:3945778/150‑1 (MQ=255)
                                                                                                                  gTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCCTCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTtgtg                                >  1:4003751/1‑150 (MQ=255)
                                                                                                                     aaGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGcc                             <  1:2568784/150‑1 (MQ=255)
                                                                                                                     aaGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGcc                             <  2:3788781/150‑1 (MQ=255)
                                                                                                                      aGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGccc                            <  1:1073178/150‑1 (MQ=255)
                                                                                                                      aGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGccc                            <  1:1018643/150‑1 (MQ=255)
                                                                                                                      aGGCGACCTTCTGAACCGCCGGGGCCGCAATGCGGCCCCAAACGCCGGCGCAACCCTCAAGACTGTTTGAGGTATTCCACAAGACACAGCAGAATGGCATCACAGCCGCGCAAAAGCAAGCCAATGACAAAATCATCCGGCTTGTGGccc                            >  1:2162422/1‑150 (MQ=255)
                                                                                                                        gCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTg                          >  2:3212224/1‑150 (MQ=255)
                                                                                                                          gACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGt                        <  2:3269268/150‑1 (MQ=255)
                                                                                                                              ttCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCAt                    >  2:1911252/1‑150 (MQ=255)
                                                                                                                              ttCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCAt                    <  1:3093684/150‑1 (MQ=255)
                                                                                                                               tCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCAtg                   >  2:315637/1‑150 (MQ=255)
                                                                                                                                 tGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCCtgct                 >  2:1374223/1‑150 (MQ=255)
                                                                                                                                  gAACCGCCGGGGCCCCACTGCGGCCCCAACCGCCCGCGCAACCCACAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCACCAAGGCAACACAGCCGCGCAACAGCTCGACACAGACAAACTCATCCGGCGTGTGGCCCAGGTcccagcag                >  1:3839156/1‑144 (MQ=255)
                                                                                                                                   aaCCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCAtgctgc               >  1:1050512/1‑150 (MQ=255)
                                                                                                                                   aaCCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCGGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCAtgctgc               >  2:2659560/1‑150 (MQ=255)
                                                                                                                                      cgccgGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCAg            >  1:4249107/1‑150 (MQ=255)
                                                                                                                                         cgGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCAGGGc         >  2:1221284/1‑150 (MQ=255)
                                                                                                                                          ggggCCGCAGTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCAGGGcc        >  2:1060131/1‑150 (MQ=255)
                                                                                                                                          ggggCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCAGGGcc        <  1:2643470/150‑1 (MQ=255)
                                                                                                                                          ggggCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCAGGGcc        <  1:2272614/150‑1 (MQ=255)
                                                                                                                                          ggggCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCAGGGcc        <  1:1983207/150‑1 (MQ=255)
                                                                                                                                            ggCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCAGGGCCGc      <  1:3003552/150‑1 (MQ=255)
                                                                                                                                             gCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCAGGGCCGca     <  1:3794129/150‑1 (MQ=255)
                                                                                                                                                gCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGACCCTGGTCCATGCTGCCAGGGCCGCacac  >  1:3006532/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TACGGCTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCAGGGCCGCACAC  >  NC_002947/4053071‑4053364

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