Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,832,188 A→T 15.0% intergenic (+22/+31) 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,1880AT15.0% 298.4 / 46.2 120intergenic (+22/+31)yicH/yicIputative inner membrane‑anchored periplasmic AsmA family protein/putative alpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base A (48/54);  new base T (10/8);  total (58/62)
Fisher's exact test for biased strand distribution p-value = 6.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.43e-02

Only 100 of 120 total aligned reads displayed.

AGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGC  >  NC_000913/3832090‑3832285
                                                                                                  |                                                                                                 
agcagcAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGAc                                                                                                  <  2:976665/100‑1 (MQ=255)
agcagcAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGAc                                                                                                  >  1:1397789/1‑100 (MQ=255)
 gcagcaAGTGACTCAAACCATGAATGCCGGCATCGTTGCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGAc                                                                                                  <  2:1008767/99‑1 (MQ=255)
   agcaAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCa                                                                                               >  1:1127914/1‑100 (MQ=255)
   agcaAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGGGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCTCCGGGTACGACTCa                                                                                               >  1:448788/1‑100 (MQ=255)
     caAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGa                                                                                             <  2:126758/100‑1 (MQ=255)
     caAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGa                                                                                             <  2:1048066/100‑1 (MQ=255)
     caAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACAGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGa                                                                                             >  2:1432801/1‑100 (MQ=255)
       aGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACg                                                                                           >  2:1446103/1‑100 (MQ=255)
       aGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACg                                                                                           <  1:1319037/100‑1 (MQ=255)
        gTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGt                                                                                          <  2:1963627/100‑1 (MQ=255)
        gTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGt                                                                                          >  1:2136620/1‑100 (MQ=255)
          gACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAc                                                                                        >  1:1638103/1‑100 (MQ=255)
          gACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGATACGACTCAGACGTAc                                                                                        >  2:1305389/1‑100 (MQ=255)
            cTCAATTCATGACTGTTGGCATCGTTTCCAGCGGGGAAGTTTCATTGACGGAGCGGGTGCGGTAAACATCAATCTAGCCGGGAACGACTCCGACGt                                                                                          >  2:445480/1‑96 (MQ=255)
            cTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAccc                                                                                      <  1:2087275/100‑1 (MQ=255)
            cTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAccc                                                                                      >  1:2207434/1‑100 (MQ=255)
            cTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAcac                                                                                      >  2:1918704/1‑98 (MQ=255)
             tCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCg                                                                                     <  2:618390/100‑1 (MQ=255)
             tCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCg                                                                                     >  2:2010098/1‑100 (MQ=255)
              cAAACCATGAATGCTGGCATCGTTTCCAGCGTTAAAGATACATCGCCGTAGCCGGTGCGGTATACTTCAATCTCGCCGGTTACGACTCAGACGTACCCgg                                                                                    <  2:1936229/100‑1 (MQ=255)
              cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCgg                                                                                    <  1:1585821/100‑1 (MQ=39)
              cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCgg                                                                                    <  2:243222/100‑1 (MQ=39)
              cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCgg                                                                                    <  1:422694/100‑1 (MQ=39)
              cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCgg                                                                                    >  2:1628782/1‑100 (MQ=255)
              cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCgg                                                                                    <  1:1198987/100‑1 (MQ=255)
              cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCgg                                                                                    <  2:825088/100‑1 (MQ=255)
               aaaCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGc                                                                                   <  2:1970769/100‑1 (MQ=255)
                aaCCATGAATGCTGGCATCGTTTCCATCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCGGc                                                                                   <  2:1782447/99‑1 (MQ=38)
                aaCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCa                                                                                  <  2:1654383/100‑1 (MQ=255)
                 aCCAGGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCAt                                                                                 <  2:173133/100‑1 (MQ=255)
                    aTGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTcc                                                                              <  2:1127911/100‑1 (MQ=255)
                    aTGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTcc                                                                              >  1:1640921/1‑100 (MQ=255)
                     tGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCa                                                                             >  1:1691132/1‑100 (MQ=255)
                        aTGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCa                                                                          <  2:316573/100‑1 (MQ=255)
                         tgttGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCa                                                                          <  1:53215/96‑1 (MQ=255)
                               cATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGAtt                                                                   >  1:1003686/1‑100 (MQ=255)
                                aTCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTa                                                                  >  1:1970343/1‑100 (MQ=255)
                                aTCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTa                                                                  >  1:1295997/1‑100 (MQ=255)
                                aTCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAACCGGTATTAAATCAATAGATTa                                                                  >  1:735788/1‑100 (MQ=255)
                                   gTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACaa                                                               <  1:238881/100‑1 (MQ=255)
                                    tttCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAAc                                                              <  1:172657/100‑1 (MQ=255)
                                    tttCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAAc                                                              >  2:1346377/1‑100 (MQ=255)
                                       ccAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTa                                                           <  2:2207432/100‑1 (MQ=255)
                                         aGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAAt                                                         >  2:983838/1‑100 (MQ=255)
                                         aGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACATACCCGGCATTCCATCAATAGATTACAACGTAAt                                                         >  2:1179577/1‑100 (MQ=255)
                                           cGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTg                                                       <  2:2231827/100‑1 (MQ=255)
                                           cGGTGAAGTTACATCGACGGAGCCGGTGCAGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTg                                                       <  2:1286898/100‑1 (MQ=255)
                                            ggTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGt                                                      >  2:1404022/1‑100 (MQ=255)
                                                 aGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAgcg                                                 <  1:1287862/100‑1 (MQ=255)
                                                 aGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAgcg                                                 <  1:643224/100‑1 (MQ=255)
                                                    tACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCAt                                              >  2:343010/1‑100 (MQ=255)
                                                    tACATCGACGGAGCCGGTGCGGTAAACATCAAGCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCAt                                              <  1:2017347/100‑1 (MQ=255)
                                                     aCATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCAtt                                             <  2:1638092/100‑1 (MQ=255)
                                                      cATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTc                                            >  2:1116153/1‑100 (MQ=255)
                                                        tCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTccc                                          <  2:870337/100‑1 (MQ=255)
                                                        tCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTccc                                          >  2:970591/1‑100 (MQ=255)
                                                        tCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACAACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTccc                                          >  2:1763980/1‑100 (MQ=255)
                                                          gACGGAGCCGGTGCGGTAAACATCAATGTCGCCGAGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCt                                         >  2:2050261/1‑99 (MQ=255)
                                                           aCGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCtt                                        <  1:913820/99‑1 (MQ=255)
                                                           aCGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTg                                       <  2:328588/100‑1 (MQ=255)
                                                               aGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGc                                   <  1:1922046/100‑1 (MQ=255)
                                                                gCCTGTGCTGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCt                                  <  1:1582058/100‑1 (MQ=255)
                                                                gCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCt                                  >  1:1883162/1‑100 (MQ=255)
                                                                  cGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTc                                <  2:390074/100‑1 (MQ=255)
                                                                   ggTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCCTCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGCGGCTTca                               >  1:771709/1‑100 (MQ=255)
                                                                   ggTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTca                               >  1:170185/1‑100 (MQ=255)
                                                                      gCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTcacca                            >  1:1783338/1‑100 (MQ=255)
                                                                       cGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTcaccac                           >  2:1121940/1‑100 (MQ=255)
                                                                        ggTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCaccacc                          >  1:1434099/1‑100 (MQ=255)
                                                                        ggTAAACATCAATCTCGCAAGGTACGACTCAGACGTAATCAGCATATCATTAATCGATTACAACGTAATTTTCAGCGCATTCCCTTGATGCTTca                               >  2:1720979/1‑95 (MQ=255)
                                                                            aaCATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGcc                      <  1:1432811/100‑1 (MQ=255)
                                                                            aaCATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGcc                      <  1:780002/100‑1 (MQ=255)
                                                                            aaCATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGcc                      >  2:1975540/1‑100 (MQ=255)
                                                                              cATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGcccc                    >  2:874872/1‑100 (MQ=255)
                                                                               aTCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCt                   <  1:1305395/100‑1 (MQ=255)
                                                                                tcattttctccGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATGACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGcccc                    <  1:1772826/89‑1 (MQ=255)
                                                                                 cAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGc                 <  2:181811/100‑1 (MQ=255)
                                                                                   aTCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTc               <  2:868644/100‑1 (MQ=255)
                                                                                   aTCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTc               >  2:1483620/1‑100 (MQ=255)
                                                                                    tctcGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCa              <  2:1849129/100‑1 (MQ=255)
                                                                                     ctcGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCAc             >  1:1458734/1‑100 (MQ=255)
                                                                                      atgCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACt            >  1:1479947/3‑100 (MQ=38)
                                                                                       tgCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACtt           >  1:771076/2‑100 (MQ=39)
                                                                                       tgCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACtt           >  1:499156/2‑100 (MQ=255)
                                                                                       cGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACtt           <  1:1537167/100‑1 (MQ=255)
                                                                                        gCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACttt          >  2:2065081/1‑100 (MQ=39)
                                                                                         ccGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTc         >  2:2223150/1‑100 (MQ=39)
                                                                                         ccGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTc         >  2:247805/1‑100 (MQ=255)
                                                                                          cGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCa        >  1:1071101/1‑100 (MQ=255)
                                                                                          cGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCa        <  1:1483630/100‑1 (MQ=255)
                                                                                           gggTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAg       <  1:1011424/100‑1 (MQ=255)
                                                                                           gggTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAg       >  1:113600/1‑100 (MQ=255)
                                                                                           gggTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAg       <  2:69386/100‑1 (MQ=255)
                                                                                           gggTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAg       <  2:605682/100‑1 (MQ=255)
                                                                                           gggTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAg       >  2:2146948/1‑100 (MQ=255)
                                                                                           gggTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAg       <  2:1003684/100‑1 (MQ=255)
                                                                                            ggTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGc      <  1:724771/100‑1 (MQ=255)
                                                                                             gTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAgac     >  2:1053514/1‑98 (MQ=255)
                                                                                                cGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGc  <  1:1379613/100‑1 (MQ=255)
                                                                                                  |                                                                                                 
AGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGC  >  NC_000913/3832090‑3832285

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: