Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,232,481 C→T 19.2% intergenic (+27/+19) umuC → / ← dsbB translesion error‑prone DNA polymerase V subunit; DNA polymerase activity/oxidoreductase that catalyzes reoxidation of DsbA protein disulfide isomerase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,232,4810CT19.2% 191.5 / 39.7 94intergenic (+27/+19)umuC/dsbBtranslesion error‑prone DNA polymerase V subunit; DNA polymerase activity/oxidoreductase that catalyzes reoxidation of DsbA protein disulfide isomerase I
Reads supporting (aligned to +/- strand):  ref base C (45/30);  new base T (15/3);  total (61/33)
Fisher's exact test for biased strand distribution p-value = 9.83e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.33e-02

GCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGA  >  NC_000913/1232385‑1232579
                                                                                                |                                                                                                  
gCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCa                                                                                                 <  1:203815/100‑1 (MQ=255)
 cAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGAAACAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAg                                                                                                >  1:269454/1‑100 (MQ=255)
   gATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCatg                                                                                               <  1:2051037/99‑3 (MQ=255)
    aTGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAg                                                                                             <  1:272214/100‑1 (MQ=255)
       aaGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAgcg                                                                                           <  1:1046008/99‑1 (MQ=255)
        aGCGAGCCATGCTTTCACCACGTTATACAACGCGAATTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGcc                                                                                         <  1:856024/100‑1 (MQ=255)
         gCGAGCCATGCTTTCACCACGTTATACGACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCg                                                                                        <  1:1718419/100‑1 (MQ=255)
             gCCATGCTTTCACCACGTTATACAACGCGACGTTCTGATTTACTGAGGGTCAAATAAATATAACGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATg                                                                                    >  1:990641/1‑100 (MQ=255)
              ccATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgcttttttc                                                                                   <  1:1215345/100‑8 (MQ=255)
              ccATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATgg                                                                                   <  2:589861/100‑1 (MQ=255)
               cATGCTTTCACCACGTTATACAACGAGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAACAGAGCTCCCGCAAGAGCGCCCAATgg                                                                                   >  2:702933/1‑99 (MQ=255)
                aTGCTTTCACGACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGCAATAAAGCGCTCCTGCGGGAgcgcttttttcct                                                                                 <  2:1995530/100‑10 (MQ=37)
                 ttcTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGAtt                                                                                <  2:2100407/98‑1 (MQ=255)
                 ttcTTTCACCAAGTTATACAACGCGAAGTTTTGTTTTACTTTGGGTCAAATAAATATAGAGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGAt                                                                                 <  1:912294/97‑1 (MQ=255)
                  gCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTa                                                                               <  2:229270/100‑1 (MQ=255)
                     ttCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCg                                                                            <  1:119637/100‑1 (MQ=255)
                      tCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGa                                                                           >  1:1874660/1‑100 (MQ=255)
                       caccacGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGAAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGa                                                                           >  2:562998/1‑99 (MQ=255)
                        accacGTTATACAACGCGAAGGTCTGATTTAGTGGGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCTCGCAGGAGCGACGAATGGATTAGCGAc                                                                          >  2:1349269/1‑99 (MQ=255)
                            ggTTATACAACGCGAAGTTCTGATATACTGAGGGAGAAATAAATATAGCGGCAGGAAAAAAGGGCTCCAGCAGGAGAGACGAATGGATTAGCGACCGa                                                                       >  2:2201374/2‑98 (MQ=255)
                            cGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAAc                                                                     <  2:621037/100‑1 (MQ=255)
                             gTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACa                                                                    >  2:2669594/1‑100 (MQ=255)
                             gTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACa                                                                    <  2:1006929/100‑1 (MQ=255)
                                     aCGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATATCGCCAGGGAAAAATTGCTCCCGGAGGAGCTGCGAATGGATTAGCGACCGAACAGATCACGt                                                            >  2:1649692/1‑100 (MQ=255)
                                     aCGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGt                                                            >  2:1427077/1‑100 (MQ=255)
                                        cgAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGtttt                                                         <  1:2720734/100‑1 (MQ=255)
                                           aGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTc                                                      >  2:676215/1‑100 (MQ=255)
                                           aGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTc                                                      >  2:1554156/1‑100 (MQ=255)
                                           aGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGATCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTc                                                      >  2:675598/1‑100 (MQ=255)
                                            gTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCg                                                     <  1:1461224/100‑1 (MQ=255)
                                            gTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGATTTTTCg                                                     >  2:1901967/1‑100 (MQ=255)
                                             ttCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGTCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGc                                                    <  2:1874638/100‑1 (MQ=255)
                                             ttCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGc                                                    >  1:2625413/1‑100 (MQ=255)
                                             ttCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGc                                                    >  1:398033/1‑100 (MQ=255)
                                               cTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCtt                                                  <  2:28666/100‑1 (MQ=255)
                                               cTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATAGAATAGCGACCGAACAGATCACGTTTTTTCGCtt                                                  >  1:2832201/1‑100 (MQ=255)
                                                 gATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAACAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTa                                                >  2:1620871/1‑100 (MQ=255)
                                                 gATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAATCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACAGAACAGATCACGTTTTTTCGCTTTa                                                >  1:1151471/1‑100 (MQ=255)
                                                     tACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACg                                            >  2:487427/1‑100 (MQ=255)
                                                       cTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATGACGATTTTTCGCTTTAAACGGc                                          >  2:1195324/1‑100 (MQ=255)
                                                        tgcgggTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACCGAACACGTTTTTTCGCttt                                                 >  2:2283947/4‑92 (MQ=255)
                                                         gAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCt                                         >  2:1879649/1‑99 (MQ=255)
                                                         gAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTg                                        <  2:9772/100‑1 (MQ=255)
                                                         gAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCACAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCt                                         >  2:1346651/1‑99 (MQ=255)
                                                          aGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTgg                                       <  1:487430/100‑1 (MQ=255)
                                                            ggTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGAGCGAACAAATAACGTTTTTTCGCTTTAAATGGCttgaa                                     >  1:1206981/1‑96 (MQ=255)
                                                             gTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTTGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGaa                                    >  2:1935174/1‑100 (MQ=255)
                                                              tCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGaaa                                   <  2:347163/100‑1 (MQ=255)
                                                              tCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGaaa                                   <  2:100798/100‑1 (MQ=255)
                                                               cAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAAt                                  >  1:645458/1‑100 (MQ=255)
                                                                aaataaatATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATc                                 <  1:1411478/100‑1 (MQ=255)
                                                                aaataaatATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATc                                 <  1:3022296/100‑1 (MQ=255)
                                                                 aataaataTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATca                                >  1:2125128/1‑100 (MQ=255)
                                                                  ataaataTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATcac                               <  1:2918409/100‑1 (MQ=255)
                                                                   taaataTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATcacc                              >  1:1124229/1‑100 (MQ=255)
                                                                     aataTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATcaccac                            <  2:2544092/100‑1 (MQ=255)
                                                                      atatAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCaccacc                           >  1:424077/1‑100 (MQ=255)
                                                                      atatAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGAATTTTCGCTTTAAACGGCTGGGAAATcaccac                            >  2:1439108/1‑99 (MQ=255)
                                                                       tataGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCAccacca                          >  2:2862329/1‑100 (MQ=255)
                                                                           gCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAc                      >  1:2064800/1‑100 (MQ=255)
                                                                             ggCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTg                    >  1:733872/1‑100 (MQ=255)
                                                                             ggCAGGAAAAAAGCGCTCCCACAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTGTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTg                    <  1:643915/100‑1 (MQ=255)
                                                                              gCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGc                   >  1:1704111/1‑100 (MQ=255)
                                                                               cAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCg                  >  2:1264325/1‑100 (MQ=255)
                                                                                 tccattcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCTCTTTATACTGCTGGGAAATCACCACCAACACTGCGAc                >  2:1333540/9‑100 (MQ=37)
                                                                                 ggAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGAc                <  1:2266184/100‑1 (MQ=255)
                                                                                 ggAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGAc                >  2:1114100/1‑100 (MQ=255)
                                                                                 ggAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGAc                >  1:1206522/1‑100 (MQ=255)
                                                                                  gAAAAAAGCGCTCCCGCAGGAGCGCCGACTGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACa               >  1:445889/1‑100 (MQ=255)
                                                                                  gAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACa               >  1:1357706/1‑100 (MQ=255)
                                                                                   gatctggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACaa              <  1:1857973/94‑1 (MQ=255)
                                                                                   cattcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACaa              >  1:1326521/7‑100 (MQ=255)
                                                                                   cattcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACaa              >  1:2074289/7‑100 (MQ=255)
                                                                                   cattcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGACTGGGAAATAACCACCAGAACTGCGACaa              >  1:2603269/7‑100 (MQ=255)
                                                                                   aaaaaaGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACaa              >  2:37123/1‑100 (MQ=255)
                                                                                    attcggcgcTCCTGCGGGAGCGCTGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAAt             >  1:2744964/6‑100 (MQ=255)
                                                                                    attcggcgcTCCTGCGGGAGCGCCGAATTGATTAGCGACCGAACAGATCACGTTTTTTCGCAATAAACAACAAGAAAAACACCACCAGCACTGCGACaa              >  1:1058340/6‑99 (MQ=255)
                                                                                    attcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAATTGCTGGGAAATCACCACCAGCACTGCGaga               >  1:661885/6‑96 (MQ=255)
                                                                                    attcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAAt             >  1:2334292/6‑100 (MQ=255)
                                                                                    attcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAAt             >  1:2439125/6‑100 (MQ=255)
                                                                                    attcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAAACACCACAAGCACTGCGACAAt             >  1:1741839/6‑100 (MQ=255)
                                                                                    attcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCCCTTTAAACGGCTGGGAAATCACCACCAACACTGCGACaa              >  1:2229220/6‑99 (MQ=255)
                                                                                       cggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGTTTTAAACTGCTGGGAAATCACAACCAGCACTGCGACAATCa           >  2:2361605/3‑99 (MQ=255)
                                                                                       aaGCGCTCCCGCAGGAGCGCTGAATGGATTAGTGTCCGAACAGGTCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACCGCGACAATc            <  1:1194902/98‑1 (MQ=255)
                                                                                        ggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACaaagatc         >  2:1857952/2‑95 (MQ=255)
                                                                                        ggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAgg         >  1:946715/2‑100 (MQ=255)
                                                                                        ggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAgg         >  2:1984179/2‑100 (MQ=255)
                                                                                         gcgcTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTCCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGt        >  2:1330423/1‑100 (MQ=255)
                                                                                          cgcTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTa       <  1:1099269/100‑1 (MQ=255)
                                                                                          cgcTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTa       >  1:1811035/1‑100 (MQ=255)
                                                                                          cgcTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACAACCAGCACTGCGACAATCAGGTa       >  2:575924/1‑100 (MQ=255)
                                                                                          cgcTCCCGCAGGAGCACCGAATAGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTAGGAAAACACCACCAACACTGCGACAATCAGGTa       >  2:83968/1‑100 (MQ=255)
                                                                                           gcTCCCGCAGTAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTCTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAgg         <  2:315152/97‑1 (MQ=255)
                                                                                           gcTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTaa      >  1:2222724/1‑100 (MQ=255)
                                                                                              cccGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCg   >  1:232493/1‑100 (MQ=255)
                                                                                              cccGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCg   >  2:832301/1‑100 (MQ=255)
                                                                                              cccGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGATTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGAGACAATCAGGTaa      >  2:1772819/1‑97 (MQ=255)
                                                                                               ccGCAGGAGCGCCGAATGGATTAGCGACCGAAGAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGa  >  1:2885974/1‑100 (MQ=255)
                                                                                                |                                                                                                  
GCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGA  >  NC_000913/1232385‑1232579

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