Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,961,681 C→G 7.0% A174G (GCC→GGC)  lysA‑II → diaminopimelate decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,961,6810CG 7.0% 214.1 / 6.0 70A174G (GCC→GGC) lysA‑IIdiaminopimelate decarboxylase
Reads supporting (aligned to +/- strand):  ref base C (34/31);  new base G (2/3);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

CTGGAGCGCCTGCAAGTAGTCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCTGCTGGTGCTGGTCGAC  >  NC_002947/5961548‑5961826
                                                                                                                                     |                                                                                                                                                 
cTGGAGCGCCTGCAAGTAGTCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAg                                                                                                                                   <  2:1814110/150‑1 (MQ=255)
           gCAAGTAGTCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGAACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCtacgt                                                                                                                        >  2:560887/1‑150 (MQ=255)
                 agtCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACgcgc                                                                                                                  <  1:635351/150‑1 (MQ=255)
                   tCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTATCCACCCGTACATCTCCACCGACCTTAAAGAAAACAAGTTCGGCATCGCCCTTTCGTAGCCCGAGGCCATCTACGTACGcgccg                                                                                                                <  1:1020760/150‑1 (MQ=255)
                   tCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGCTCGGCATCGGCAAtggcgatgccgaac                                                                                                                                  <  2:924099/132‑3 (MQ=255)
                   tCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGCTCGGCATCGGCAAtggcgatgccgaac                                                                                                                                  >  1:924099/1‑130 (MQ=255)
                       cgccgAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGcaca                                                                                                            <  2:793431/150‑1 (MQ=255)
                           gAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATtgc                                                                                                                                             <  2:639623/113‑1 (MQ=255)
                           gAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATtgc                                                                                                                                             >  1:639623/1‑113 (MQ=255)
                                 ggCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAAc                                                                                                  <  2:1843731/150‑1 (MQ=255)
                                       gTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGGTGGCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGaa                                                                                            <  1:509773/150‑1 (MQ=255)
                                            cccGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGt                                                                                       <  1:915464/150‑1 (MQ=255)
                                             ccGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTc                                                                                      >  1:1300735/1‑150 (MQ=255)
                                             ccGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTc                                                                                      <  1:1139102/150‑1 (MQ=255)
                                              cGGTGTCGCTGCGGGTCAACCCGGCCGACGGTGAAGGTCACCAGCTCTACGTTACCGGTGCCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGAGCCGCCCAGCTACCGTACCGGGAAGTGGTCg                                                                                     <  2:426514/150‑1 (MQ=255)
                                              cGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCg                                                                                     <  2:1716085/150‑1 (MQ=255)
                                                       tGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACt                                                                            >  2:1246735/1‑150 (MQ=255)
                                                              cAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGc                                                                          <  2:1159893/145‑1 (MQ=255)
                                                              cAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGc                                                                          >  1:1159893/1‑145 (MQ=255)
                                                                 cccGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCGTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGa                                                                              >  2:1405681/1‑138 (MQ=255)
                                                                 cccGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGa                                                                              <  1:1405681/138‑1 (MQ=255)
                                                                    ggACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTAcc                                                                                                      >  1:272260/1‑111 (MQ=255)
                                                                    ggACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTAcc                                                                                                      <  2:272260/111‑1 (MQ=255)
                                                                        gTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTc                                                                                      <  2:416707/123‑1 (MQ=255)
                                                                        gTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTc                                                                                      >  1:416707/1‑123 (MQ=255)
                                                                        gTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAg                                                           >  2:1916627/1‑150 (MQ=255)
                                                                            aTGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGa                                                       <  1:1443246/150‑1 (MQ=255)
                                                                            aTGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGa                                                       >  2:667671/1‑150 (MQ=255)
                                                                                cGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCt                                                         <  2:1345827/144‑1 (MQ=255)
                                                                                cGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCt                                                         >  1:1345827/1‑144 (MQ=255)
                                                                                     cccacccGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCg                                              >  2:1598112/1‑150 (MQ=255)
                                                                                      ccacccGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGgtg                                                                                  <  1:185561/113‑1 (MQ=255)
                                                                                      ccacccGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGgtg                                                                                  >  2:185561/1‑113 (MQ=255)
                                                                                        acccGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGt                                                <  2:1579228/145‑1 (MQ=255)
                                                                                        acccGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGt                                                >  1:1579228/1‑145 (MQ=255)
                                                                                         cccGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACGGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGcc                                          >  1:589924/1‑150 (MQ=255)
                                                                                          ccGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGaa                                                                                            >  2:208840/1‑99 (MQ=255)
                                                                                          ccGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGaa                                                                                            <  1:208840/99‑1 (MQ=255)
                                                                                          ccGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGccc                                         >  2:791494/1‑150 (MQ=255)
                                                                                                aTCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTc                                   <  1:667671/150‑1 (MQ=255)
                                                                                                 tctcCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTTGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCg                                  >  1:348643/1‑150 (MQ=255)
                                                                                                 tctcCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGccc                                         <  2:30393/143‑1 (MQ=255)
                                                                                                 tctcCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGccc                                         >  1:30393/1‑143 (MQ=255)
                                                                                                  ctcCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCg                                                                                     <  1:1015726/98‑1 (MQ=255)
                                                                                                  ctcCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCg                                                                                     >  2:1015726/1‑98 (MQ=255)
                                                                                                   tcCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTCCGCGCCGCACAGCTTCCGAACCTGTAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGAc                                <  1:233872/150‑1 (MQ=255)
                                                                                                         ggCCTTAAAGAAAACAAGTTCGGGATGGCCATTGCCGATTGCGAGGCCATCTACGTACGCGCCGAACCGCTAACGAACCTGGAAGTGGTCGGTGTCGATCTCCACAACGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTc                          <  2:1726873/150‑1 (MQ=255)
                                                                                                         ggCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTc                          >  1:1730660/1‑150 (MQ=255)
                                                                                                             ttAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGAcc                      >  1:1391033/1‑150 (MQ=255)
                                                                                                              tAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCg                     >  2:1106375/1‑150 (MQ=255)
                                                                                                              tAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGAACCGTTAGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCg                     >  1:1573925/1‑150 (MQ=255)
                                                                                                               aaaGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGc                    <  2:1730660/150‑1 (MQ=255)
                                                                                                                  gAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCctg                 >  1:72221/1‑150 (MQ=255)
                                                                                                                  gAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCctg                 >  2:493329/1‑150 (MQ=255)
                                                                                                                   aaaaCAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCctgc                >  1:1267590/1‑150 (MQ=255)
                                                                                                                     aaCAAGTTCGGCATCGCCATTGCCGATGCCGAGGCGATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCt                           <  2:1300583/137‑1 (MQ=255)
                                                                                                                     aaCAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCt                           >  1:1300583/1‑137 (MQ=255)
                                                                                                                       cAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCt                                        >  1:599451/1‑122 (MQ=255)
                                                                                                                       cAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCt                                        <  2:599451/122‑1 (MQ=255)
                                                                                                                       cAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCtgctggt            >  2:1673238/1‑150 (MQ=255)
                                                                                                                        aaGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCtgctggtg           >  2:1203928/1‑150 (MQ=255)
                                                                                                                            tctgCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCtgctggtgctgg       <  1:212310/147‑1 (MQ=255)
                                                                                                                            tcggcatcggcaATGGCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCtgctggtgctgg       >  2:1714853/1‑150 (MQ=255)
                                                                                                                            tCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCtgctggtgctgg       <  1:131705/150‑1 (MQ=255)
                                                                                                                             cGGCATCGGCAATGGCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTg                                                                                         <  1:141125/67‑1 (MQ=38)
                                                                                                                             cGGCATCGGCAATGGCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTg                                                                                         >  2:141125/1‑67 (MQ=37)
                                                                                                                             cGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGa                                                       >  1:879229/1‑101 (MQ=255)
                                                                                                                             cGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGa                                                       <  2:879229/101‑1 (MQ=255)
                                                                                                                              ggCATCGGCAATGGCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCTGCTGGTGCTGGTc     <  1:1714853/150‑1 (MQ=255)
                                                                                                                                cATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTTCTCGACTCCCTCGACCGCCTGCTGGGGCTGGTAGa   >  2:212310/1‑150 (MQ=255)
                                                                                                                                 aTCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCTGCTGGTGCTGGTCGAc  >  1:1857054/1‑150 (MQ=255)
                                                                                                                                     |                                                                                                                                                 
CTGGAGCGCCTGCAAGTAGTCGCCGCCGAAATGGGCAAGGTTGCCCCGGTGTCGCTGCGGGTCAACCCGGACGTCGATGCCGGTACCCACCCGTACATCTCCACCGGCCTTAAAGAAAACAAGTTCGGCATCGCCATTGCCGATGCCGAGGCCATCTACGTACGCGCCGCACAGCTACCGAACCTGGAAGTGGTCGGTGTCGACTGCCACATCGGCTCGCAGCTGACCACCGTCGAGCCCTTCCTCGACGCCCTCGACCGCCTGCTGGTGCTGGTCGAC  >  NC_002947/5961548‑5961826

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: