Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,832,191 C→G 15.2% intergenic (+25/+28) 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,1910CG15.2% 317.3 / 49.9 119intergenic (+25/+28)yicH/yicIputative inner membrane‑anchored periplasmic AsmA family protein/putative alpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base C (47/53);  new base G (10/8);  total (58/61)
Fisher's exact test for biased strand distribution p-value = 6.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.79e-01

Only 100 of 119 total aligned reads displayed.

AGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAAC  >  NC_000913/3832093‑3832289
                                                                                                  |                                                                                                  
agcaAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCa                                                                                                   >  1:1127914/1‑100 (MQ=255)
agcaAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGGGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCTCCGGGTACGACTCa                                                                                                   >  1:448788/1‑100 (MQ=255)
  caAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGa                                                                                                 <  2:1048066/100‑1 (MQ=255)
    aGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACTGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACg                                                                                               <  2:285441/100‑1 (MQ=255)
    aGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACg                                                                                               <  1:2013092/100‑1 (MQ=255)
    aGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACg                                                                                               >  2:1443100/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)
       gACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAc                                                                                            >  1:1638103/1‑100 (MQ=255)
       gACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGATACGACTCAGACGTAc                                                                                            >  2:1305389/1‑100 (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                                                                                        <  1:422694/100‑1 (MQ=39)
           cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCgg                                                                                        <  1:66094/100‑1 (MQ=39)
           cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCgg                                                                                        <  2:243222/100‑1 (MQ=39)
           cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCgg                                                                                        <  2:825088/100‑1 (MQ=255)
           cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCgg                                                                                        >  2:1628782/1‑100 (MQ=255)
            aaaCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCgg                                                                                        <  2:870166/99‑1 (MQ=255)
             aaCCATGAATGCTGGCATCGTTTCCATCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCGGc                                                                                       <  2:1782447/99‑1 (MQ=38)
             aaCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCGGcg                                                                                      <  1:861676/100‑2 (MQ=38)
             aaCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCa                                                                                      <  2:1654383/100‑1 (MQ=255)
             aaCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCa                                                                                      >  1:422688/1‑100 (MQ=255)
              aCCAGGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCAt                                                                                     <  2:173133/100‑1 (MQ=255)
                 aTGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTcc                                                                                  <  2:1127911/100‑1 (MQ=255)
                 aTGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTcc                                                                                  >  1:1640921/1‑100 (MQ=255)
                 aTGAATGCGGGCATCGTTTCCAGCGGTGAAGTTACATGGACGGAGGCGGTGCGGTAAACATCAATCTCGCCGGGTACGATTCAGACGTACCCGGCATTcc                                                                                  <  2:1271314/100‑1 (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:1295997/1‑100 (MQ=255)
                             aTCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTa                                                                      >  1:1970343/1‑100 (MQ=255)
                             aTCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTAACCGGTATTAAATCAATAGATTa                                                                      >  1:735788/1‑100 (MQ=255)
                                gTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACaa                                                                   <  1:136568/100‑1 (MQ=255)
                                gTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACaa                                                                   <  1:238881/100‑1 (MQ=255)
                                 tttCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAAc                                                                  >  2:1346377/1‑100 (MQ=255)
                                 tttCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAAc                                                                  <  1:172657/100‑1 (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:643224/100‑1 (MQ=255)
                                              aGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAgcg                                                     <  1:1287862/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)
                                                  aCATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCAt                                                  >  1:1755610/1‑99 (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)
                                                               cGGTTCGGTAAACATCAATCCCGCCGGAAACAACTAAGACGTACCCGGCAGTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCATCAGGCtc                                     >  2:1970874/1‑98 (MQ=255)
                                                               cGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTc                                    <  2:390074/100‑1 (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                          >  2:1975540/1‑100 (MQ=255)
                                                                         aaCATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGcc                          <  1:780002/100‑1 (MQ=255)
                                                                          aCATCAATCTCGCCGGGTACGACACAGACGTACCCAACATTCCATCAATAGATTACAACGTAATTGTAAGCGCATTCCCTTGAGGCTTCACCACCAGcc                          >  2:1095033/1‑99 (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                   <  1:1628793/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)
                                                                                   tcGCCGGGTACGACTCAGACGTACCCGGCCTTCCATCCATCGCTTACAACGTAATTGTCAGCGCATTCACTTGAGGCTTCACCACCAGCCCCTGCTCACt                <  2:448789/100‑1 (MQ=255)
                                                                                   atgCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACt                >  1:1479947/3‑100 (MQ=38)
                                                                                   atgCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACt                >  1:1171527/3‑100 (MQ=38)
                                                                                    tgCCGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACtt               >  1:771076/2‑100 (MQ=39)
                                                                                    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)
                                                                                       cGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCa            >  1:146108/1‑100 (MQ=255)
                                                                                        gggTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAg           >  1:113600/1‑100 (MQ=255)
                                                                                        gggTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAg           <  1:1011424/100‑1 (MQ=255)
                                                                                        gggTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAg           <  2:69386/100‑1 (MQ=255)
                                                                                        gggTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAg           <  2:605682/100‑1 (MQ=255)
                                                                                        gggTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAg           <  2:1003684/100‑1 (MQ=255)
                                                                                         ggTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGc          <  1:724771/100‑1 (MQ=255)
                                                                                             cGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGc      <  2:2207419/100‑1 (MQ=255)
                                                                                                cTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGaa   >  2:1840992/1‑100 (MQ=255)
                                                                                                  |                                                                                                  
AGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAAC  >  NC_000913/3832093‑3832289

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: