Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913349,8010CT34.1% 92.7 / 68.1 88intergenic (+229/‑211)prpB/prpC2‑methylisocitrate lyase/2‑methylcitrate synthase
Reads supporting (aligned to +/- strand):  ref base C (26/32);  new base T (14/16);  total (40/48)
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.23e-01

CGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTG  >  NC_000913/349706‑349899
                                                                                               |                                                                                                  
cgccgcTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCt                                                                                                <  2:564040/100‑1 (MQ=33)
                              ggTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCAc                                                                  <  2:1197976/100‑1 (MQ=32)
                               gTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACt                                                                 >  1:681410/1‑100 (MQ=33)
                                tGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTg                                                                <  1:1652508/100‑1 (MQ=32)
                                 gCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCCCTGa                                                               >  2:817176/1‑100 (MQ=17)
                                  cACTGCCCTGTACGTCGGACAAGCCGCCTTTCGTCGCATCCGACAGCCGATGCCTGATGCGTCGCTGACGCTACTTATCAGGCCTACGGGGTGCACt                                                                 <  1:1462172/97‑1 (MQ=11)
                                  cACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGaa                                                              >  2:1167592/1‑100 (MQ=32)
                                   aCTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAAc                                                             <  2:938773/100‑1 (MQ=33)
                                   aCTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAAc                                                             <  2:397711/100‑1 (MQ=32)
                                   aCTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATCCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAAc                                                             >  1:1034439/1‑100 (MQ=17)
                                    cTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACt                                                            <  1:629489/100‑1 (MQ=33)
                                      gAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGt                                                          >  2:484185/1‑100 (MQ=32)
                                       atcTGTAGGTCGGATAAGACGGATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGt                                                          <  2:1062363/97‑1 (MQ=18)
                                       aaCTGTAGGTCGGATAAGACGGATAGCGGCGCATCCGAGAACCGCTGCCTGATGTGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACt                                                            <  1:796465/97‑1 (MQ=11)
                                       aaCTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTa                                                         <  1:753586/100‑1 (MQ=32)
                                             aGGTCTGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCg                                                    <  1:257187/99‑1 (MQ=17)
                                             aGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCgg                                                   <  1:1439797/100‑1 (MQ=33)
                                              ggTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTAAACTGTAGGTCGGa                                                  <  2:979524/100‑1 (MQ=18)
                                               gTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATACGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGAt                                                 >  2:1047766/1‑100 (MQ=17)
                                                  ggATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAg                                              >  2:1264045/1‑100 (MQ=32)
                                                   gATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGa                                             >  1:391238/1‑100 (MQ=33)
                                                   gATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGa                                             >  1:692249/1‑100 (MQ=32)
                                                   gATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGa                                             <  1:1037747/100‑1 (MQ=33)
                                                    aTAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGAc                                            <  2:136916/100‑1 (MQ=33)
                                                       aGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCATGCCTACGGGGTGCACTCAACTGTAGGTCGAATAAGACGCa                                         >  2:853038/1‑100 (MQ=255)
                                                       aGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCa                                         >  1:1383009/1‑100 (MQ=32)
                                                          cgtatggcttCGCATCCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAg                                      <  2:624313/91‑1 (MQ=1)
                                                          cGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATa                                       >  2:1697294/1‑99 (MQ=33)
                                                            cATAGCGTCGCATCCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCg                                    <  1:1793643/100‑1 (MQ=12)
                                                            cATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCg                                    >  2:1610603/1‑100 (MQ=33)
                                                            cATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCg                                    >  2:290866/1‑100 (MQ=32)
                                                            cATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCg                                    <  1:992293/100‑1 (MQ=33)
                                                             atggCGTCGCATCCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGt                                   <  2:1657052/97‑1 (MQ=255)
                                                               aGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCg                                 >  1:1672103/1‑100 (MQ=33)
                                                                gCGTCGCATCCGATAGCCGATGCCGGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGc                                <  2:1285482/100‑1 (MQ=11)
                                                                gCGTCGCATCCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGc                                >  1:237170/1‑100 (MQ=12)
                                                                 cGTCGCATTCGACAGCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCg                                 <  1:473477/98‑1 (MQ=12)
                                                                  gTCGCATCCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCAt                              <  2:420793/100‑1 (MQ=12)
                                                                  gTCGCATCCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCAt                              <  1:1576093/100‑1 (MQ=12)
                                                                  gTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCAt                              <  2:373196/100‑1 (MQ=32)
                                                                  gTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCAt                              >  1:1205794/1‑100 (MQ=32)
                                                                   tCGCATCCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATc                             <  2:1297948/100‑1 (MQ=255)
                                                                    cGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGGCGCATAGCGTCGCATcc                            >  1:415270/1‑100 (MQ=17)
                                                                     gCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCg                           >  2:292739/1‑100 (MQ=32)
                                                                      cATCCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCCCTGACCTGTAGGTCGGCTAAGACGCATAGCGTCGCATCCGa                          >  2:306588/1‑100 (MQ=11)
                                                                      cATCCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGa                          >  2:713806/1‑100 (MQ=255)
                                                                      cATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGa                          <  1:1734831/100‑1 (MQ=33)
                                                                       aTCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGAc                         <  1:506479/100‑1 (MQ=33)
                                                                          cGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAAc                      >  2:1037757/1‑100 (MQ=255)
                                                                          cGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCACATCCGACAAc                      >  2:531756/1‑100 (MQ=18)
                                                                           gACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGCTCGGATAAGACGCATAGCGTCGCATCCGACAAc                      >  1:958299/1‑99 (MQ=11)
                                                                           gACAAGCGATGTCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAAcc                     <  2:188703/100‑1 (MQ=16)
                                                                           gACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAAcc                     <  2:1260921/100‑1 (MQ=32)
                                                                            acagccGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCg                    <  1:830088/96‑1 (MQ=9)
                                                                            aCAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCg                    <  2:681414/100‑1 (MQ=32)
                                                                              agccGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCACATCCGACAACCGAt                  >  2:1439808/3‑100 (MQ=255)
                                                                               gccGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTATGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATg                 >  1:1440477/2‑100 (MQ=255)
                                                                               aCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAAc                      >  1:1621225/1‑95 (MQ=32)
                                                                               aCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGAt                  >  2:1549197/1‑99 (MQ=33)
                                                                               aCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAGATg                 >  2:677727/1‑100 (MQ=18)
                                                                                ccGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGc                >  2:239077/1‑100 (MQ=18)
                                                                                ccGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGc                <  2:480097/100‑1 (MQ=17)
                                                                                ccGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGc                <  2:985707/100‑1 (MQ=17)
                                                                                ccGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGc                >  1:1830/1‑100 (MQ=18)
                                                                                ccGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAGGACGCATAGCGTCGCATCCGACAAGCGAt                  >  1:1591335/1‑98 (MQ=16)
                                                                                ccGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGc                >  2:1684727/1‑100 (MQ=33)
                                                                                ccGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGc                <  2:237171/100‑1 (MQ=32)
                                                                                 cGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGcc               <  1:1403524/100‑1 (MQ=17)
                                                                                 cGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGTTAAGACGCATAGCGTCGCATCCGACAACCGATGcc               <  1:888560/100‑1 (MQ=17)
                                                                                 cGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGcc               <  2:1776503/100‑1 (MQ=33)
                                                                                 cGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGcc               <  2:1282035/100‑1 (MQ=33)
                                                                                 cGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGcc               <  1:239076/100‑1 (MQ=32)
                                                                                  gATGCCTGATGCGCCGCTGACGCGACTTATCAGGTCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCt              <  2:166085/100‑1 (MQ=17)
                                                                                   aTGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTg             <  2:504421/100‑1 (MQ=33)
                                                                                     gCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGAt           <  2:1669720/100‑1 (MQ=18)
                                                                                     gCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGAt           >  2:75068/1‑100 (MQ=18)
                                                                                     gCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGAt           <  1:1809365/100‑1 (MQ=33)
                                                                                      ccTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATg          <  1:1152401/100‑1 (MQ=18)
                                                                                      ccTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATg          <  1:713802/100‑1 (MQ=33)
                                                                                       cTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATg          >  1:135488/1‑99 (MQ=17)
                                                                                        tGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATgcg        >  2:1153112/1‑100 (MQ=17)
                                                                                         gATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATgcgc       <  2:1391286/100‑1 (MQ=17)
                                                                                         gATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATgcg        >  1:136916/1‑99 (MQ=18)
                                                                                         gATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATgc         >  2:509788/1‑98 (MQ=33)
                                                                                          aTGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGcgcc      <  1:540808/100‑1 (MQ=18)
                                                                                          aTGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGcgcc      <  1:1436759/100‑1 (MQ=32)
                                                                                           tGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATCggccg     >  1:1776490/1‑100 (MQ=255)
                                                                                           tGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGcgccg     >  2:315753/1‑100 (MQ=33)
                                                                                            gcgtCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGcgccgc    >  1:1319784/1‑100 (MQ=18)
                                                                                            gcgcCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGcgccgc    <  1:1591800/100‑1 (MQ=33)
                                                                                            gcgcCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGcgccgc    >  2:868057/1‑100 (MQ=33)
                                                                                             cgccgcTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCt   >  2:1212775/1‑100 (MQ=32)
                                                                                              gccgcTGACGCGACTTATCAGGCCTAAGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTg  >  2:1103720/1‑100 (MQ=18)
                                                                                               |                                                                                                  
CGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTG  >  NC_000913/349706‑349899

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: