Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,364,888 C→G 6.7% M199I (ATG→ATC lysA‑I ← diaminopimelate decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,364,8880CG 6.7% 184.0 / 4.8 60M199I (ATG→ATClysA‑Idiaminopimelate decarboxylase
Reads supporting (aligned to +/- strand):  ref base C (25/31);  new base G (2/2);  total (27/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.44e-01

GCTTCGATCTCCTTGCGCGCCACCGCCCAGGTCTGCGCATAACGCTGCAGGTCGACCGGTTTGTCGGTGCTGCGGTAAGGGGTGGAAAGGCCGCCGCCGATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTGGTAGTCCACCCCGGAGCCAATGTGCATGTGCACGCCCACCAGGTGCAGGCCATGGGCCTTGATGCACGCCAGCGCCTCGTCCAGCTCTTCGTGCCAGATGCCAT  >  NC_002947/2364742‑2365021
                                                                                                                                                  |                                                                                                                                       
gCTTCGATCTCCTTGCGCGCCACCGCCCAGGTCTGCGCATAACGCTGCAGGTCGACCGGTTTGTCGGTGCTGCGGTAAGGGGTGGAAAGGCCGCCGCCGATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGatcatc                                                                                                                                      <  1:217862/150‑1 (MQ=255)
gCTTCGATCTCCTTGCGCGCCACCGCCCAGGTCTGCGCATAACGCTGCAGGTCGACCGGTTTGTCGGTGCTGCGGTAAGGGGTGGAAAGGCCGCCGCCGATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATAAGCTCGatcatc                                                                                                                                      <  1:51079/150‑1 (MQ=255)
    cGATCTCCTTGCGCGCCACCGCCCAGGTCTGCGCATAACGCTGCAGGTCGACCGGTTTGTCGGTGCTGCGGTAAGGGGTGGAAAGGCCGCCGCCGATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAga                                                                                                                                  <  2:421978/150‑2 (MQ=255)
     gATCTCCTTGCGCGCCACCGCCCAGGTCTGCGCATAACGCTGCAGGTCGACCGGTTTGTCGGTGCTGCGGTAAGGGGTGGAAAGGCCGCCGCCGATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCg                                                                                                                                 >  2:673208/1‑150 (MQ=255)
      aTCTCCTTGCGCGCCACCGCCCAGGTCTGCGCATAACGCTGCAGGTCGACCGGTTTGTCGGTGCTGCGGTAAGGGGTGGAAAGGCCGCCGCCGATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCgg                                                                                                                                >  2:342530/1‑150 (MQ=255)
           cTTGCGCGCCACCGCCCAGGTCTGCGCATAACGCTGCAGGTCGACCGGTTTGTCGGTGCTGCGGTAAGGGGTGGAAAGGCCGCCGCCGATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCAc                                                                                                                           >  2:754278/1‑150 (MQ=255)
           cTTGCGCGCCACCGCCCAGGTCTGCGCATAACGCTGCAGGTCGACCGGTTTGTCGGTGCTGCGGTAAGGGGTGGAAAGGCCGCCGCCGATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCAc                                                                                                                           <  1:260913/150‑1 (MQ=255)
                  gcCACCGCCCAGGTCTGCGCATAACGCTGCAGGTCGACCGGTTTGTCGGTGCTGCGGTAAGGGGTGGAAAGGCCGCCGCCGATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTcc                                                                                                                    >  1:331652/1‑150 (MQ=255)
                   ccACCGCCCAGGTCTGCGCATAACGCTGCAGGTCGACCGGTTTGTCGGTGCTGCGGTAAGGGGTGGAAAGGCCGCCGCCGATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCa                                                                                                                   >  1:929496/1‑150 (MQ=255)
                    cACCGCCCAGGTCTGCGCATAACGCTGCAGGTCGACCGGTTTGTCGGTGCTGCGGTAAGGGGTGGAAAGGCCGCCGCCGATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAg                                                                                                                  >  1:676482/1‑150 (MQ=255)
                     aCCGCCCAGGTCTGCGCATAACGCTGCAGGTCGACCGGTTTGTCGGTGCTGCGGTAAGGGGTGGAAAGGCCGCCGCCGATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAg‑‑g                                                                                                               <  1:342530/150‑1 (MQ=255)
                         cccAGGTCTGCGCATAACGCTGCAGGTCGACCGGTTTGTCGGTGCTGCGGTAAGGGGTGGAAAGGCCGCCGCCGATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCt                                                                                                           <  1:696131/150‑1 (MQ=255)
                         cccAGGTCTGCGCATAACGCTGCAGGTCGACCGGTTTGTCGGTGCTGCGGTAAGGGGTGGAAAGGCCGCCGCCGATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCt                                                                                                           <  2:676482/150‑1 (MQ=255)
                          ccAGGTCTGCGCATAACGCTGCAGGTCGACCGGTTTGTCGGTGCTGCGGTAAGGGGTGGAAAGGCCGCCGCCGATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTg                                                                                                          <  1:174829/150‑1 (MQ=255)
                          ccAGGTCTGCGCATAACGCTGCAGGTCGACCGGTTTGTCGGTGCTGCGGTAAGGGGTGGAAAGGCCGCCGCCGATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTg                                                                                                          <  1:98747/150‑1 (MQ=255)
                           cAGGTCTGCGCATAACGCTGCAGGTCGACCGGTTTGTCGGTGCTGCGGTAAGGGGTGGAAAGGCCGCCGCCGATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTgg                                                                                                         <  1:673208/150‑1 (MQ=255)
                                cTGCGCATAACGCTGCAGGTCGACCGGTTTGTCGGTGCTGCGGTAAGGGGTGGAAAGGCCGCCGCCGATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGTCCACCTGTTCCAG‑‑GTGCTGGTAGTc                                                                                                    <  2:991260/150‑1 (MQ=255)
                                 tGCGCATAACGCTGCAGGTCGACCGGTTTGTCGGTGCTGCGGTAAGGGGTGGAAAGGCCGCCGCCGATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTGGTAGTcc                                                                                                   >  1:672296/1‑150 (MQ=255)
                                  gcgcATAACGCTGCAGGTCGACCGGTTTGTCGGTGCTGCGGTAAGGGGTGGAAAGGCCGCCGCCGATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTGGTAGTCCa                                                                                                  <  2:179771/150‑1 (MQ=255)
                                   cgcATAACGCTGCAGGTCGACCGGTTTGTCGGTGCTGCGGTAAGGGGTGGAAAGGCCGCCGCCGATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTGGTAGTCCAc                                                                                                 >  1:1034206/1‑150 (MQ=255)
                                     cATAACGCTGCAGGTCGACCGGTTTGTCGGTGCTGCGGTAAGGGGTGGAAAGGCCGCCGCCGATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTGGTAGTCCAccc                                                                                               <  1:609077/150‑1 (MQ=255)
                                      aTAACGCTGCAGGTCGACCGGTTTGTCGGTGCTGCGGTAAGGGGTGGAAAGGCCGCCGCCGATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATGATCGAGctgat                                                                                                                              >  1:928399/1‑116 (MQ=255)
                                      aTAACGCTGCAGGTCGACCGGTTTGTCGGTGCTGCGGTAAGGGGTGGAAAGGCCGCCGCCGATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATGATCGAGctgat                                                                                                                              <  2:928399/120‑5 (MQ=255)
                                         aCGCTGCAGGTCGACCGGTTTGTCGGTGCTGCGGTAAGGGGTGGAAAGGCCGCCGCCGATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTGGTAGTCCACCCCGGa                                                                                           <  1:665118/150‑1 (MQ=255)
                                                aGGTCGACCGGTTTGTCGGTGCTGCGGTAAGGGGTGGAAAGGCCGCCGCCGATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTGGTAGTCCACCCCGGAGCCAatg                                                                                    <  2:690114/150‑1 (MQ=255)
                                                 ggTCGACCGGTTTGTCGGTGCTGCGGTAAGGGGTGGAAAGGCCGCCGCCGATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTGGTAGTCCACCCCGGAGCCAatgt                                                                                   <  2:1014204/150‑1 (MQ=255)
                                                 ggTCGACCGGTTTGTCGGTGCTGCGGTAAGGGGTGGAAAGGCCGCCGCCGATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTGGTAGTCCACCCCGGAGCCAatgt                                                                                   <  2:936822/150‑1 (MQ=255)
                                                   gggACGGGTTTGTCGGTGCTGCGGTAAGGGGTGGAAAGGCCGCCGCCGATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGGGCGGGCCACCTGTTCCAG‑‑GTGCTGGTAGTCCACCCCGGAGCCAatgtgc                                                                                 <  2:825700/148‑1 (MQ=255)
                                                     gACCGGTTTGTCGGTGCTGCGGTAAGGGGTGGAAAGGCCGCCGCCGATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTGGTAGTCCACCCCGGAGCCAatgtgcat                                                                               >  2:93107/1‑150 (MQ=255)
                                                        cGGTTTGTCGGTGCTGCGGTAAGGGGTGGAAAGGCCGCCGCCGATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTGGTAGTCCACCCCGGAGCCAatgtgcatgtg                                                                            <  1:429253/150‑1 (MQ=255)
                                                                   tgctgcGGTAAGGGGTGGAAAGGCCGCCGCCGATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTGGTAGTCCACCCCGGAGCCAATGTGCATGTGCACGCccacca                                                                 <  2:40219/150‑1 (MQ=255)
                                                                   tgctgcGGTAAGGGGTGGAAAGGCCGCCGCCGATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTGGTAGTCCACCCCGGAGCCAATGTGCATGTGCACGCccacca                                                                 <  1:406449/150‑1 (MQ=255)
                                                                    gctgcGGTAAGGGGTGGAAAGGCCGCCGCCGATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTGGTAGTCCACCCCGGAGCCAATGTGCATGTGCACGCCCACCAg                                                                >  2:668227/1‑150 (MQ=255)
                                                                          gTAAGGGGTGGAAAGGCCGCCGCCGATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTGGTAGTCCACCCCGGAGCCAATGTGCATGTGCACGCCCACCAGGTGCAg                                                          <  2:812935/150‑1 (MQ=255)
                                                                               gggTGGAAAGGCCGCCGCCGATGGAGAACGCATCGCTGTCCACGCCCAGGGGTCCGATGAGATCGATCATCTAGCGG‑‑CACCTGCTGTAGTCGATTTGTTATTCCACCCCGGAGCCAATGTGCATGTGCACGACCAACAGTTGCAGGCCAt                                                     >  1:252407/1‑150 (MQ=255)
                                                                                     aaaGGCCGCCGCCGATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTGGTAGTCCACCCCGGAGCCAATGTGCATGTGCACGCCCACCAGGTGCAGGCCATGGGCCt                                               <  1:93107/150‑1 (MQ=255)
                                                                                     aaaGGCCGCCGCCGATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTGGTAGTCCACCCCGGAGCCAATGTGCATGTGCACGCCCACCAGGTGCAGGCCATGGGCCt                                               >  1:824449/1‑150 (MQ=255)
                                                                                       aGGCCGCCGCCGATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTGGTAGTCCACCCCGGAGCCAATGTGCATGTGCACGCCCACCAGGTGCAGGCCATGGGCCTTg                                             >  1:390781/1‑150 (MQ=255)
                                                                                       aGGCCGCCGCCGATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTGGTAGTCCACCCCGGAGCCAATGTGCATGTGCACGCCCACCAGGTGCAGGCCATGGGCCTTg                                             <  2:331652/150‑1 (MQ=255)
                                                                                       aGGCCGCCGCCGATGGAGAACGACTGGTAGGGCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTGGTAGTCCACCCCGGAGCCAATGTGCATGTGCACGCCCACCAGGTGCAGGCCATGGGCCTTg                                             >  1:578555/1‑150 (MQ=255)
                                                                                         gccgccgccGATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTGGTAGTCCACCCCGGAGCCAATGTGCATGTGCACGCCCACCAGGTGCAGGCCATGGGCCTTGAt                                           <  1:754278/150‑1 (MQ=255)
                                                                                         gccgccgccGATGGAGAACGCCTCGAAGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTGGTAGTCCACCCCGGAGCCAATGTGCATGTGCACGCCCACCAGGTGCAGGCCATGGGCCTTGAt                                           <  2:1127485/150‑1 (MQ=255)
                                                                                           cgccgccgATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTGGTAGTCCACCCCGGAGCCAATGTGCATGTGCACGCCCACCAGGTGCAGGCCATGGGCCTTGATGc                                         >  2:1032756/1‑150 (MQ=255)
                                                                                            gccgccgATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTGGTAGTCCACCCCGGAGCCAATGTGCATGTGCACGCCCACCAGGTGCAGGCCATGGGCCTTGATGCa                                        >  1:341224/1‑150 (MQ=255)
                                                                                             ccgccgATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTGGTAGTCCACCCCGGAGCCAatgtgcat                                                                               <  2:548229/110‑1 (MQ=255)
                                                                                             ccgccgATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTGGTAGTCCACCCCGGAGCCAatgtgcat                                                                               >  1:548229/1‑110 (MQ=255)
                                                                                             ccgccgATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTGGTAGTCCACCCCGGAGCCAATGTGCATGTGCACGCCCACCAGGTGCAGGCCATGGGCCTTGATGCAc                                       <  1:771690/150‑1 (MQ=255)
                                                                                              cgccgATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTGGTAGTCCACCCCGGAGCCAATGTGCATGTGCACGcc                                                                      <  1:590501/118‑1 (MQ=255)
                                                                                              cgccgATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTGGTAGTCCACCCCGGAGCCAATGTGCATGTGCACGcc                                                                      >  2:590501/1‑118 (MQ=255)
                                                                                               gccgATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTGGTAGTCCACCCCGGAGCCAATGTGCATGTGCACGCCCACCAGGTGCAGGCCATGGGCCTTGATGCACGc                                     >  2:1022140/1‑150 (MQ=255)
                                                                                                  gATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTGGTAGTCCACCCCGGAGCCAATGTGCATGTGCACGCCCACCAGGTGCAGGCCATGGGCCTTGATg                                          <  2:907054/142‑1 (MQ=255)
                                                                                                  gATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTGGTAGTCCACCCCGGAGCCAATGTGCATGTGCACGCCCACCAGGTGCAGGCCATGGGCCTTGATg                                          >  1:907054/1‑142 (MQ=255)
                                                                                                            gCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTGGTAGTCCACCCCGGAGCCAATGTGCATGTGCACGCCCACCAGGTGCAGGCCATGGGCCTTGATGCACGCCAGCGCCTCGTcc                        >  1:320612/1‑150 (MQ=255)
                                                                                                            gCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTGGTAGTCCACCCCGGAGCCAATGTGCATGTGCACGCCCACCAGGTGCAGGCCATGGGCCTTGATGCACGCCAGCGCCTCGTcc                        >  2:837185/1‑150 (MQ=255)
                                                                                                                 gATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTGGTAGTCCACCCCGGAGCCAATGTGCATGTGCACGCCCACCAGGTGCAGGCCATGGGCCTTGATGCACGCCAGCGCCTCGTCCAGctc                   >  1:526080/1‑150 (MQ=255)
                                                                                                                       cACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTGGTAGTCCACCCCGGAGCCAATGTGCATGTGCACGCCCACCAGGTGCAGGCCATGGGCCTTGATGCACGCCAGCGCCTCGTCCAGCTCTTCGTg             >  1:1009975/1‑150 (MQ=255)
                                                                                                                        aCGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTGGTAGTCCACCCCGGAGCCAATGTGCATGTGCACGCCCACCAGGTGCAGGCCATGGGCCTTGATGCACGCCAGCGCCTCGTCCAGCTCTTCGTGc            <  2:821626/150‑1 (MQ=255)
                                                                                                                          gCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTGGTAGTCCACCCCGGAGCCAATGTGCATGTGCACGCCCACCAGGTGCAGGCCATGGGCCTTGATGCACGCCAGCGCCTCGTCCAGCTCTTCGTGCCa          <  2:1009975/150‑1 (MQ=255)
                                                                                                                            ccAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTGGTAGTCCACCCCGGAGCCAATGTGCATGTGCACGCCCACCAGGTGCAGGCCATGGGCCTTGATGCACGCCAGCGCCTCGTCCAGCTCTTCGTGCCAGa        >  2:1055679/1‑150 (MQ=255)
                                                                                                                                  ggcCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTGGTAGTCCACCCCGGAGCCAATGTGCATGTGCACGCCCACCAGGTGCAGGCCATGGGCCTTGATGCACGCCAGCGCCTCGTCCAGCTCTTCGTGCCAGATGCCAt  <  1:1032756/150‑1 (MQ=255)
                                                                                                                                      ggcccgCTCGATGATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTGGTAGTCCACCCCGGAGCCAATGTGCATGTGCACGCCCACCAGGTGCAGGCCATGGGCCTTGATGCAc                                       >  2:644772/6‑109 (MQ=255)
                                                                                                                                      ggcccgCTCGATGATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTGGTAGTCCACCCCGGAGCCAATGTGCATGTGCACGCCCACCAGGTGCAGGCCATGGGCCTTGATGCAc                                       <  1:644772/104‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                       
GCTTCGATCTCCTTGCGCGCCACCGCCCAGGTCTGCGCATAACGCTGCAGGTCGACCGGTTTGTCGGTGCTGCGGTAAGGGGTGGAAAGGCCGCCGCCGATGGAGAACGCCTCGATGTCCACGCCCAGGCGGCCGATCAGCTCGATCATCGAGCGGGCCACCTGTTCCAG‑‑GTGCTGGTAGTCCACCCCGGAGCCAATGTGCATGTGCACGCCCACCAGGTGCAGGCCATGGGCCTTGATGCACGCCAGCGCCTCGTCCAGCTCTTCGTGCCAGATGCCAT  >  NC_002947/2364742‑2365021

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