Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,647,873 C→T 16.0% intergenic (‑23/‑61) dicC ← / → dicA Qin prophage; DNA‑binding transcriptional regulator for DicB/Qin prophage; putative regulator for DicB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,647,8730CT16.0% 213.2 / 27.7 89intergenic (‑23/‑61)dicC/dicAQin prophage; DNA‑binding transcriptional regulator for DicB/Qin prophage; putative regulator for DicB
Reads supporting (aligned to +/- strand):  ref base C (40/34);  new base T (11/3);  total (51/38)
Fisher's exact test for biased strand distribution p-value = 1.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.98e-01

AAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGG  >  NC_000913/1647775‑1647972
                                                                                                  |                                                                                                   
aaGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTcn                                                                                                    <  2:165351/100‑2 (MQ=255)
 aGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCAt                                                                                                   <  1:135808/100‑1 (MQ=255)
 aGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCAt                                                                                                   <  1:1628421/100‑1 (MQ=255)
    cAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCAtaac                                                                                                >  2:431418/1‑100 (MQ=255)
     aaaCGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGCTAGTCAtaact                                                                                               >  2:2388743/1‑100 (MQ=255)
      aaCGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCAtaacta                                                                                              >  1:717111/1‑100 (MQ=255)
        cGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCAtaactaac                                                                                            >  1:2533148/1‑100 (MQ=255)
          aaTACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACaa                                                                                          <  2:516813/100‑1 (MQ=255)
           aTACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTTATGACTAACaat                                                                                         <  2:1386134/100‑2 (MQ=255)
             aCCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAagat                                                                                       <  2:932291/100‑1 (MQ=255)
             aCCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAagat                                                                                       <  2:654711/100‑1 (MQ=255)
              ccTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAagata                                                                                      >  2:1814571/1‑100 (MQ=255)
                 gctgctTGGGCAAGTTTTGTTTGTGAACCGAAATACAAAAGAGCGTCAGGTTGAAGCATTTGGAACACCTTTAGTGTTAGTCATAACTAACaa                                                                                          <  2:1963088/93‑1 (MQ=255)
                  ctgctTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATg                                                                                  <  1:1997870/100‑1 (MQ=255)
                   tgctTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGt                                                                                 >  1:538350/1‑100 (MQ=255)
                   tgctTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGt                                                                                 >  2:2106892/1‑100 (MQ=255)
                    gctTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGtt                                                                                <  1:1560440/100‑1 (MQ=255)
                     ctTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGCCATAACTAACAAGATAGATGTTa                                                                               <  2:1706720/100‑1 (MQ=255)
                      ttGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTTATGACTAACAATAAAGGTGtttt                                                                              <  1:883312/100‑3 (MQ=37)
                      ttGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTaa                                                                              <  1:528738/100‑1 (MQ=255)
                       tgtgCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAAc                                                                             <  1:1575450/100‑1 (MQ=255)
                       tgtgCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCGTTATTGTTAGTCATAACTAACAAGATAGATGTTAAc                                                                             >  1:1798614/1‑100 (MQ=255)
                        gtgCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACa                                                                            <  2:1830077/100‑1 (MQ=255)
                        gtgCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACa                                                                            <  2:144580/100‑1 (MQ=255)
                           cAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACaaaa                                                                         <  2:1524203/100‑1 (MQ=255)
                            aaGTTTTGTTTTTGAAAAGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACAACTTTATTGTTAGTCATAACTAACAAGATAAATATTAACaaaaa                                                                        <  2:1582392/100‑1 (MQ=255)
                              gttttgttttTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTTATGACTAACAAAATAGATGTTAAAAAAAACa                                                                      >  1:2136408/1‑100 (MQ=255)
                              gttttgttttTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACa                                                                      <  1:548343/100‑1 (MQ=255)
                               ttttgttttTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACAt                                                                     >  1:1450445/1‑100 (MQ=255)
                                 ttgttttTGAACCGAAATACAAAAGAGCGTCATTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAg                                                                   >  2:2341867/1‑100 (MQ=255)
                                 ttgttttTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAg                                                                   <  1:1413120/100‑1 (MQ=255)
                                 ttgttttTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAg                                                                   <  2:1507734/100‑1 (MQ=255)
                                  tgttttTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGt                                                                  <  1:1125964/100‑1 (MQ=255)
                                    tttttGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCa                                                                >  1:2501250/1‑100 (MQ=255)
                                      tttGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAAt                                                              <  1:899283/100‑1 (MQ=255)
                                      tttGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAAt                                                              <  2:1596745/100‑1 (MQ=255)
                                       ttGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATa                                                             >  1:920552/1‑100 (MQ=255)
                                        tGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATAc                                                            >  1:2497697/1‑100 (MQ=255)
                                        tGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATAc                                                            <  1:791864/100‑1 (MQ=255)
                                          aaCCGAAATACAAAAGAGAGTCAGTTTTAAGCAATAAAAACAACTTTATTGTTAGTGCTATCTTAAAATATAAAAAAAAAAAAAAACATAGTCAATACGa                                                          <  2:1459613/100‑1 (MQ=255)
                                           aCCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATATTAACAAAAAGATAGTGCATACTAt                                                         >  2:2462724/1‑100 (MQ=255)
                                            ccGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGAtt                                                        >  1:1872221/1‑100 (MQ=255)
                                            ccGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGAtt                                                        >  2:328009/1‑100 (MQ=255)
                                              gAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTa                                                      >  1:1982928/1‑100 (MQ=255)
                                                aaTACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGc                                                    <  2:2044309/100‑1 (MQ=255)
                                                 aTACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCa                                                   <  1:2271281/100‑1 (MQ=255)
                                                   aCAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCAtt                                                 <  1:1324428/100‑1 (MQ=255)
                                                            cGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTa                                        <  2:1872207/100‑1 (MQ=255)
                                                            cGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTa                                        >  2:1506148/1‑100 (MQ=255)
                                                             gTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTa                                        >  2:2123335/1‑99 (MQ=255)
                                                             gTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTAt                                       >  2:1769501/1‑100 (MQ=255)
                                                              tCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATg                                      >  1:502730/1‑100 (MQ=255)
                                                               cAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATgg                                     <  1:2106906/100‑1 (MQ=255)
                                                                aGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGa                                    <  1:1192059/100‑1 (MQ=255)
                                                                    ttAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACa                                <  1:1375632/100‑1 (MQ=255)
                                                                     tAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACaa                               >  1:1153480/1‑100 (MQ=255)
                                                                     tAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACaa                               >  1:2472432/1‑100 (MQ=255)
                                                                      aaGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGGTAACAAGAACATAGTCAATACGATTTAGCATTAGCTAACTATGGaa                                   <  2:573419/95‑1 (MQ=255)
                                                                       aGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACaaaa                             >  1:1202373/1‑100 (MQ=255)
                                                                       aGCATTTAAAACACCTGTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTAAATACGATTTAGCATTAGCTAACTATGGAAACaaaa                             >  2:755113/1‑100 (MQ=255)
                                                                        gCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACaaaaa                            >  2:8890/1‑100 (MQ=255)
                                                                          aTTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAAt                          >  2:339306/1‑100 (MQ=255)
                                                                             tAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTa                       >  2:1890904/1‑100 (MQ=255)
                                                                             tAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTa                       <  1:973979/100‑1 (MQ=255)
                                                                                  catctatcttGTTAGTTATGACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTAt                  >  1:432864/9‑100 (MQ=38)
                                                                                   atctatcttGTTAGTTATGACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATc                 >  2:821328/8‑100 (MQ=38)
                                                                                    ccTTTATTGTTAGTCATAACTAACAAGATAGCTGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCg                >  2:1385816/1‑100 (MQ=255)
                                                                                    ccTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCg                <  2:1153472/100‑1 (MQ=255)
                                                                                     ctatcttGTTAGTTATGACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCgg               >  2:894581/6‑100 (MQ=255)
                                                                                     ctatcttGTTAGTTATGACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCgg               >  1:541003/6‑100 (MQ=255)
                                                                                     ctatcttGTTAGTTATGACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCgg               >  2:1479175/6‑100 (MQ=255)
                                                                                     cTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCgg               >  1:848538/1‑100 (MQ=255)
                                                                                     cTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCgg               <  1:2386981/100‑1 (MQ=255)
                                                                                      tttATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGc              >  1:2201388/1‑100 (MQ=255)
                                                                                      tattttGTTAGTTATGACTAACAAGATAGATGTTAACAAAAACATTGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGc              >  1:1477066/5‑100 (MQ=255)
                                                                                      tatcttGTTAGTTATGACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGc              >  1:807971/5‑100 (MQ=255)
                                                                                        tcttGTTAGTTATGACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGa            >  2:2025448/3‑100 (MQ=255)
                                                                                        tcttGTTAGTTATGACTAACAAGATAGATGTTAACAAAAACATAGTCAATACAATTTAGCATTAACTAACTATGGAAACAAAAAATTTAACTATCAGcaa            >  2:285052/3‑98 (MQ=255)
                                                                                        tATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCAAACTATAGAAACAAAAAATTTAACTATCGGCGa            >  2:1155512/1‑100 (MQ=255)
                                                                                         cttGTTAGTTATGACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGaa           >  2:640534/2‑100 (MQ=255)
                                                                                         aTTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGaa           <  1:80034/100‑1 (MQ=255)
                                                                                          ttGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAAc          >  2:2119928/1‑100 (MQ=255)
                                                                                           tGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACg         >  1:148354/1‑100 (MQ=255)
                                                                                            gTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGc        >  2:659707/1‑100 (MQ=255)
                                                                                            gTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGc        >  1:1663017/1‑100 (MQ=255)
                                                                                            gTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGc        >  1:1822300/1‑100 (MQ=255)
                                                                                             ttAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCa       >  1:1782408/1‑100 (MQ=255)
                                                                                               aGTTATGACTAACAAGATAGATGTTAACAACAACATAGTCAACACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATc     <  1:640541/100‑1 (MQ=255)
                                                                                                gTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCa    >  1:1817204/1‑100 (MQ=255)
                                                                                                gTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCa    >  1:2408475/1‑100 (MQ=255)
                                                                                                  cATAACTAACAAGATAGATGTTAACAAAAACATATTCAATACGATTTAGCATTAGCTAACTACGGAAACAAAAACTTTAACTATCGGCGAACGCATCa    >  2:1823463/1‑98 (MQ=255)
                                                                                                  cATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAgg  >  1:1407333/1‑100 (MQ=255)
                                                                                                  cATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAgg  <  1:1962114/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
AAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGG  >  NC_000913/1647775‑1647972

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