Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 654,560 A→C 16.4% intergenic (+18/+23) citB → / ← dcuC response regulator in two‑component regulatory system with CitA/anaerobic C4‑dicarboxylate transport

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913654,5600AC16.4% 468.6 / 82.7 194intergenic (+18/+23)citB/dcuCresponse regulator in two‑component regulatory system with CitA/anaerobic C4‑dicarboxylate transport
Reads supporting (aligned to +/- strand):  ref base A (67/93);  new base C (18/14);  total (85/109)
Fisher's exact test for biased strand distribution p-value = 1.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

Only 100 of 202 total aligned reads displayed.

TGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAA  >  NC_000913/654462‑654659
                                                                                                  |                                                                                                   
tGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGaa                                                                                                    <  1:2540510/100‑1 (MQ=255)
 gcgcCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGaaa                                                                                                   >  1:497575/1‑100 (MQ=255)
   gcCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACt                                                                                                 <  1:1248611/100‑1 (MQ=255)
   gcCAGCCGCCATCTGATTATTGCTGAAATTGTTCACAGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACt                                                                                                 >  2:376553/1‑100 (MQ=255)
     cAGCCGCCATCTAATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTg                                                                                               >  1:2709194/1‑100 (MQ=255)
       gccgccATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGcc                                                                                             <  2:604288/100‑1 (MQ=255)
       gccgccATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGcc                                                                                             >  2:3539187/1‑100 (MQ=255)
       gccgccATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGcc                                                                                             >  1:128086/1‑100 (MQ=255)
       gccgccATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGcc                                                                                             <  2:2448113/100‑1 (MQ=255)
          gccATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTcc                                                                                          >  2:2504834/1‑100 (MQ=255)
            cATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCgg                                                                                        <  2:3554563/100‑1 (MQ=255)
             aTCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGGTAAAGCCGGAGGAAACTTGCCTCCGGc                                                                                       >  1:3046589/1‑100 (MQ=255)
             aTCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCGGc                                                                                       <  1:2143113/100‑1 (MQ=255)
             aTCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGc                                                                                       <  2:1423136/100‑1 (MQ=255)
             aTCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGc                                                                                       <  1:755909/100‑1 (MQ=255)
              tCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCGGct                                                                                      <  1:1045733/100‑2 (MQ=255)
              tCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAACTTGCCTCCGGCa                                                                                      >  1:4076216/1‑100 (MQ=255)
               cTGATTATTTCTGAAATTGTTCACGACAAAGTTGGCAGACCACAACGCATATACCACGGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCa                                                                                      <  2:2588261/99‑1 (MQ=255)
                tGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCAtt                                                                                    >  1:154993/1‑100 (MQ=255)
                 gactattGCTTGAATTGTTCCCGACAAAGGGGGCAGACCACACCTCATAGACCACAGTGGGTGAAATGATAAAGCCGGAGGAAACTTGCCTCCGGCAtt                                                                                    <  2:4128872/96‑1 (MQ=255)
                 gATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTa                                                                                   >  2:3314266/1‑100 (MQ=255)
                  attattGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTAc                                                                                  <  1:1790593/100‑1 (MQ=255)
                   ttattGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCGGCTTTatc                                                                                 <  1:15762/100‑3 (MQ=38)
                       tGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAATGCCGGAGGCAAGTTTCCTCCGGCCTTatcatgt                                                                             <  1:1704483/100‑7 (MQ=25)
                        gCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTa                                                                            <  2:3404280/100‑1 (MQ=255)
                        gCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTa                                                                            >  2:1935758/1‑100 (MQ=255)
                         cTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATACAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTAc                                                                           >  2:3576078/1‑100 (MQ=255)
                         cTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTAc                                                                           <  1:3988215/100‑1 (MQ=255)
                         cTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTAc                                                                           <  2:3147160/100‑1 (MQ=255)
                         cTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTAc                                                                           <  1:2629335/100‑1 (MQ=255)
                         cTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTAc                                                                           <  2:3772057/100‑1 (MQ=255)
                          tGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCGGCTTTatcatgtcac                                                                          <  1:4116606/100‑10 (MQ=25)
                          tGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTATTATTTACt                                                                          >  1:798268/1‑100 (MQ=255)
                          tGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACt                                                                          >  2:3459583/1‑100 (MQ=255)
                            aaaTTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTg                                                                        <  1:2047078/100‑1 (MQ=255)
                             aaTTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGc                                                                       <  1:1159738/100‑1 (MQ=255)
                                  ttCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCtgtg                                                                  <  1:2878157/100‑1 (MQ=255)
                                  ttCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCtgtg                                                                  <  1:1641949/100‑1 (MQ=255)
                                    cACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGAc                                                                >  2:3481237/1‑100 (MQ=255)
                                    cACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGAc                                                                <  1:1964199/100‑1 (MQ=255)
                                        gCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGTCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCgct                                                            >  1:3090773/1‑100 (MQ=255)
                                          aaaGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCgctgc                                                          <  2:1346880/100‑1 (MQ=255)
                                          aaaGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCgctgc                                                          <  2:1279656/100‑1 (MQ=255)
                                           aaGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCgctgct                                                         <  1:3381833/100‑1 (MQ=255)
                                           aaGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCgctgct                                                         <  1:3433674/100‑1 (MQ=255)
                                             gTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGctgctg                                                        <  2:252769/99‑1 (MQ=255)
                                              ttGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGcc                                                      <  1:229519/100‑1 (MQ=255)
                                               tGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCg                                                     >  2:2008883/1‑100 (MQ=255)
                                                ggCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGt                                                    >  2:3819649/1‑100 (MQ=255)
                                                  cAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTc                                                  >  1:521670/1‑100 (MQ=255)
                                                   aGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTcc                                                 <  1:1841938/100‑1 (MQ=255)
                                                    gACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCt                                                >  1:3871236/1‑100 (MQ=255)
                                                     aCCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCt                                                <  2:2798389/99‑1 (MQ=255)
                                                     aCCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTg                                               <  1:324714/100‑1 (MQ=255)
                                                     aCCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTg                                               <  2:3912596/100‑1 (MQ=255)
                                                        acaACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGTCGTTCCTGGCa                                            <  1:1376410/100‑1 (MQ=255)
                                                         caACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCAc                                           >  2:3544446/1‑100 (MQ=255)
                                                         caACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCAc                                           <  1:695788/100‑1 (MQ=255)
                                                          agcGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCAcc                                          <  2:3654785/98‑1 (MQ=255)
                                                               ctataCCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAg                                     <  2:1142641/99‑1 (MQ=255)
                                                                tataCCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGc                                    <  2:1467121/100‑1 (MQ=255)
                                                                 ataCCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAgag                                   >  1:1142648/1‑98 (MQ=255)
                                                                  taCCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGctc                                  <  1:249437/100‑1 (MQ=255)
                                                                    ccACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTg                                >  2:3746986/1‑100 (MQ=255)
                                                                    ccACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTg                                <  2:3037198/100‑1 (MQ=255)
                                                                      acaGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTa                              <  2:3221697/100‑1 (MQ=255)
                                                                       aaGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAg                             <  2:792545/99‑1 (MQ=255)
                                                                        aGTGGGTGACATGATAAAGCCGGAGGAACCTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGc                            <  2:2045305/100‑1 (MQ=255)
                                                                         gTGGCTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCa                           <  1:1219525/100‑1 (MQ=255)
                                                                           gggTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAAc                         >  2:2264727/1‑100 (MQ=255)
                                                                               gACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATa                     <  2:3440558/100‑1 (MQ=255)
                                                                                 cATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAAc                   >  1:3191274/1‑100 (MQ=255)
                                                                                 cATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGAACGCTGCTGCCGTTCCTGGCACCATCAGCGCTGTAGCAACGATAAc                   >  2:1513394/1‑100 (MQ=255)
                                                                                  aTGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACa                  <  1:2239202/100‑1 (MQ=255)
                                                                                   ttaTAATGCCGGAGGTAACTTGCCTCCGGCATTACTATTTACTTGCCTTTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACa                  <  2:308413/97‑1 (MQ=255)
                                                                                   tGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACaa                 >  2:3620462/1‑100 (MQ=255)
                                                                                     gtAATGCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATc               >  1:2081774/2‑100 (MQ=38)
                                                                                       aatgCCGGAGGCAAGTTTCCTCCGGCTTTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCCGCTCTGTAGCAACGATAACAATCac             >  1:3947737/4‑100 (MQ=37)
                                                                                       aatgCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCac             >  2:3833476/4‑100 (MQ=38)
                                                                                       aatgCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCac             >  1:3311069/4‑100 (MQ=38)
                                                                                       aaaGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCac             >  1:1417421/1‑100 (MQ=255)
                                                                                       aaaGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCa              <  2:509765/99‑1 (MQ=255)
                                                                                        ntgCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCacc            >  1:1297324/3‑100 (MQ=255)
                                                                                        atgCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCacc            >  2:2916467/3‑100 (MQ=39)
                                                                                        atgCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCacc            >  1:1328959/3‑100 (MQ=255)
                                                                                        aaGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCacc            <  1:2943440/100‑1 (MQ=255)
                                                                                        aaGCCGGAGGAAACTTGCCTCCGACATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTACTGGAACAATAAGCTCTGTAGCAACGATAACAATCac             >  1:1194609/1‑99 (MQ=255)
                                                                                         tgCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCacca           >  1:3777320/2‑100 (MQ=255)
                                                                                          gCCGGAGGCACGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCACTGTAGCAACGATAACAATCacca           >  2:1208454/1‑99 (MQ=255)
                                                                                          gCCGGAGGCAAGTTTGCTCCGGCATTACTATTTACTTGCCTGGGACCGCTGCTGCCATTCCTAGACCCACAAGCTATTTAGCAAAAACAACACTCgccaa          >  2:740796/1‑100 (MQ=255)
                                                                                          gCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaa          >  1:3534226/1‑100 (MQ=255)
                                                                                          gCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaa          >  1:1746979/1‑100 (MQ=255)
                                                                                          gCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaa          >  1:45672/1‑100 (MQ=255)
                                                                                          gCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACACTCacca           >  2:2351028/1‑99 (MQ=255)
                                                                                          gCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaa          >  1:1311466/1‑100 (MQ=255)
                                                                                          gCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaa          <  1:1680043/100‑1 (MQ=255)
                                                                                          gCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCacc            <  2:2776555/98‑1 (MQ=255)
                                                                                           ccgcaagaAACTTGCCTCCGCCATTACTATTTACCTGCCTGTGACCGCTGCTGCCGTTCCTTGCACCATCAGCTCTGTAGCAACGATAACAATCaccaaa         <  2:4071559/94‑1 (MQ=255)
                                                                                           ccGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaaa         <  2:4049559/100‑1 (MQ=255)
                                                                                                  aaaCTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACcaacaa  >  1:294035/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
TGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAA  >  NC_000913/654462‑654659

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: