Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,647,876 A→G 16.3% intergenic (‑26/‑58) 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,8760AG16.3% 188.3 / 30.3 86intergenic (‑26/‑58)dicC/dicAQin prophage; DNA‑binding transcriptional regulator for DicB/Qin prophage; putative regulator for DicB
Reads supporting (aligned to +/- strand):  ref base A (38/34);  new base G (11/3);  total (49/37)
Fisher's exact test for biased strand distribution p-value = 8.60e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.49e-01

CAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTA  >  NC_000913/1647779‑1647974
                                                                                                 |                                                                                                  
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:144580/100‑1 (MQ=255)
                    gtgCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACa                                                                              <  2:1830077/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                                                                     <  2:1507734/100‑1 (MQ=255)
                             ttgttttTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAg                                                                     <  1:1413120/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:2472432/1‑100 (MQ=255)
                                                                 tAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACaa                                 >  1:1153480/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                         <  1:973979/100‑1 (MQ=255)
                                                                         tAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTa                         >  2:1890904/1‑100 (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                 >  1:541003/6‑100 (MQ=255)
                                                                                 ctatcttGTTAGTTATGACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCgg                 >  2:894581/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          >  1:1663017/1‑100 (MQ=255)
                                                                                        gTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGc          >  1:1822300/1‑100 (MQ=255)
                                                                                        gTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGc          >  2:659707/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:1962114/100‑1 (MQ=255)
                                                                                              cATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAgg    >  1:1407333/1‑100 (MQ=255)
                                                                                                taactaacAAGATAGATGTTAACAAATACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTa  >  1:1561591/1‑100 (MQ=255)
                                                                                                taactaacAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTa  <  2:388332/100‑1 (MQ=255)
                                                                                                taactaacAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTa  <  2:91258/100‑1 (MQ=255)
                                                                                                taactaacAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTa  >  2:1264749/1‑100 (MQ=255)
                                                                                                 |                                                                                                  
CAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAACACCTTTATTGTTAGTCATAACTAACAAGATAGATGTTAACAAAAACATAGTCAATACGATTTAGCATTAGCTAACTATGGAAACAAAAAATTTAACTATCGGCGAACGCATCAGGTA  >  NC_000913/1647779‑1647974

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: