Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,184,594 G→C 16.7% intergenic (‑204/‑208) 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,5940GC16.7% 215.3 / 33.3 83intergenic (‑204/‑208)mltF/purLmembrane‑bound lytic murein transglycosylase F/phosphoribosylformylglycinamidine synthase
Reads supporting (aligned to +/- strand):  ref base G (35/34);  new base C (7/7);  total (42/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

TCTGGTAGCGGGTTGTGTGGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGAAAGCCA  >  NC_002947/1184447‑1184742
                                                                                                                                                   |                                                                                                                                                     
tCTGGTAGCGGGTTGTGTGGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCgcgggc                                                                                                                                                     <  1:1683959/150‑1 (MQ=255)
    gTAGCGGGTTGTGTGGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGgcggcg                                                                                                                                                 >  2:1275892/1‑150 (MQ=255)
      aGCGGGTTGTGTGGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCt                                                                                                                                               <  1:2136510/150‑1 (MQ=255)
         gggTTGTGTGGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                                            <  2:2231146/150‑1 (MQ=255)
          ggTTGTGTGGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGAt                                                                                                                                           <  1:1940180/150‑1 (MQ=255)
           tttGTGTGGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGTCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGAtc                                                                                                                                          <  1:1008499/149‑1 (MQ=37)
           gTTGTGTGGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCTGGCGGCGCTCGaac                                                                                                                                          >  1:1696611/1‑148 (MQ=255)
           gTTGTGTGGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGAtc                                                                                                                                          <  1:2280414/150‑1 (MQ=255)
              gtgtgGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTca                                                                                                                                       <  2:1076085/150‑1 (MQ=255)
                        gggAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGa                                                                                                                             >  1:1337650/1‑150 (MQ=255)
                            ggggggATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAAtgt                                                                                                                         >  2:2061427/1‑150 (MQ=255)
                                 gATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGa                                                                                                                    <  2:1621883/150‑1 (MQ=255)
                                 gATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGa                                                                                                                    <  2:1621877/150‑1 (MQ=255)
                                       aTCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGa                                                                                                              >  1:1578981/1‑150 (MQ=255)
                                       aTCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCAGGTGTCCGGCTCGACACATTCTGCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTcata                                                                                                                                     <  2:1956008/126‑3 (MQ=11)
                                       aTCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCAGGTGTCCGGCTCGACACATTCGGCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTcata                                                                                                                                     >  1:1956008/1‑124 (MQ=11)
                                          aGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGgcgc                                                                                                                                                <  2:975286/113‑1 (MQ=255)
                                          aGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGgcgc                                                                                                                                                >  1:975286/1‑113 (MQ=255)
                                            ggCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCt                                                                                                         <  1:2061427/150‑1 (MQ=255)
                                                   aGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGgcgc                                                                                                 >  1:1183584/1‑150 (MQ=255)
                                                   aGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGgcgc                                                                                                 >  1:119740/1‑150 (MQ=255)
                                                     gTGGCCGGCTTGAGATATTCAGCGCCTGTGCGATCGTGTCGCGCAAGCTGTCCGCGTTAAGATATTCAGCGCGTATCGGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTc                                                                                               <  1:2341238/150‑1 (MQ=37)
                                                          cGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGcc                                                                                                                 <  2:1310212/128‑1 (MQ=255)
                                                          cGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGcc                                                                                                                 >  1:1310212/1‑128 (MQ=255)
                                                           ggCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGAtctc                                                                                         >  1:324485/1‑150 (MQ=255)
                                                            gCTTGAGATATTCAGCGCCTGTGAGATCGAGAGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                                                              >  1:1116300/1‑97 (MQ=255)
                                                            gCTTGAGATATTCAGCGCCTGTGAGATCGAGAGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                                                              <  2:1116300/97‑1 (MQ=255)
                                                              ttGAGATATTCGGCTGCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATAGTCAGCGCCTATGAGTTCGAGCGCCGCTCGGGCGGCGCTCTATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCTCTCGATCTCGtc                                                                                      <  2:119740/150‑1 (MQ=38)
                                                              ttGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGtc                                                                                      >  2:456040/1‑150 (MQ=255)
                                                                     atTCAGCGCCTGTGAGATCGAGAGCCGCCAGGTGTCCGCCTTAAGATATGCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGcc                                                                               >  2:1835716/1‑150 (MQ=255)
                                                                      ttCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGATATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                              >  2:1532788/1‑87 (MQ=21)
                                                                      ttCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCg                                                                              <  2:2273142/150‑1 (MQ=255)
                                                                      ttCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                              <  1:1532788/87‑1 (MQ=255)
                                                                       tCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCTGCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGa                                                                                                                                            <  1:880386/88‑1 (MQ=17)
                                                                       tCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCTGCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGa                                                                                                                                            >  2:880386/1‑88 (MQ=14)
                                                                        cAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACCCCTGCGGCGCTCGATCTCGTCTCAAGCCGGa                                                                            <  1:821138/150‑1 (MQ=255)
                                                                        cAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGa                                                                            <  2:1578981/150‑1 (MQ=255)
                                                                               tgtgAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACctct                                                                     <  2:1337650/150‑1 (MQ=255)
                                                                                    gaTCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGaaaa                                                                >  2:2154952/1‑150 (MQ=255)
                                                                                     aTCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAAc                                                               >  2:2029848/1‑150 (MQ=255)
                                                                                        gAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGa                                                            <  1:2154952/150‑1 (MQ=255)
                                                                                         aGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGAc                                                           <  1:1275892/150‑1 (MQ=255)
                                                                                          gcgcCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGAcc                                                          <  1:1972810/150‑1 (MQ=255)
                                                                                              cgccAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGa                                                      >  1:1364968/1‑150 (MQ=255)
                                                                                                  aGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCCCGCGGACGGCGCTCGAGCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGATACCTCTGAAAACTGACCCCGAAGTc                                                  >  2:1191313/1‑150 (MQ=255)
                                                                                                    gtgtCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCg                                                <  2:510085/150‑1 (MQ=255)
                                                                                                       tCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCg                                                                                                                >  1:1336694/1‑84 (MQ=255)
                                                                                                       tCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCg                                                                                                                <  2:1336694/84‑1 (MQ=255)
                                                                                                         cGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTcc                                           <  2:1364968/150‑1 (MQ=255)
                                                                                                         cGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTcc                                           >  2:1511385/1‑150 (MQ=255)
                                                                                                         cGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTATCACGCCGGACACCTCTGAAAACTCACCCCGAAGTCCGGTTcc                                           >  2:501138/1‑150 (MQ=255)
                                                                                                         cGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACAGGCGGCGCTCGATCTCGGCTCAAGCCGTACACCTCTGAAAACTGACCCAGAAGTCCGGTTcc                                           >  2:924729/1‑150 (MQ=38)
                                                                                                          gCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTccg                                          >  1:2092733/1‑150 (MQ=255)
                                                                                                           ccTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGgcgc                                                                                                 >  1:916347/1‑94 (MQ=255)
                                                                                                           ccTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGgcgc                                                                                                 <  2:916347/94‑1 (MQ=255)
                                                                                                           ccTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCGGGACCCCTCTGAAAACTGACCCCGAAGTCCGGTTccgc                                         >  1:890066/1‑150 (MQ=255)
                                                                                                             ttAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGc                                                                                >  2:2266063/1‑109 (MQ=255)
                                                                                                             ttAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGc                                                                                <  1:2266063/109‑1 (MQ=255)
                                                                                                               aaGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCgccgcc                                     >  2:152880/1‑150 (MQ=255)
                                                                                                                    atTCTGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                                                              <  2:640711/41‑1 (MQ=255)
                                                                                                                    atTCTGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                                                              >  1:640711/1‑41 (MQ=255)
                                                                                                                      tCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTATCAAGCCGGACACCTCTGAAAACTGACTCCGAAGTCCGGTTCCGCCGCCCGGCAgg                              <  1:924729/150‑1 (MQ=255)
                                                                                                                         gcgcCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGt                           <  2:324485/150‑1 (MQ=255)
                                                                                                                           gcCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTgg                         <  1:2029848/150‑1 (MQ=255)
                                                                                                                            ccTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCATAGGCGCTGAATATCTTAAGGCGGACACCTGGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTgg                         <  1:1969500/150‑1 (MQ=17)
                                                                                                                             cTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGaaaa                                                                >  2:1805102/1‑109 (MQ=255)
                                                                                                                             cTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGaaaa                                                                <  1:1805102/109‑1 (MQ=255)
                                                                                                                               aTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTa                     <  1:152880/150‑1 (MQ=255)
                                                                                                                                tGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGc                                                                                <  1:2377273/90‑1 (MQ=255)
                                                                                                                                tGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGc                                                                                >  2:2377273/1‑90 (MQ=255)
                                                                                                                                tGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTag                    >  2:441107/1‑150 (MQ=255)
                                                                                                                                     tCGAGCGCCGCCCGCGCGGCGCTCGATCCCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAg               >  2:239581/1‑150 (MQ=255)
                                                                                                                                      cGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAgg              >  1:1784222/1‑150 (MQ=255)
                                                                                                                                       gAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCaa                                                                                  <  1:1772689/81‑1 (MQ=255)
                                                                                                                                       gAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCaa                                                                                  >  2:1772689/1‑81 (MQ=255)
                                                                                                                                       gAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCa                >  1:1693610/1‑147 (MQ=255)
                                                                                                                                       gAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCa                <  2:1693610/147‑1 (MQ=255)
                                                                                                                                       gAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGt             >  2:2346469/1‑150 (MQ=255)
                                                                                                                                       gAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGt             >  2:2233835/1‑150 (MQ=255)
                                                                                                                                         gcgcCGCCCGCGCGGCGCTCGATCTCATAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACctct                                                                     >  2:27128/1‑92 (MQ=255)
                                                                                                                                         gcgcCGCCCGCGCGGCGCTCGATCTCATAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACctct                                                                     <  1:27128/92‑1 (MQ=25)
                                                                                                                                           gccgcGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCg         >  2:1040585/1‑150 (MQ=255)
                                                                                                                                              gcgcgGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGaaa      <  2:1183584/150‑1 (MQ=255)
                                                                                                                                              gcgcgGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGaaa      >  1:502830/1‑150 (MQ=255)
                                                                                                                                                 cgggcggCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGgcgc                                                                                                 <  2:1101188/56‑1 (MQ=255)
                                                                                                                                                 cgggcggCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGgcgc                                                                                                 >  1:1101188/1‑56 (MQ=255)
                                                                                                                                                  gggcggCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGAAAGCCa  <  1:1040585/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TCTGGTAGCGGGTTGTGTGGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGAAAGCCA  >  NC_002947/1184447‑1184742

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