Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,832,191 C→G 20.0% 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,1910CG20.0% 275.2 / 63.0 115intergenic (+25/+28)yicH/yicIputative inner membrane‑anchored periplasmic AsmA family protein/putative alpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base C (43/48);  new base G (9/14);  total (53/62)
Fisher's exact test for biased strand distribution p-value = 6.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

Only 100 of 117 total aligned reads displayed.

CAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACC  >  NC_000913/3832092‑3832290
                                                                                                   |                                                                                                   
 agcaAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGACGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCa                                                                                                    <  2:717521/100‑1 (MQ=255)
 agcaAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCa                                                                                                    <  2:599949/100‑1 (MQ=255)
     aGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACg                                                                                                >  1:242143/1‑100 (MQ=255)
     aGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACg                                                                                                <  2:1877660/100‑1 (MQ=255)
        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                                                                                         <  1:1427392/100‑1 (MQ=39)
            cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCgg                                                                                         <  2:1355357/100‑1 (MQ=39)
            cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCgg                                                                                         <  1:2559615/100‑1 (MQ=39)
            cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCgg                                                                                         <  1:1318947/100‑1 (MQ=39)
            cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCgg                                                                                         <  1:956317/100‑1 (MQ=255)
            cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCgg                                                                                         >  2:2446020/1‑100 (MQ=255)
            cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCgg                                                                                         <  1:2319500/100‑1 (MQ=255)
            cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCgg                                                                                         <  1:930777/100‑1 (MQ=255)
            cAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCgg                                                                                         <  1:2062920/100‑1 (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)
                ccATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCggagat                                                                                     >  2:1920404/1‑96 (MQ=38)
                ccATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCGGcgag                                                                                     <  1:121881/100‑4 (MQ=38)
                  ttccgatcTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTcc                                                                                   <  2:12263/93‑1 (MQ=255)
                   tGAATGGTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGTCTGAGTCGTACCCgg                                                                                         <  2:1078132/93‑1 (MQ=38)
                    gAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCa                                                                                       <  2:2330090/94‑1 (MQ=255)
                     aaTGCTGGCATCGTTTCCAGCTGTGAAGTTACATCGACGGAGCTGGTGCGGTAATCATCAATCTCGCCGTGTACGACTCAGACGTACCCGGCATTCCATc                                                                                <  2:1463848/100‑1 (MQ=255)
                     aaTGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATc                                                                                <  2:25346/100‑1 (MQ=255)
                     aaTGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATc                                                                                <  2:1703964/100‑1 (MQ=255)
                     aaTGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATc                                                                                <  1:1731849/100‑1 (MQ=255)
                       tGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGCCGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCa                                                                               <  2:1425505/99‑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                                                                          >  2:366220/1‑100 (MQ=255)
                           ggCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGa                                                                          >  1:1517831/1‑100 (MQ=255)
                            gCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGAt                                                                         <  1:628072/100‑1 (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)
                                              aaGTTACATCGACGGAGCTGGTGCGGTAAACATCAATATCGCGGGGTACGAGTGAGACGTAACCGGCATTCCATCAATAGATTACAACGTAATTGTCAgc                                                       >  2:2434971/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                                                >  1:2291731/1‑100 (MQ=255)
                                                     aTCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTcc                                                >  2:1936382/1‑100 (MQ=255)
                                                       cGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCt                                              >  2:695731/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)
                                                                   tGCGGTAAACATAAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACTTAATTGTCCGCGCATTCCCTTGAGGCTTca                                    >  2:289200/1‑98 (MQ=255)
                                                                      ggTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCaccacc                               >  2:278217/1‑100 (MQ=255)
                                                                      ggTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCaccacc                               <  1:2508291/100‑1 (MQ=255)
                                                                      ggTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCaccacc                               >  2:276693/1‑100 (MQ=255)
                                                                       gTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTCCAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCAccacca                              >  1:2177407/1‑100 (MQ=255)
                                                                       gTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCAccacca                              >  2:2140030/1‑100 (MQ=255)
                                                                        tAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAg                             <  2:1690407/100‑1 (MQ=255)
                                                                            cATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGcccc                         <  1:136573/100‑1 (MQ=255)
                                                                            cATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGcccc                         <  1:444805/100‑1 (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)
                                                                                     cGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACtt                <  1:1682074/100‑1 (MQ=255)
                                                                                     cGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACtt                >  1:2504773/1‑100 (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)
                                                                                       ccGGGTACGTCTGAGTCGTACCCGGCATTCCCTCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTc              >  2:106590/1‑100 (MQ=39)
                                                                                       ccGGGTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTc              >  1:1697378/1‑100 (MQ=39)
                                                                                         gggTACGTCTGAGTCGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTAACCACCAGCCCCTGCTCACTTTCa             >  2:2527663/1‑99 (MQ=255)
                                                                                          ggTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGc           >  1:1955219/1‑100 (MQ=255)
                                                                                          ggTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGc           <  2:2077460/100‑1 (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    >  1:646671/1‑100 (MQ=255)
                                                                                                 cTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGaa    >  2:2302673/1‑100 (MQ=255)
                                                                                                  tCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAAc   <  1:1672000/100‑1 (MQ=255)
                                                                                                   cAGACGTACCCGGCATTCCATCAATAGATTACAACATAATTGTCAGCGCATTCCCTTGAGGCTTCACCAAAAGCACCTGCTCACTTTCAGCCTGCGAAcc  >  1:1492946/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
CAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACC  >  NC_000913/3832092‑3832290

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: