Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,832,194 A→T 19.0% 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,1940AT19.0% 253.7 / 53.8 116intergenic (+28/+25)yicH/yicIputative inner membrane‑anchored periplasmic AsmA family protein/putative alpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base A (47/46);  new base T (8/14);  total (55/61)
Fisher's exact test for biased strand distribution p-value = 2.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.15e-01

Only 100 of 116 total aligned reads displayed.

AGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGT  >  NC_000913/3832097‑3832292
                                                                                                 |                                                                                                  
   gACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAc                                                                                               >  2:1129230/1‑100 (MQ=255)
    aCTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTAcc                                                                                              <  2:1317100/100‑1 (MQ=255)
     cTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAccc                                                                                             >  1:1699843/1‑100 (MQ=255)
     cTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAccc                                                                                             <  1:2016904/100‑1 (MQ=255)
       cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCgg                                                                                           <  2:1355357/100‑1 (MQ=39)
       cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCgg                                                                                           <  1:1427392/100‑1 (MQ=39)
       cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCgg                                                                                           <  1:930777/100‑1 (MQ=255)
       cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCgg                                                                                           <  1:956317/100‑1 (MQ=255)
       cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCgg                                                                                           <  1:2319500/100‑1 (MQ=255)
       cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCgg                                                                                           <  1:2062920/100‑1 (MQ=255)
       cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCgg                                                                                           >  2:2446020/1‑100 (MQ=255)
       cAAACCACGAATGCTGGCAGCGTTTCCAGCAGTTAAGTTACATCGACGGAGCCTGTGCGCTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCgg                                                                                           <  2:1516406/100‑1 (MQ=37)
       atctccATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCgg                                                                                           <  1:1920418/96‑1 (MQ=38)
        aaaCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCGGc                                                                                          <  2:1463047/100‑1 (MQ=39)
        aaaCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCGGc                                                                                          <  1:1711116/100‑1 (MQ=39)
        aaaCCATGAAGGCTGTCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCgg                                                                                           <  2:1750266/99‑1 (MQ=255)
         aaCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCa                                                                                         <  2:490695/100‑1 (MQ=255)
         aaCCATGAATGCTGGCATCGTTGCCAGCGGTGTAGTTACGTCGACTGAGCCGGTGCGGTAATCATCAAGCTCGCCGGGTACGTCTGAGTCGTACCCgg                                                                                           <  2:251220/98‑1 (MQ=37)
           ccATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCggagat                                                                                       >  2:1920404/1‑96 (MQ=38)
           ccATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCGGcgag                                                                                       <  1:121881/100‑4 (MQ=38)
           ccATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGCACGCCTCAGACGTACCCGGCAt                                                                                        >  2:2200672/1‑99 (MQ=255)
             ttccgatcTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTcc                                                                                     <  2:12263/93‑1 (MQ=255)
              tGAATGGTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCgg                                                                                           <  2:1078132/93‑1 (MQ=38)
              tGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCa                                                                                    <  1:1257870/100‑1 (MQ=255)
               gAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCa                                                                                         <  2:2330090/94‑1 (MQ=255)
                aaTGCTGGCATCGTTTCCAGCTGTGAAGTTACATCGACGGAGCTGGTGCGGTAATCATCAATCTCGCCGTGTACGACTCAGACGTACCCGGCATTCCATc                                                                                  <  2:1463848/100‑1 (MQ=255)
                aaTGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATc                                                                                  <  2:1703964/100‑1 (MQ=255)
                aaTGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATc                                                                                  <  2:25346/100‑1 (MQ=255)
                aaTGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATc                                                                                  <  1:1731849/100‑1 (MQ=255)
                   gCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAAt                                                                               <  2:1961581/100‑1 (MQ=255)
                     tGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAg                                                                             <  2:1841315/100‑1 (MQ=255)
                      ggCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCagatcggaa                                                                            >  1:12264/1‑92 (MQ=255)
                      ggCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGa                                                                            >  1:1517831/1‑100 (MQ=255)
                      ggCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGa                                                                            >  2:366220/1‑100 (MQ=255)
                          tCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTACACATCAATCCCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGAtt                                                                          <  2:1848272/98‑1 (MQ=255)
                            gTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACaa                                                                      >  2:2320770/1‑100 (MQ=255)
                            gTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACaa                                                                      >  1:32289/1‑100 (MQ=255)
                            gTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAAATTACa                                                                       >  2:504688/1‑99 (MQ=255)
                             tttCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAAc                                                                     <  1:1057009/100‑1 (MQ=255)
                              ttCGAGCGGTGACGTTACATCGAAGGAGCCGGTGCGGTAAACATCAAGGTCACCGCGATGGACTAAGACCTACCCGTCATTAAATCAACAGATTACAACg                                                                    >  2:246152/1‑100 (MQ=255)
                              ttCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACg                                                                    <  1:1236387/100‑1 (MQ=255)
                               tCCAGCGGTGAAGTGACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCAAGAATAGATTACAACg                                                                    >  2:2103348/1‑99 (MQ=255)
                                     ggTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGt                                                             >  2:1468592/1‑100 (MQ=255)
                                         aaGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAgc                                                         <  1:2284275/100‑1 (MQ=255)
                                          aGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAgcg                                                        >  1:2221859/1‑100 (MQ=255)
                                              aCATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCAtt                                                    >  1:1861569/1‑100 (MQ=255)
                                               cATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTc                                                   >  2:864413/1‑100 (MQ=255)
                                                aTCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTcc                                                  >  2:1936382/1‑100 (MQ=255)
                                                aTCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTcc                                                  >  1:2291731/1‑100 (MQ=255)
                                                   gACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCCATAGATTACAACGTAATTGTCAGCGCATTCCCtt                                               <  2:873509/100‑1 (MQ=255)
                                                    aCGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTg                                              >  2:49858/1‑100 (MQ=255)
                                                      ggAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAg                                            <  1:2477244/100‑1 (MQ=255)
                                                      ggAGCCGGTGCGGTAAACATCAATCTCGCCGGCTACGACTCAGCCGTACCCGGCATTCCATCAATCGATTACAACGTAATTGCCAGCGCATTCCCTTGAg                                            <  1:117073/100‑1 (MQ=255)
                                                       gAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAgg                                           <  1:2375002/100‑1 (MQ=255)
                                                       gAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAgg                                           >  2:2103694/1‑100 (MQ=255)
                                                        aGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGc                                          >  1:22261/1‑100 (MQ=255)
                                                        aGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGc                                          >  1:760954/1‑100 (MQ=255)
                                                           cGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCttn                                       <  2:920606/100‑2 (MQ=255)
                                                           cGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTc                                       <  1:2446034/100‑1 (MQ=255)
                                                            ggTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTca                                      <  1:1400437/100‑1 (MQ=255)
                                                             gTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTcac                                     <  2:1911570/100‑1 (MQ=255)
                                                              tGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTcacc                                    <  1:534456/100‑1 (MQ=255)
                                                              tGCGGTAAACATAAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACTTAATTGTCCGCGCATTCCCTTGAGGCTTca                                      >  2:289200/1‑98 (MQ=255)
                                                                 ggTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCaccacc                                 >  2:276693/1‑100 (MQ=255)
                                                                 ggTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCaccacc                                 <  1:2508291/100‑1 (MQ=255)
                                                                  gTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTCCAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCAccacca                                >  1:2177407/1‑100 (MQ=255)
                                                                  gTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCAccacca                                >  2:2140030/1‑100 (MQ=255)
                                                                       cATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGcccc                           <  1:136573/100‑1 (MQ=255)
                                                                       cATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGcccc                           <  1:444805/100‑1 (MQ=255)
                                                                        aTCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCt                          >  2:1943158/1‑100 (MQ=255)
                                                                          cAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGc                        >  1:32744/1‑100 (MQ=255)
                                                                           aaTCTCGCTGGGTACGACTCAGACGTACCCGGCATGCCATCGATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGc                        <  2:1172041/99‑1 (MQ=255)
                                                                           aaTCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCt                       >  1:1207344/1‑100 (MQ=255)
                                                                             tctcGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCa                     <  1:278221/100‑1 (MQ=255)
                                                                              ctcGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCAc                    <  2:198130/100‑1 (MQ=255)
                                                                                tgCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGCTTACAACGTAATTGTGCGCGCATTCACTTGAGGCTTCACCACCAGCCCCTGCTCACtt                  >  2:565064/2‑100 (MQ=37)
                                                                                tgCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACtt                  >  1:2395536/2‑100 (MQ=39)
                                                                                cGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACtt                  >  1:2504773/1‑100 (MQ=255)
                                                                                cGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACtt                  <  1:1682074/100‑1 (MQ=255)
                                                                                cGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACtn                  <  2:896727/100‑2 (MQ=255)
                                                                                 gCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCCCttt                 >  2:894022/1‑100 (MQ=38)
                                                                                 gCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACttt                 >  2:114292/1‑100 (MQ=255)
                                                                                 gCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACt                   <  1:2435878/98‑1 (MQ=255)
                                                                                  ccGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTc                >  1:1697378/1‑100 (MQ=39)
                                                                                   cGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTAAGACGGAATTGTCAGCGAATTCGCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCa               >  2:1011158/1‑100 (MQ=255)
                                                                                    gggTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTAACCACCAGCCCCTGCTCACTTTCa               >  2:2527663/1‑99 (MQ=255)
                                                                                     ggTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGc             >  1:1955219/1‑100 (MQ=255)
                                                                                      gTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGcc            >  2:810608/1‑100 (MQ=255)
                                                                                      gTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGcc            >  1:2545760/1‑100 (MQ=255)
                                                                                         cGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGc         >  2:1624338/1‑100 (MQ=255)
                                                                                         cGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGc         <  1:475725/100‑1 (MQ=255)
                                                                                          gACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCg        <  1:818171/100‑1 (MQ=255)
                                                                                           aCTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGa       >  2:1611083/1‑100 (MQ=255)
                                                                                            cTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCAGCACCAGCCGCTGCTCACTTTCAGCCTGCGaa      >  2:1725969/1‑100 (MQ=255)
                                                                                            cTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGaa      >  2:2302673/1‑100 (MQ=255)
                                                                                            cTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGaa      >  1:646671/1‑100 (MQ=255)
                                                                                             tCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAAc     <  1:1672000/100‑1 (MQ=255)
                                                                                              cAGACGTACCCGGCATTCCATCAATAGATTACAACATAATTGTCAGCGCATTCCCTTGAGGCTTCACCAAAAGCACCTGCTCACTTTCAGCCTGCGAAcc    >  1:1492946/1‑100 (MQ=255)
                                                                                               aGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCg   <  1:1588564/100‑1 (MQ=255)
                                                                                                gACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGt  >  2:1724473/1‑100 (MQ=255)
                                                                                                 |                                                                                                  
AGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGT  >  NC_000913/3832097‑3832292

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