Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,961,684 T→A 6.7% I175N (ATT→AAT)  lysA‑II → diaminopimelate decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,961,6840TA 6.7% 270.5 / 9.2 89I175N (ATT→AAT) lysA‑IIdiaminopimelate decarboxylase
Reads supporting (aligned to +/- strand):  ref base T (45/38);  new base A (3/3);  total (48/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 = 3.62e-01

ACCGACGAGCTGGAGCGCCTGCAAGTAGTCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCTGCTGGTGCTGGTCGACC  >  NC_002947/5961539‑5961827
                                                                                                                                                 |                                                                                                                                               
aCCGACGAGCTGGAGCGCCTGCAAGTAGTCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCCGGCCTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATtgcc                                                                                                                                             >  1:759507/1‑150 (MQ=255)
     cgaGCTGGAGCGCCTGCAAGTAGTCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATtgccgatgc                                                                                                                                        <  2:243187/150‑1 (MQ=255)
         cTGGAGCGCCTGCAAGTAGTCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGCCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAg                                                                                                                                    <  2:1045012/150‑1 (MQ=255)
            gAGCGCCTGCAAGTAGTCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGGCAAtggcgatgccgaactt                                                                                                                                 <  2:1555890/150‑5 (MQ=255)
              gcgcCTGCAAGTAGTCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCAt                                                                                                                               >  2:809179/1‑150 (MQ=255)
                 ccTGCAAGTAGTCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCta                                                                                                                            >  1:76060/1‑150 (MQ=255)
                  cTGCAAGTAGTCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCtac                                                                                                                           >  2:238169/1‑150 (MQ=255)
                  cTGCAAGTAGTCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCtac                                                                                                                           >  1:1859880/1‑150 (MQ=255)
                  cTGCAAGTAGTCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCtac                                                                                                                           <  1:1695271/150‑1 (MQ=255)
                   tGCAAGTAGTCGCCGCCGAAAAGCTCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCAGCGGTACCTGTTGAACGGCCCCAAAGAAAACAAGTTCGGCATCGCCATTGCCGTTGCCGAGGCCATCtact                                                                                                                          >  1:2089872/1‑149 (MQ=255)
                    gCAAGTAGTCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCtacgt                                                                                                                         >  2:2135751/1‑150 (MQ=255)
                     cAAGTAGTCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCtacgta                                                                                                                        >  1:831994/1‑150 (MQ=255)
                          agtCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACgcgc                                                                                                                   <  1:238169/150‑1 (MQ=255)
                            tCGCCGGCGAGATGGGCAAGGTTGCCCCGGTATCGCTGCGGGGCAACCCGGGCGTCGATGCCGGTACCCACCCGTAGATCTCCACCGGCCTAAAAGGAAACAAGTTCGTCATCGCCATTGCCGATGCCTAGACCATCTACGTACGcgccg                                                                                                                 >  2:650278/1‑150 (MQ=255)
                                  ccgAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGc                                                                                                           <  1:473837/150‑1 (MQ=255)
                                    gAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTa                                                                                                         >  2:2376690/1‑150 (MQ=255)
                                     aaaTGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTAc                                                                                                        <  2:939538/150‑1 (MQ=255)
                                     aaaTGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTAc                                                                                                        >  2:1753270/1‑150 (MQ=255)
                                         gtgCAAGGTTGCCCCGGTGTGGCCGCGGGTCAACCCGGACGTCGGTGCCCGTGCCGACACGTACCTCTCACCAGGCCTTAAACAAAACAAGTACGTCACCGCCATTCCCGTTGCCGAGGCCAACTCCGTACCCACCGCACAGCTACCGaa                                                                                                    >  2:788058/3‑150 (MQ=255)
                                          ggCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAAc                                                                                                   <  1:809179/150‑1 (MQ=255)
                                            cAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCAtt                                                                                                                                                >  1:2124960/1‑103 (MQ=255)
                                            cAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCAtt                                                                                                                                                <  2:2124960/103‑1 (MQ=255)
                                              accttGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGcgcc                                                                                                                  <  2:402986/128‑1 (MQ=255)
                                              accttGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGcgcc                                                                                                                  >  1:402986/4‑131 (MQ=255)
                                                gTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGaa                                                                                             <  2:418010/150‑1 (MQ=255)
                                                gTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGaa                                                                                             <  1:1225308/150‑1 (MQ=255)
                                                gTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGaa                                                                                             >  1:915365/1‑150 (MQ=255)
                                                    ccccGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACAACTCCACCGGCCTTAAAGAAAACAAGGTCGGCATCGCCCTTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTgg                                                                                         >  2:405892/1‑150 (MQ=255)
                                                      ccGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTc                                                                                       >  1:11532/1‑150 (MQ=255)
                                                         gtgtCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGgt                                                                                    <  2:1433036/150‑1 (MQ=255)
                                                           gtCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGgtgt                                                                                  <  2:1859880/150‑1 (MQ=255)
                                                           agagCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGcac                                                                                                              >  2:2290506/4‑122 (MQ=255)
                                                           agagCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGcac                                                                                                              <  1:2290506/119‑1 (MQ=255)
                                                            tCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTc                                                                                 >  2:1432285/1‑150 (MQ=255)
                                                                   gggTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGcc                                                                          >  1:2216997/1‑150 (MQ=255)
                                                                    ggTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATtgc                                                                                                                                              <  2:1087624/81‑1 (MQ=255)
                                                                    ggTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATtgc                                                                                                                                              >  1:1087624/1‑81 (MQ=255)
                                                                    ggTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGt                                                                                        >  1:684204/1‑135 (MQ=255)
                                                                    ggTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGt                                                                                        <  2:684204/135‑1 (MQ=255)
                                                                      tCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGt                                                                                           <  1:697229/130‑1 (MQ=255)
                                                                      tCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGt                                                                                           >  2:697229/1‑130 (MQ=255)
                                                                      tCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCcaca                                                                       >  1:1840093/1‑150 (MQ=255)
                                                                      tCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCcaca                                                                       >  1:234700/1‑150 (MQ=255)
                                                                           ccGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCATACCTACAGAACCTAGAAGTGGTCGGTGTCGACTGCCACATCGGc                                                                  >  2:1344946/1‑150 (MQ=255)
                                                                           ccGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATc                                                                     <  1:143972/147‑1 (MQ=255)
                                                                           ccGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATc                                                                     >  2:143972/1‑147 (MQ=255)
                                                                           acGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGc                                                                  >  2:1900797/2‑150 (MQ=255)
                                                                            cGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGc                                                                           <  1:1301238/140‑1 (MQ=255)
                                                                            cGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGc                                                                           >  2:1301238/1‑140 (MQ=255)
                                                                                     aTGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGa                                                        >  2:143196/1‑150 (MQ=255)
                                                                                             acccacccGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTc                                                <  2:234700/150‑1 (MQ=255)
                                                                                               ccacccGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGa                                              >  1:1114001/1‑150 (MQ=255)
                                                                                                cacccGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAg                                             <  1:1723113/150‑1 (MQ=255)
                                                                                                 acccGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGc                                            >  2:1888901/1‑150 (MQ=255)
                                                                                                 acccATACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGACGCCCTCTTCTTACGCGCCGCACAGCTACCTAACCTGCAACTGCTCGGTGTCGACTGCCCCATTGGCTCGCAGCTGACCACCGTCGAGc                                            >  2:982877/1‑150 (MQ=255)
                                                                                                    cGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCt                                         >  2:2404134/1‑150 (MQ=255)
                                                                                                      tACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTc                                       <  2:809538/150‑1 (MQ=255)
                                                                                                       aCATCTCCACCGGCCTTAAGGACCACAACGTCGACATCGCCATTGCCGATGCCGAGGCCATCTACTTCCTCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTcc                                      <  1:1251082/150‑1 (MQ=255)
                                                                                                            tcCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTc                                                                                 <  1:1959954/102‑1 (MQ=255)
                                                                                                            tcCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTc                                                                                 >  2:1959954/1‑102 (MQ=255)
                                                                                                                ccGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGccc                             <  1:143196/150‑1 (MQ=255)
                                                                                                                 cGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGAAATCTACGTAAGCGTACCACAGCTACTGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGAGCACCGGCGAGCCATTCCTCGACGCCCt                            >  2:405299/1‑150 (MQ=255)
                                                                                                                    ccTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGa                         <  2:1114001/150‑1 (MQ=255)
                                                                                                                     cTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGAc                        >  1:1037229/1‑150 (MQ=255)
                                                                                                                     cTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGAc                        <  1:1900797/150‑1 (MQ=255)
                                                                                                                      ttAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGAcc                       <  2:2089872/150‑1 (MQ=255)
                                                                                                                        aaaGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTAc                                                                                                        <  2:1792066/67‑1 (MQ=255)
                                                                                                                        aaaGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTAc                                                                                                        >  1:1792066/1‑67 (MQ=255)
                                                                                                                          aGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCct                   >  1:1334909/1‑150 (MQ=255)
                                                                                                                           gAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCt                                     <  2:95358/131‑1 (MQ=255)
                                                                                                                           gAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCt                                     >  1:95358/1‑131 (MQ=255)
                                                                                                                             aaaCAAGTTCGGCATTGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGctggct                <  2:1677873/150‑6 (MQ=255)
                                                                                                                             aaaCAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCctgct                <  2:2216997/150‑1 (MQ=255)
                                                                                                                              aaCAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCctgctg               <  1:2376690/150‑1 (MQ=255)
                                                                                                                              aaCAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCCCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCctgctg               >  2:827088/1‑150 (MQ=255)
                                                                                                                               aCAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCtgctgg              <  1:1432285/150‑1 (MQ=255)
                                                                                                                                  aGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCtgctggtgc           <  2:915365/150‑1 (MQ=255)
                                                                                                                                    ttCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCtgctggtgctg         <  1:2135751/150‑1 (MQ=255)
                                                                                                                                    ttCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCtgctggtgctg         >  2:1940220/1‑150 (MQ=255)
                                                                                                                                    ctcggcatcggcaATGGCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGgt                                                                                    >  1:549269/2‑75 (MQ=37)
                                                                                                                                    ctcggcatcggcaATGGCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGgt                                                                                    <  2:549269/74‑1 (MQ=37)
                                                                                                                                    ctcggcatcggcaATGGCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGAc                                 >  1:1037938/2‑126 (MQ=255)
                                                                                                                                    ctcggcatcggcaATGGCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGAc                                 <  2:1037938/125‑1 (MQ=255)
                                                                                                                                    ctcggcatcggcaATGGCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCtgctggtgctg         >  2:2255515/2‑150 (MQ=255)
                                                                                                                                     tCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCGCGACGCCCTCGACCGCCtgctggtgctgg        <  1:925384/150‑1 (MQ=255)
                                                                                                                                       ggCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCTGCTGGTGCTGGTc      <  1:405892/150‑1 (MQ=255)
                                                                                                                                          aTCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACGGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCTGCTGGTGCTGGTCgcc   >  2:1148155/1‑148 (MQ=255)
                                                                                                                                          aTCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCACAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCTGCTGGTGCTGGTCGAc   >  1:628472/1‑150 (MQ=255)
                                                                                                                                           tCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTc                                       >  2:2404407/1‑113 (MQ=255)
                                                                                                                                           tCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCTGCTGGTGCTGGTCGAcc  >  1:1382937/1‑150 (MQ=255)
                                                                                                                                           tCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCCTc                                       <  1:2404407/113‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
ACCGACGAGCTGGAGCGCCTGCAAGTAGTCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCTGCTGGTGCTGGTCGACC  >  NC_002947/5961539‑5961827

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: