Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 654,563 C→G 15.5% intergenic (+21/+20) 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,5630CG15.5% 553.3 / 86.3 200intergenic (+21/+20)citB/dcuCresponse regulator in two‑component regulatory system with CitA/anaerobic C4‑dicarboxylate transport
Reads supporting (aligned to +/- strand):  ref base C (71/98);  new base G (17/14);  total (88/112)
Fisher's exact test for biased strand distribution p-value = 2.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.51e-01

Only 100 of 203 total aligned reads displayed.

CGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGC  >  NC_000913/654464‑654661
                                                                                                   |                                                                                                  
cgcCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAAc                                                                                                    <  1:3226097/100‑1 (MQ=255)
 gcCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACt                                                                                                   <  1:1248611/100‑1 (MQ=255)
 gcCAGCCGCCATCTGATTATTGCTGAAATTGTTCACAGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACt                                                                                                   >  2:376553/1‑100 (MQ=255)
     gccgccATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGcc                                                                                               >  1:3361563/1‑100 (MQ=255)
     gccgccATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGcc                                                                                               <  2:604288/100‑1 (MQ=255)
     gccgccATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGcc                                                                                               >  1:128086/1‑100 (MQ=255)
     gccgccATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGcc                                                                                               >  2:3539187/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)
           aTCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCGGc                                                                                         <  2:1264376/100‑1 (MQ=255)
           aTCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGc                                                                                         <  1:755909/100‑1 (MQ=255)
            tCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCGGct                                                                                        <  1:1045733/100‑2 (MQ=255)
            tCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCGGct                                                                                        <  1:3512537/100‑2 (MQ=255)
             cTGATTATTTCTGAAATTGTTCACGACAAAGTTGGCAGACCACAACGCATATACCACGGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCa                                                                                        <  2:2588261/99‑1 (MQ=255)
             cTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACCTGATAACGCCGGAGGAAACTTGCCTCCgg                                                                                          >  2:1508488/1‑97 (MQ=255)
              tGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCAtt                                                                                      >  1:154993/1‑100 (MQ=255)
              tGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCAtt                                                                                      >  2:1159731/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)
                      gCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTAttt                                                                               <  1:3398260/99‑1 (MQ=255)
                      gCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTa                                                                              <  2:3404280/100‑1 (MQ=255)
                       cTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATACAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTAc                                                                             >  2:3576078/1‑100 (MQ=255)
                       cTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTAc                                                                             <  2:3147160/100‑1 (MQ=255)
                       cTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTAc                                                                             <  2:3772057/100‑1 (MQ=255)
                        tGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTATTATTTACt                                                                            >  1:798268/1‑100 (MQ=255)
                        tGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACt                                                                            >  2:3459583/1‑100 (MQ=255)
                         gAAATTGTTCACGGCAAAGTTGGCAGACCACAAAGCATATACCACAGCCGGTGGCATTATAACGCCGGAGGAAACTTGCCTCCGGCAGTACTATTTAc                                                                             >  2:213457/1‑98 (MQ=255)
                          aaaTTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTg                                                                          <  1:2047078/100‑1 (MQ=255)
                           aaTTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGc                                                                         <  1:1159738/100‑1 (MQ=255)
                               gTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCtgt                                                                     >  1:2887035/1‑100 (MQ=255)
                                ttCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCtgtg                                                                    <  1:3955009/100‑1 (MQ=255)
                                ttCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCtgtg                                                                    <  1:1641949/100‑1 (MQ=255)
                                 tCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGa                                                                   >  1:3029751/1‑100 (MQ=255)
                                  cACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGAc                                                                  >  1:4095151/1‑100 (MQ=255)
                                  cACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGAc                                                                  <  1:2720180/100‑1 (MQ=255)
                                  cACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGAc                                                                  <  1:1964199/100‑1 (MQ=255)
                                  cACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGAc                                                                  >  2:3481237/1‑100 (MQ=255)
                                    cGGCAAAGTTGGCAGACCACAACGCATATTCCACAGTTGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGAc                                                                  <  2:1442961/98‑1 (MQ=255)
                                           gTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGctgctg                                                          <  2:252769/99‑1 (MQ=255)
                                            ttGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGcc                                                        <  1:229519/100‑1 (MQ=255)
                                             tGGCAGACCACCACGCATATACCACTGTGGCTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCg                                                       <  2:2350567/100‑1 (MQ=255)
                                              ggCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGt                                                      >  2:3819649/1‑100 (MQ=255)
                                               gCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGtt                                                     <  1:382791/100‑1 (MQ=255)
                                               gCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGtt                                                     >  1:3642688/1‑100 (MQ=255)
                                                 aGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTcc                                                   <  1:30958/100‑1 (MQ=255)
                                                 aGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTcc                                                   <  1:1841938/100‑1 (MQ=255)
                                                   aCCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCt                                                  <  2:2798389/99‑1 (MQ=255)
                                                   aCCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTg                                                 <  2:3912596/100‑1 (MQ=255)
                                                     cacaACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGc                                               <  1:2940608/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)
                                                              tataCCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGc                                      <  2:1943505/100‑1 (MQ=255)
                                                               ataCCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAgag                                     >  1:1142648/1‑98 (MQ=255)
                                                               ataCCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGct                                     <  1:4000612/100‑1 (MQ=255)
                                                                taCCACAGTGGGTGACATGATAACGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGctc                                    <  1:2662012/100‑1 (MQ=255)
                                                                taCCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGctc                                    <  1:249437/100‑1 (MQ=255)
                                                                taCCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGctc                                    <  1:3910734/100‑1 (MQ=255)
                                                                  ccACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGctcag                                  >  2:2501307/1‑98 (MQ=255)
                                                                  ccACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTg                                  >  2:3746986/1‑100 (MQ=255)
                                                                  ccACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTg                                  >  1:532086/1‑100 (MQ=255)
                                                                  ccACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTg                                  >  1:517187/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)
                                                                       gTGGCTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCa                             <  1:1219525/100‑1 (MQ=255)
                                                                            tGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGAt                        >  1:3263020/1‑100 (MQ=255)
                                                                             gACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATa                       <  2:3440558/100‑1 (MQ=255)
                                                                              aaatagtAATGCCGTAGGCAAGTTTCCTCCGGCTTTACTATGTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCACGTCTGTAGCAACGATaa                      >  2:2408597/7‑100 (MQ=21)
                                                                                aTGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACa                    <  2:1617333/100‑1 (MQ=255)
                                                                                aTGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACa                    <  1:2239202/100‑1 (MQ=255)
                                                                                 ttaTAATGCCGGAGGTAACTTGCCTCCGGCATTACTATTTACTTGCCTTTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACa                    <  2:308413/97‑1 (MQ=255)
                                                                                 tGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACaa                   >  2:3620462/1‑100 (MQ=255)
                                                                                   gtAATGCCGGAGGCAAGTTTCCTCCGGCTTTACTATTTACTTGCCTGTGACCGCTGCTGCCCTTCCTGGCACCATCAGCTATGTAGCCACGACAACAATc                 >  2:1338214/2‑100 (MQ=25)
                                                                                   gtAATGCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATc                 >  1:2081774/2‑100 (MQ=38)
                                                                                    tAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCa                >  2:2026696/1‑100 (MQ=255)
                                                                                    tAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCa                >  1:2569238/1‑100 (MQ=255)
                                                                                     aatgCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCac               >  2:3833476/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)
                                                                                      aaGCCGGAGGAAACTTGCCTCCGACATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTACTGGAACAATAAGCTCTGTAGCAACGATAACAATCac               >  1:1194609/1‑99 (MQ=255)
                                                                                        gCCGGAGGCAAGTTTGCTCCGGCATTACTATTTACTTGCCTGGGACCGCTGCTGCCATTCCTAGACCCACAAGCTATTTAGCAAAAACAACACTCgccaa            >  2:740796/1‑100 (MQ=255)
                                                                                        gCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaa            >  1:1746979/1‑100 (MQ=255)
                                                                                        gCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaa            >  1:45672/1‑100 (MQ=255)
                                                                                        gCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaa            <  1:1680043/100‑1 (MQ=255)
                                                                                        gCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaa            >  1:1311466/1‑100 (MQ=255)
                                                                                        gCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCacc              <  2:2776555/98‑1 (MQ=255)
                                                                                         ccgcaagaAACTTGCCTCCGCCATTACTATTTACCTGCCTGTGACCGCTGCTGCCGTTCCTTGCACCATCAGCTCTGTAGCAACGATAACAATCaccaaa           <  2:4071559/94‑1 (MQ=255)
                                                                                         ccGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaaa           <  2:4049559/100‑1 (MQ=255)
                                                                                           ggaggaAACTTGCCTCCGGCATTACTATTGACTTGCCTGTGACCGCTGCTGCCGATCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaaac          <  2:1354173/99‑1 (MQ=255)
                                                                                                 aaCTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAg   >  2:1124295/1‑100 (MQ=255)
                                                                                                  aCTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGc  <  1:442963/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
CGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGC  >  NC_000913/654464‑654661

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: