Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,695,585 (G)7→8 intergenic (‑42/‑216) purL ← / → mltF phosphoribosylformyl‑glycineamide synthetase/membrane‑bound lytic transglycosylase F, murein hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,695,5781.G89.9% 228.6 / 21.1 69intergenic (‑35/‑223)purL/mltFphosphoribosylformyl‑glycineamide synthetase/membrane‑bound lytic transglycosylase F, murein hydrolase
Reads supporting (aligned to +/- strand):  ref base . (4/3);  new base G (29/33);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-02

GGACATACTCGGCGTAAATATTGTGAACCGGGAGCCTGGCAGCCTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCAACGCGTGCCGCCAATGGCTGTAATAAGTTGCCATCTGGCGCAGGTTTACGCAAAATGCC  >  NC_000913/2695429‑2695723
                                                                                                                                                       |                                                                                                                                                   
ggACATACTCGGCGTAAATATTGTGAACCGGGAGCCTGGCAGCCTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCaa                                                                                                                                                      <  1:408921/150‑1 (MQ=255)
    aTACTCGGCGTAAATATTGTGAACCGGGAGCCTGGCAGCCTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAggggg                                                                                                                                               <  1:901999/151‑1 (MQ=255)
      aCTCGGCGTAAATATTGTGAACCGGGAGCCTGGCAGCCTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAggggtg                                                                                                                                              <  1:694565/150‑3 (MQ=255)
       cTCGGCGTAAATATTGTGAACCGGTAGCCTGGCAGCCTGAAAACGTGCCAGCAGTTTGTGGATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCGAggggga                                                                                                                                              <  1:748500/149‑2 (MQ=255)
         cGGCGTAAATATTGTGAACCGGGAGCCTGGCAGCCTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAgggggggga                                                                                                                                          >  1:32064/1‑148 (MQ=255)
         cGGCGTAAATATTGTGAACCGGGAGCCTGGCAGCCTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGgaa                                                                                                                                         <  1:403392/151‑1 (MQ=255)
          ggCGTAAATATTGTGAACCGGGAGCCTGGCAGCCTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGgaa                                                                                                                                         >  2:1080218/1‑150 (MQ=255)
           gCGTAAATATTGTGAACCGGGAGCCTGGCAGCCTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGgaaac                                                                                                                                       <  2:261873/151‑3 (MQ=255)
           gCGTAAATATTGTGAACCGGGAGCCTGGCAGCCTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGgaaac                                                                                                                                       >  2:598924/1‑149 (MQ=255)
             gTAAATATTGTGAACCGGGAGCCTGGCAGCCTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACgg                                                                                                                                      <  1:82106/151‑1 (MQ=255)
             gTAAATATTGTGAACCGGGAGCCTGGCAGCCTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACg                                                                                                                                       <  2:1268447/150‑1 (MQ=255)
               aaaTATTGTGAACCGGGAGCCTGGCAGCCTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGc                                                                                                                                    <  2:1052128/151‑1 (MQ=255)
                aaTATTGTGAACCGGGAGCCTGGCAGCCTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCg                                                                                                                                   >  2:1092659/1‑151 (MQ=255)
                aaTATTGTGAACCGGGAGCCTGGCAGCCTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCg                                                                                                                                   >  1:761608/1‑151 (MQ=255)
                         aaCCGGGAGCCTGGCAGCCTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATaa                                                                                                                          <  1:367374/151‑1 (MQ=255)
                          aCCGGGAGCCTGGCAGCCTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATaaa                                                                                                                         <  2:1227921/151‑1 (MQ=255)
                                aGCCTGGCAGCCTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAAGGGGGGGAAACGGGCGTCATTATAAAGAATCt                                                                                                                   >  1:953690/1‑151 (MQ=255)
                                aGCCTGGCAGCCTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCt                                                                                                                   <  1:491688/151‑1 (MQ=255)
                                        aGCCTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCt                                                                                                           <  1:276322/151‑1 (MQ=255)
                                        aGCCTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCt                                                                                                           <  2:1211159/151‑1 (MQ=255)
                                         gCCTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTg                                                                                                          >  2:683575/1‑151 (MQ=255)
                                          ccTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTg                                                                                                          <  1:418426/150‑1 (MQ=255)
                                          ccTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGa                                                                                                         <  2:930983/151‑1 (MQ=255)
                                             gAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGaa                                                                                                        <  2:1261152/149‑2 (MQ=255)
                                             gAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACg                                                                                                       <  2:477152/150‑1 (MQ=255)
                                              aaaaCGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGaa                                                                                                     <  2:361755/151‑1 (MQ=255)
                                               aaaCGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGaaa                                                                                                    <  2:121224/151‑1 (MQ=255)
                                                       cagcagTTTGTTGATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGc                                                                                            <  2:281575/151‑1 (MQ=255)
                                                             tttGTTGATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGa                                                                                       >  2:876063/1‑150 (MQ=255)
                                                                 ttgATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATa                                                                                   <  1:549520/150‑1 (MQ=255)
                                                                  tgATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATaa                                                                                  >  1:598860/1‑150 (MQ=255)
                                                                  tgATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATaa                                                                                  >  1:922737/1‑150 (MQ=255)
                                                                      tCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATc                                                                              <  2:1012859/150‑1 (MQ=255)
                                                                       cGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCAc                                                                            >  1:709436/1‑151 (MQ=255)
                                                                        ggAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCAcc                                                                           >  2:105695/1‑151 (MQ=255)
                                                                         gAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCAcc                                                                           >  2:679408/1‑150 (MQ=255)
                                                                            tGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCg                                                                       >  2:544905/1‑151 (MQ=255)
                                                                             gCCGACAGTGCATGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGTGAAACGGGAGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTGTGCGTGGAAATAAAATCACCATc                                                                        <  1:1291549/149‑1 (MQ=255)
                                                                              ccGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTg                                                                     >  1:706839/1‑151 (MQ=255)
                                                                                gccAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTg                                                                     >  2:1199862/3‑149 (MQ=255)
                                                                                   aGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAAtt                                                                 >  2:1074887/1‑150 (MQ=255)
                                                                                    gTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAg                                                               <  1:1279754/151‑1 (MQ=255)
                                                                                      gCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGc                                                              >  1:945883/1‑150 (MQ=255)
                                                                                      gCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGc                                                              >  2:796944/1‑150 (MQ=255)
                                                                                      gCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCa                                                             <  1:421319/151‑1 (MQ=255)
                                                                                       cAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCaa                                                            <  2:815877/151‑1 (MQ=255)
                                                                                       cAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCa                                                             >  1:413984/1‑150 (MQ=255)
                                                                                         ggCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCAAcg                                                          <  1:1092618/151‑1 (MQ=255)
                                                                                           cGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCAAcgc                                                         >  1:987362/1‑150 (MQ=255)
                                                                                           cGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCAAcgc                                                         >  2:305995/1‑150 (MQ=255)
                                                                                            gAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCAACGCGt                                                       <  1:827891/151‑1 (MQ=255)
                                                                                             aaCCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCAACGCGTg                                                      >  1:1014860/1‑151 (MQ=255)
                                                                                             aaCCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCAACGCGTg                                                      >  1:415877/1‑151 (MQ=255)
                                                                                                cACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCAACGCGTgccg                                                   <  2:321088/151‑1 (MQ=255)
                                                                                                 aCGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCAACGCGTgccgc                                                  >  1:375514/1‑151 (MQ=255)
                                                                                                     aGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCAACGCGTGCCGCCAAt                                              <  2:923404/151‑1 (MQ=255)
                                                                                                     aGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCAACGCGTGCCGCCAAt                                              >  2:244665/1‑151 (MQ=255)
                                                                                                         tttCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCAACGCGTGCCGCCAATGGCt                                          >  2:1230187/1‑151 (MQ=255)
                                                                                                            ccATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCAACGCGTGCCGCCAATGGCTg                                         >  1:642370/1‑149 (MQ=255)
                                                                                                              atcatAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCTCTGACGAAACCGTTTGCGTGGAAATAAAATCGCCATCGTGAATTAGCAACGGGTGCCGCCAATGGCTGtaa                                      <  1:344024/150‑1 (MQ=255)
                                                                                                                   aaGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCAACGCGTGCCGCCAATGGCTGTAATAAGtt                                >  1:33801/1‑151 (MQ=255)
                                                                                                                   aaGTCGCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCAACGCGTGCCGCCAATGGCTGTAATAAGtt                                <  1:1283893/151‑1 (MQ=255)
                                                                                                                    aGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCAACGCGTGCCGCCAATGGCTGTAATAAGTTg                               <  2:819482/151‑1 (MQ=255)
                                                                                                                          tcGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCAACGCGTGCCGCCAATGGCTGTAATAAGTTGCCATCt                         >  1:1014302/1‑151 (MQ=255)
                                                                                                                           cGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCAACGCGTGCCGCCAATGGCTGTAATAAGTTGCCATCTg                        <  1:897066/151‑1 (MQ=255)
                                                                                                                               ttctCAGCTTTCGGTGTACCCAAGCGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGAGTGGAAATAAAATCACCATCGTGAATTAGCAACCCGTCCCGCCAATGGCTGTAATAAGTTGCCATCTGgcgc                    >  1:465743/1‑151 (MQ=255)
                                                                                                                                  tAAGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCAACGCGTGCCGCCAATGGCTGTAATAAGTTGCCATCTGGCGCAgg                 >  2:735612/1‑151 (MQ=255)
                                                                                                                                   aaGCTTTCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCAACGCGTGCCGCCAATGGCTGTAATAAGTTGCCATCTGGCGCAGGt                <  2:1162400/151‑1 (MQ=255)
                                                                                                                                         tCGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCAACGCGTGCCGCCAATGGCTGTAATAAGTTGCCATCTGGCGCAGGTTTACGc          >  1:552586/1‑151 (MQ=255)
                                                                                                                                          cGGTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCAACGCGTGCCGCCAATGGCTGTAATAAGTTGCCATCTGGCGCAGGTTTACGCa         <  1:120300/151‑1 (MQ=255)
                                                                                                                                           ggTGTACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCAACGCGTGCCGCCAATGGCTGTAATAAGTTGCCATCTGGCGCAGGTTTACGCaa        <  1:1216504/151‑1 (MQ=255)
                                                                                                                                              gtACCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCAACGCGTGCCGCCAATGGCTGTAATAAGTTGCCATCTGGCGCAGGTTTACGCAAAAt     <  1:598893/151‑1 (MQ=255)
                                                                                                                                                aCCCAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCAACGCGTGCCGCCAATGGCTGTAATAAGTTGCCATCTGGCGCAGGTTTACGCAAAATGc   >  2:294594/1‑151 (MQ=255)
                                                                                                                                                 cccAAGGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCAACGCGTGCCGCCAATGGCTGTAATAAGTTGCCATCTGGCGCAGGTTTACGCAAAATGcc  <  2:1098488/151‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                   
GGACATACTCGGCGTAAATATTGTGAACCGGGAGCCTGGCAGCCTGAAAACGTGCCAGCAGTTTGTTGATTCGGAATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACCATCGTGAATTAGCAACGCGTGCCGCCAATGGCTGTAATAAGTTGCCATCTGGCGCAGGTTTACGCAAAATGCC  >  NC_000913/2695429‑2695723

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