Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,184,611 C→T 6.6% intergenic (‑221/‑191) mltF ← / → purL membrane‑bound lytic murein transglycosylase F/phosphoribosylformylglycinamidine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,184,6110CT 6.6% 267.5 / 5.1 90intergenic (‑221/‑191)mltF/purLmembrane‑bound lytic murein transglycosylase F/phosphoribosylformylglycinamidine synthase
Reads supporting (aligned to +/- strand):  ref base C (40/44);  new base T (3/3);  total (43/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

TGTGGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGAAAGCCATCGAATTAGGCTAGA  >  NC_002947/1184462‑1184757
                                                                                                                                                     |                                                                                                                                                  
tgtgGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTcac                                                                                                                                                    >  2:1234441/1‑150 (MQ=255)
tgtgGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTcac                                                                                                                                                    >  2:870909/1‑150 (MQ=255)
  tgGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAg                                                                                                                                                  >  2:391894/1‑150 (MQ=255)
    gTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGgc                                                                                                                                                >  2:1340535/1‑150 (MQ=255)
           gAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAAt                                                                                                                                         <  2:182529/150‑1 (MQ=255)
             ggggggATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAAtgt                                                                                                                                       >  1:1161768/1‑150 (MQ=255)
              gggggATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAAtgtg                                                                                                                                      <  1:1371452/150‑1 (MQ=255)
               ggggATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATgtgt                                                                                                                                     <  1:686376/150‑1 (MQ=255)
                gggATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAAtgt                                                                                                                                       <  1:1643266/147‑1 (MQ=255)
                gggATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAAtgt                                                                                                                                       >  2:1643266/1‑147 (MQ=255)
                gggATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTc                                                                                                                                    >  2:847623/1‑150 (MQ=255)
                     tGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGcc                                                                                                                               <  2:922552/150‑1 (MQ=255)
                        aTCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGACTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCgg                                                                                                                             <  2:199836/149‑1 (MQ=255)
                        aTCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGACTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCgg                                                                                                                             >  1:199836/1‑149 (MQ=255)
                             ggCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCt                                                                                                                       >  2:165154/1‑150 (MQ=255)
                              gCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTg                                                                                                                      <  2:46299/150‑1 (MQ=255)
                               cTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACTCCTgc                                                                                                                     <  1:391894/150‑1 (MQ=255)
                                  ccAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTgcggc                                                                                                                  >  1:462163/1‑150 (MQ=255)
                                  ccAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTgcggc                                                                                                                  >  1:462069/1‑150 (MQ=255)
                                   cAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGgcgc                                                                                                                                              <  2:1742870/121‑1 (MQ=255)
                                   cAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGgcgc                                                                                                                                              >  1:1742870/1‑121 (MQ=255)
                                    aGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGgcgc                                                                                                                <  2:224610/150‑1 (MQ=255)
                                      gTGGCCGGGCAGAAATAGGCAGCGCCTGTGAGATTGAGCGCCGCATGGTGTCGGCGACACCAAAAACCGCGCGGGTGAGACCGAGCTCCGCGCGGGCGCCGCTCGATCTTACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTc                                                                                                              <  1:1592590/150‑1 (MQ=255)
                                         gCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGAt                                                                                                           <  1:1890336/150‑1 (MQ=255)
                                         gCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGAt                                                                                                           <  2:1161768/150‑1 (MQ=255)
                                          ccGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGAtc                                                                                                          >  2:1742661/1‑150 (MQ=255)
                                              cTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGt                                                                                                      <  1:847623/150‑1 (MQ=255)
                                                  agaTATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCa                                                                                                  >  1:1881518/1‑150 (MQ=255)
                                                   gaTATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCaa                                                                                                 <  1:1269814/150‑1 (MQ=255)
                                                    atatTCAGAGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAg                                                                                                >  2:355055/1‑150 (MQ=255)
                                                     tatTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGc                                                                                               >  2:313474/1‑150 (MQ=255)
                                                      atTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTTGAGCCGGACACCTGCGGAGCTCGATCTCGGCTCAAgaa                                                                                              >  2:707434/1‑148 (MQ=255)
                                                       ttCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAg                                                                                                                                 >  1:152769/1‑114 (MQ=255)
                                                       ttCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAg                                                                                                                                 <  2:152769/114‑1 (MQ=255)
                                                       ttCAGCGCCTGTGACATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCg                                                                                             <  1:881876/150‑1 (MQ=255)
                                                        tCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCgg                                                                                            <  1:623991/150‑1 (MQ=255)
                                                           gcgcCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATGGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCTTACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGaca                                                                                         >  2:192797/1‑150 (MQ=255)
                                                              ccTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACct                                                                                      <  2:1881518/150‑1 (MQ=255)
                                                               cTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCg                                                                                             <  1:1551560/142‑1 (MQ=255)
                                                               cTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCg                                                                                             >  2:1551560/1‑142 (MQ=255)
                                                                  tgAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGa                                                                                  <  1:870909/150‑1 (MQ=255)
                                                                  tgAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGa                                                                                  >  2:37353/1‑150 (MQ=255)
                                                                     gaTCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGaaaa                                                                               >  1:1776246/1‑150 (MQ=255)
                                                                         gAGCGCCGCCCGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGa                                                                           >  2:697913/1‑150 (MQ=255)
                                                                         gAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACctct                                                                                    >  2:1397661/1‑141 (MQ=255)
                                                                         gAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACctct                                                                                    <  1:1397661/141‑1 (MQ=255)
                                                                          aGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGAc                                                                          <  1:1516212/150‑1 (MQ=255)
                                                                           gcgcCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGAACACGTGCCTCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGAcc                                                                         <  1:37353/150‑1 (MQ=255)
                                                                             gccgccAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGAcccc                                                                       >  2:641422/1‑150 (MQ=255)
                                                                                gccAGGTGTCCGCCTTACGATACTATGCAGCTCTCAGCTAGACAAGCGCGGCGGCTGCGGTGGAGCGCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGaa                                                                    <  2:908646/150‑1 (MQ=255)
                                                                                     gtgtCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCg                                                               >  2:944558/1‑150 (MQ=255)
                                                                                       gcccGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACTCCTCTGAAAACTGACCCCGAAGTCCGGt                                                             <  1:1043049/148‑1 (MQ=255)
                                                                                        tCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGtt                                                            >  2:1061274/1‑150 (MQ=255)
                                                                                         ccGCCTTAAGATGTTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACCCCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTc                                                           <  2:1769162/150‑1 (MQ=255)
                                                                                          cGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTcc                                                          >  2:114739/1‑150 (MQ=255)
                                                                                          cGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTcc                                                          >  1:60785/1‑150 (MQ=255)
                                                                                             cTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCGGGTTccgcc                                                       >  1:759898/1‑150 (MQ=255)
                                                                                                  gATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCg                                                  >  1:1728727/1‑150 (MQ=255)
                                                                                                    tatTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGc                                                <  1:313474/150‑1 (MQ=255)
                                                                                                    tatTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGc                                                <  1:114739/150‑1 (MQ=255)
                                                                                                        cAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAATTCCGGTTCCGCCGCCCGGCAGGc                                            >  2:1253110/1‑150 (MQ=255)
                                                                                                          gcgcCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGt                                          >  1:1049061/1‑150 (MQ=255)
                                                                                                           cgcCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTg                                         <  2:1728727/150‑1 (MQ=255)
                                                                                                            cccTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTgg                                        <  1:1234441/149‑1 (MQ=255)
                                                                                                              ctgtGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCATAGGCGCTGAATATCTTAAGGCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGc                                       >  2:786928/4‑149 (MQ=17)
                                                                                                              ctgtGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCATAGGCGCTGAATATCTTAAGGCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGc                                       <  1:786928/146‑1 (MQ=17)
                                                                                                              cTATGAGATCGAGCGCCGCCCGGGCGCCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCt                                      <  1:707434/150‑1 (MQ=255)
                                                                                                                gtGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCATAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCCATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCTGGCGTGCCTTa                                    >  2:437466/2‑150 (MQ=37)
                                                                                                                aTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTa                                    <  2:95685/150‑1 (MQ=255)
                                                                                                                 tGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCaa                                                                                                 <  2:1911358/88‑1 (MQ=255)
                                                                                                                 tGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCaa                                                                                                 >  1:1911358/1‑88 (MQ=255)
                                                                                                                  gagaTCGAGCGGAGCTCTGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTaga                                  <  2:49207/150‑1 (MQ=255)
                                                                                                                    gaTCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGc                                <  1:165154/150‑1 (MQ=255)
                                                                                                                    gaTCGAGCGCCGCCCGCGCGGCGCTCGATCTCATAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGc                                >  2:1907095/1‑150 (MQ=255)
                                                                                                                      tCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAg                              >  2:298566/1‑150 (MQ=255)
                                                                                                                      tCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAg                              <  1:1602044/150‑1 (MQ=255)
                                                                                                                      tCGAGCCCCGCGCGGGCGGCGCTCGATCTCATAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAg                              <  2:1885420/150‑1 (MQ=255)
                                                                                                                        gAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGt                                          >  1:1176607/1‑136 (MQ=255)
                                                                                                                        gAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGt                                          <  2:1176607/136‑1 (MQ=255)
                                                                                                                         aGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCAGAATGTGTc                                                                                                                                    <  2:1364553/45‑1 (MQ=11)
                                                                                                                         aGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCAGAATGTGTc                                                                                                                                    >  1:1364553/1‑45 (MQ=11)
                                                                                                                           cgccgcGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGcc                         <  1:641422/150‑1 (MQ=255)
                                                                                                                                cgcgGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGAAAg                    <  1:298566/150‑1 (MQ=255)
                                                                                                                                   gggcggCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAgg                                             >  2:457354/1‑122 (MQ=255)
                                                                                                                                   gggcggCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAgg                                             <  1:457354/122‑1 (MQ=255)
                                                                                                                                           cTCGATCTCATAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGAAAGCCATCGAATTa         <  1:1907095/150‑1 (MQ=255)
                                                                                                                                              gATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGAAAGCCATCGAATTAGGc      <  1:1742661/150‑1 (MQ=255)
                                                                                                                                              gATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCACGGCCGGGCAGTCGGGGCTTGGAGCAGGTGCGGAAAGCAATGGAATAAgga      >  1:408173/1‑149 (MQ=255)
                                                                                                                                               aTCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGaa                                                                    <  2:428195/87‑1 (MQ=255)
                                                                                                                                               aTCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGaa                                                                    >  1:428195/1‑87 (MQ=255)
                                                                                                                                                  tcACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGAAAGCCATCGAATTAGGCTAGa  >  2:1073697/1‑150 (MQ=255)
                                                                                                                                                  tcACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGAAAGCCATCGAATTAGGCTAGa  >  2:1770114/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TGTGGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGAAAGCCATCGAATTAGGCTAGA  >  NC_002947/1184462‑1184757

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: