Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 654,563 C→G 15.6% 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.6% 546.2 / 81.7 197intergenic (+21/+20)citB/dcuCresponse regulator in two‑component regulatory system with CitA/anaerobic C4‑dicarboxylate transport
Reads supporting (aligned to +/- strand):  ref base C (79/87);  new base G (15/16);  total (94/103)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 201 total aligned reads displayed.

CGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCA  >  NC_000913/654464‑654662
                                                                                                   |                                                                                                   
cgcCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAAc                                                                                                     <  2:4147063/100‑1 (MQ=255)
 gcCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACt                                                                                                    >  1:887707/1‑100 (MQ=255)
  ccAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGGATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAAc                                                                                                     >  1:1652792/1‑98 (MQ=255)
        gccATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTtcctcc                                                                                             <  1:1780480/100‑1 (MQ=255)
        gccATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTcc                                                                                             <  2:2106443/100‑1 (MQ=255)
           aTCTGATTATTGCTGAAATTGTTTACGGTAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCGGc                                                                                          <  1:1699587/100‑1 (MQ=255)
           aTCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCGGc                                                                                          <  1:3362134/100‑1 (MQ=255)
           aTCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCGGc                                                                                          <  1:1715313/100‑1 (MQ=255)
            tCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACTACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCGGct                                                                                         <  1:3440446/100‑2 (MQ=255)
            tCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCGGct                                                                                         <  1:2639623/100‑2 (MQ=255)
            tCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCGGct                                                                                         <  2:1146194/100‑2 (MQ=255)
             cTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCAt                                                                                        >  2:1061320/1‑100 (MQ=255)
              tGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCGGctt                                                                                        <  2:3981806/99‑3 (MQ=39)
               gATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTa                                                                                      >  2:910465/1‑100 (MQ=255)
                attattGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTAc                                                                                     <  1:358054/100‑1 (MQ=255)
                  tattGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTa                                                                                   >  1:2541015/1‑100 (MQ=255)
                    ttGCTGAAATTGTTCACGGCAAAGTTGGAAGACCACAACGAATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTAtt                                                                                 >  1:2960274/1‑100 (MQ=255)
                     tGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTAttt                                                                                >  2:4152696/1‑100 (MQ=255)
                     tGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTAttt                                                                                >  1:2238136/1‑100 (MQ=255)
                     tGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTAttt                                                                                <  1:2829555/100‑1 (MQ=255)
                      gCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTa                                                                               >  1:708220/1‑100 (MQ=255)
                       cTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTAc                                                                              >  2:1546453/1‑100 (MQ=255)
                       cTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTAc                                                                              >  2:3288295/1‑100 (MQ=255)
                        tGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACt                                                                             >  1:781395/1‑100 (MQ=255)
                        tGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACt                                                                             >  1:1019286/1‑100 (MQ=255)
                          aaaTTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTg                                                                           <  1:3939845/100‑1 (MQ=255)
                            aTTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGcc                                                                         <  2:2482644/100‑1 (MQ=255)
                            aTTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGcc                                                                         <  1:3226157/100‑1 (MQ=255)
                              tGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCtg                                                                       <  2:1217075/100‑1 (MQ=255)
                               gTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCtgt                                                                      <  2:28674/100‑1 (MQ=255)
                               gTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCtgt                                                                      >  2:2237181/1‑100 (MQ=255)
                               gTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCtgt                                                                      >  2:373726/1‑100 (MQ=255)
                                 tCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGa                                                                    >  1:1904793/1‑100 (MQ=255)
                                 tCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGa                                                                    <  2:3088750/100‑1 (MQ=255)
                                  cACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTACTTACTTGCCTGTGAc                                                                   >  1:2925885/1‑100 (MQ=255)
                                    cGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCg                                                                 >  1:3517803/1‑100 (MQ=255)
                                       cAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCgctg                                                              >  1:1122656/1‑100 (MQ=255)
                                         aaGTTGGCAGTCCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCgctgct                                                            <  1:3736484/100‑1 (MQ=255)
                                         aaGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCAGTACTATTTACTTGCCTGTGACCgctgct                                                            <  1:537205/100‑1 (MQ=255)
                                           gTTGGCAGACCACCACGCATCTACCCCAGTGGGTGACATGTTCAAGCCGGAGGAATCTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCgctgctc                                                           <  1:2752312/99‑2 (MQ=255)
                                           gTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCtgctgc                                                          <  2:781396/100‑1 (MQ=255)
                                           gTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCtgctgc                                                          <  2:1725650/100‑1 (MQ=255)
                                             tGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCg                                                        >  2:465638/1‑100 (MQ=255)
                                             tGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCg                                                        <  2:3229536/100‑1 (MQ=255)
                                             tGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCg                                                        <  1:4163083/100‑1 (MQ=255)
                                              ggCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGt                                                       >  1:3059839/1‑100 (MQ=255)
                                               gCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGtt                                                      >  2:2704808/1‑100 (MQ=255)
                                               gCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGtt                                                      >  2:2246989/1‑100 (MQ=255)
                                               gCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGtt                                                      >  2:597961/1‑100 (MQ=255)
                                                  gACCACAACGCATATACGACAGTGAATGACATCATAAAGCCGGGGGAAACTCGCCTCCGGCCTTACTGATTCCTTGCCTGTGCACgctgc                                                             >  1:2623991/1‑90 (MQ=255)
                                                  gACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCt                                                   <  2:1814608/100‑1 (MQ=255)
                                                  gACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCt                                                   <  2:2266614/100‑1 (MQ=255)
                                                   aCCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTg                                                  <  1:1858473/100‑1 (MQ=255)
                                                      acaACGCATATACCACAGTGGTTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCa                                               <  1:2175105/100‑1 (MQ=255)
                                                      acaACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCa                                               <  1:247739/100‑1 (MQ=255)
                                                       caACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCAc                                              <  2:2157573/100‑1 (MQ=255)
                                                       caACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCAc                                              >  1:4006550/1‑100 (MQ=255)
                                                        aaCGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGGCTGTGACCGCTGCTGCCGTTCCTGGCAc                                              >  1:4325702/1‑99 (MQ=255)
                                                         aCGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCa                                            >  2:3650376/1‑100 (MQ=255)
                                                           gCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATc                                          <  1:2491881/100‑1 (MQ=255)
                                                            cATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCa                                         >  2:2432087/1‑100 (MQ=255)
                                                            cATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCa                                         >  1:4134397/1‑100 (MQ=255)
                                                            cATATACCACAGTGGACGACATGATAAAGCCGGAAGAAACTTGCCTCCGGCATTAATATTTACTTGCCTGTGACCGATGCTGCCGTTCCTGGCACCATCa                                         >  1:418125/1‑100 (MQ=255)
                                                                taCCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGctc                                     <  1:3259812/100‑1 (MQ=255)
                                                                  ccACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTg                                   <  1:83389/100‑1 (MQ=255)
                                                                   cacaGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGt                                  >  2:3274321/1‑100 (MQ=255)
                                                                   cacaGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGt                                  >  2:3690466/1‑100 (MQ=255)
                                                                    acaGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGt                                  <  1:1354165/99‑1 (MQ=255)
                                                                      aGTGGGTGACCTGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGc                               <  1:2899478/100‑1 (MQ=255)
                                                                       gTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACAACTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCa                              >  2:556783/1‑100 (MQ=255)
                                                                         gggTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAAc                            >  2:1455649/1‑100 (MQ=255)
                                                                          ggTGACAAGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGCCCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGGAAc                            >  2:2315331/1‑99 (MQ=255)
                                                                             gACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATa                        <  2:1486466/100‑1 (MQ=255)
                                                                               cATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAAc                      >  1:2685270/1‑100 (MQ=255)
                                                                               cATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAAc                      >  1:61511/1‑100 (MQ=255)
                                                                                   gtAATGCCGGAGGCAAGTTTCCTCCGGCTTTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATc                  >  2:3596949/2‑100 (MQ=37)
                                                                                   gtAATGCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATc                  >  1:730069/2‑100 (MQ=38)
                                                                                   gtAATGCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATc                  >  2:342800/2‑100 (MQ=38)
                                                                                   aTAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATc                  <  1:204067/100‑1 (MQ=255)
                                                                                    tAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCa                 >  1:2002383/1‑100 (MQ=255)
                                                                                     aatgCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCCCCATCACCTCTGTAGCAACGATAACAATCac                >  1:2061184/4‑100 (MQ=38)
                                                                                     aatgCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCac                >  1:1944511/4‑100 (MQ=38)
                                                                                        gCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGAGCGCTGCTGCAGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCacca              >  1:3646615/1‑99 (MQ=255)
                                                                                        gCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaa             >  1:2383257/1‑100 (MQ=255)
                                                                                        gCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaa             >  2:2861138/1‑100 (MQ=255)
                                                                                        gCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaa             >  2:1919458/1‑100 (MQ=255)
                                                                                        gCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaa             <  2:2880841/100‑1 (MQ=255)
                                                                                        gCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaa             >  2:3807383/1‑100 (MQ=255)
                                                                                          cGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaaac           <  2:3021966/100‑1 (MQ=255)
                                                                                          cGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaaac           <  1:1455643/100‑1 (MQ=255)
                                                                                          cGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaaac           >  2:43951/1‑100 (MQ=255)
                                                                                          cGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaaa            <  1:3409036/99‑1 (MQ=255)
                                                                                           ggaggaAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaaacc          <  2:2812522/100‑1 (MQ=255)
                                                                                           ggaggaAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaaacc          >  2:3827150/1‑100 (MQ=255)
                                                                                           ggaggaAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaaac           <  1:4070853/99‑1 (MQ=255)
                                                                                           ggaggaAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTCCTGCCGCTCCTTGCACCATCAGCTCTGTAGCAACGATAACAATCaccaa             >  1:1080525/1‑97 (MQ=255)
                                                                                            gaggaAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaaacca         >  2:1602186/1‑100 (MQ=255)
                                                                                              ggaAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACcaac       <  1:3156999/100‑1 (MQ=255)
                                                                                              ggaAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACcaac       <  1:3803387/100‑1 (MQ=255)
                                                                                                  aCTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGc   <  2:2634162/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
CGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCA  >  NC_000913/654464‑654662

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: