Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,832,194 A→T 17.6% intergenic (+28/+25) yicH → / ← yicI putative inner membrane‑anchored periplasmic AsmA family protein/putative alpha‑glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,832,1940AT17.6% 438.4 / 82.0 187intergenic (+28/+25)yicH/yicIputative inner membrane‑anchored periplasmic AsmA family protein/putative alpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base A (79/74);  new base T (18/15);  total (98/89)
Fisher's exact test for biased strand distribution p-value = 8.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.59e-01

Only 100 of 190 total aligned reads displayed.

CAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTC  >  NC_000913/3832095‑3832293
                                                                                                   |                                                                                                   
caAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGa                                                                                                     <  2:2411770/100‑1 (MQ=255)
  aGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACg                                                                                                   >  1:1463968/1‑100 (MQ=255)
     gACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAc                                                                                                >  2:4054319/1‑100 (MQ=255)
     gACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAc                                                                                                >  1:1850971/1‑100 (MQ=255)
      aCTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCCACGCGTACGACTCAGAcatga                                                                                                >  2:3968173/1‑95 (MQ=255)
       cTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTAccc                                                                                              <  2:1223309/100‑1 (MQ=255)
       cTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAccc                                                                                              >  2:576199/1‑100 (MQ=255)
        tCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCg                                                                                             <  2:2832301/100‑1 (MQ=255)
         cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACTTCTGAGTCGTACCCgg                                                                                            <  2:2003048/100‑1 (MQ=255)
         cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCgg                                                                                            <  1:2841571/100‑1 (MQ=39)
         cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCgg                                                                                            <  2:1358655/100‑1 (MQ=39)
          aaaCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCGGc                                                                                           <  2:3701425/100‑1 (MQ=39)
           aaCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCGGcg                                                                                          <  2:1600448/100‑2 (MQ=38)
            aCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCGGcga                                                                                         <  2:2951611/100‑3 (MQ=38)
            aCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCAt                                                                                         >  1:2543283/1‑100 (MQ=255)
             ccATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCAtt                                                                                        <  2:656134/100‑1 (MQ=255)
             ccATGAATGCTGGCATCGTGTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGCCTCAGACGTACCCGGCAtt                                                                                        >  1:4196677/1‑100 (MQ=255)
              cATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCGGcgaga                                                                                       <  2:1173591/100‑5 (MQ=38)
              cATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTc                                                                                       <  2:1879562/100‑1 (MQ=255)
               aTGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTcc                                                                                      <  2:732319/100‑1 (MQ=255)
                tGAATGCTGGCATCGTTTCCAGCGGTGAATTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTcc                                                                                      <  2:4365049/99‑1 (MQ=255)
                tGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCa                                                                                     <  2:2287155/100‑1 (MQ=255)
                  naTGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATc                                                                                   >  1:3508811/2‑100 (MQ=255)
                  aaTGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATTTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATc                                                                                   <  1:2354264/100‑1 (MQ=255)
                  aaTGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATc                                                                                   <  1:3923525/100‑1 (MQ=255)
                  aaTGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATc                                                                                   <  2:2627348/100‑1 (MQ=255)
                      cTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCa                                                                                  <  1:2513762/97‑1 (MQ=255)
                         gCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGAt                                                                            >  1:1579822/1‑100 (MQ=255)
                         gCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGAt                                                                            >  2:2107510/1‑100 (MQ=255)
                             cGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGCCGTACCCGGCATTCCATCAATAGATTa                                                                          >  1:1269006/1‑98 (MQ=255)
                              gTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACaa                                                                       >  2:3336011/1‑100 (MQ=255)
                              gTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACaa                                                                       >  1:2084474/1‑100 (MQ=255)
                                   cAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTaa                                                                  <  1:3580780/100‑1 (MQ=255)
                                      cGGTGAAGTTACATCGACGGAGCCGGTGCGGTAATCATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTTCAACGTAATTg                                                               <  1:2314879/100‑1 (MQ=255)
                                      cGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTg                                                               <  2:1040125/100‑1 (MQ=255)
                                       ggTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGt                                                              >  2:534465/1‑100 (MQ=255)
                                          gAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGACGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAg                                                           <  1:440653/100‑1 (MQ=255)
                                           aaGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAgc                                                          <  1:27586/100‑1 (MQ=255)
                                              ttccaccGACGGAGCCTGTGCGGTAAACATCAAGCTCGCCGGGGACGACGCAGACGTACCCGGCATTCCATTAATAGATTACGACGTAATTGTCAGCGCa                                                       <  1:3641986/94‑1 (MQ=255)
                                              ttACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCa                                                       <  2:2614526/100‑1 (MQ=255)
                                               tACAGCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCAGTCCATCAATAGATTACAACGTAATTGTCGGCGCa                                                       <  1:1231067/99‑1 (MQ=255)
                                                aCATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCAtt                                                     >  1:1428683/1‑100 (MQ=255)
                                                 cATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTc                                                    >  1:3887415/1‑100 (MQ=255)
                                                   tCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTccc                                                  >  2:1664149/1‑100 (MQ=255)
                                                    cGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTcc                                                   >  1:1683491/1‑98 (MQ=255)
                                                    cGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCt                                                 >  1:1104289/1‑100 (MQ=255)
                                                      aCGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTg                                               >  2:602966/1‑100 (MQ=255)
                                                        ggAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAg                                             >  1:620205/1‑100 (MQ=255)
                                                        ggAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAg                                             <  2:2817602/100‑1 (MQ=255)
                                                            ccGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCtt                                         >  2:1905618/1‑100 (MQ=255)
                                                              ggTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTca                                       <  2:1426230/100‑1 (MQ=255)
                                                                tGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTcacc                                     <  2:279075/100‑1 (MQ=255)
                                                                tGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTcacc                                     >  1:1482107/1‑100 (MQ=255)
                                                                 gCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTcacca                                    <  2:893756/100‑1 (MQ=255)
                                                                   ggTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCaccacc                                  <  1:2439275/100‑1 (MQ=255)
                                                                   ggTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCaccacc                                  >  2:4151230/1‑100 (MQ=255)
                                                                   ggTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCaccacc                                  >  2:3202484/1‑100 (MQ=255)
                                                                      aaaCATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTACTTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGc                               <  1:713759/100‑1 (MQ=255)
                                                                       aaCATCAATCTCACCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGcc                              <  1:1803535/100‑1 (MQ=255)
                                                                         cATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGcccc                            >  2:1317183/1‑100 (MQ=255)
                                                                         cATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGcccc                            >  2:987726/1‑100 (MQ=255)
                                                                         cATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGcccc                            >  1:3734768/1‑100 (MQ=255)
                                                                           tCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTg                          <  2:408955/100‑1 (MQ=255)
                                                                           gatggaatgCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTg                          >  1:1166122/9‑100 (MQ=37)
                                                                            cAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGc                         >  2:4397756/1‑100 (MQ=255)
                                                                            cAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGc                         >  1:3173971/1‑100 (MQ=255)
                                                                             aaTCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCt                        >  2:3758074/1‑100 (MQ=255)
                                                                              ggaatgCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTc                       >  1:850643/6‑100 (MQ=38)
                                                                                ctcGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCAc                     >  1:1556866/1‑100 (MQ=255)
                                                                                aatgCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCAc                     >  2:55468/4‑100 (MQ=38)
                                                                                aatgCCGGGTACGTCTGAGTAGTACCCGGCATTCTATCAATAGATTACAACGTAATTGTAAGCGCATACCCTTGAGACTTCACCACCAGCAcctcct                        >  1:343519/4‑97 (MQ=25)
                                                                                  cGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACtt                   <  2:4203923/100‑1 (MQ=255)
                                                                                  cGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACtt                   <  1:2787602/100‑1 (MQ=255)
                                                                                   nccGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACttt                  >  1:4325826/2‑100 (MQ=39)
                                                                                   gCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACttt                  <  2:2051791/100‑1 (MQ=255)
                                                                                    ccGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTc                 >  2:513967/1‑100 (MQ=39)
                                                                                    ccGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTc                 >  2:1556073/1‑100 (MQ=255)
                                                                                      gggTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAg               >  1:3346204/1‑100 (MQ=255)
                                                                                      gggTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAg               <  2:3425237/100‑1 (MQ=255)
                                                                                      gggTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAg               <  2:3594245/100‑1 (MQ=255)
                                                                                      gggTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAg               <  2:3734792/100‑1 (MQ=255)
                                                                                       ggTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGc              >  2:1755140/1‑100 (MQ=255)
                                                                                        gTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGcc             <  2:2227961/100‑1 (MQ=255)
                                                                                        gTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGcc             >  1:1685249/1‑100 (MQ=255)
                                                                                          aCGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTg           <  2:2084482/100‑1 (MQ=255)
                                                                                          aCGACTCAGACGTACCCGGCATCCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAACCCCTGCTCACTTCAAGCCTg           >  1:381416/1‑100 (MQ=255)
                                                                                           cGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGc          >  2:2487969/1‑100 (MQ=255)
                                                                                           cGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGc          >  2:996795/1‑100 (MQ=255)
                                                                                           cGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGc          >  2:1209709/1‑100 (MQ=255)
                                                                                            gACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCg         <  1:3644364/100‑1 (MQ=255)
                                                                                              cTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGaa       >  1:3100742/1‑100 (MQ=255)
                                                                                              cTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGaa       >  1:1732935/1‑100 (MQ=255)
                                                                                                cAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAAcc     >  1:247637/1‑100 (MQ=255)
                                                                                                cAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAAcc     >  1:2446373/1‑100 (MQ=255)
                                                                                                cAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAAcc     >  2:2857037/1‑100 (MQ=255)
                                                                                                 aGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTTATCGCATTCCCTTGAGGCTTCACCACCATCCCCTGTTCACTTTTAGGCTGCGAACCg    <  1:441709/100‑1 (MQ=255)
                                                                                                 aGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTCTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAAcc     >  2:3483218/1‑99 (MQ=255)
                                                                                                  gACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCg    <  1:1376378/99‑1 (MQ=255)
                                                                                                  gACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGt   >  2:685674/1‑100 (MQ=255)
                                                                                                   aCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTc  >  2:3428062/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
CAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTC  >  NC_000913/3832095‑3832293

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