Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,831,757 C→T 15.3% A434V (GCC→GTC)  yicH → putative inner membrane‑anchored periplasmic AsmA family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,831,7570CT15.3% 445.6 / 59.5 183A434V (GCC→GTC) yicHputative inner membrane‑anchored periplasmic AsmA family protein
Reads supporting (aligned to +/- strand):  ref base C (84/71);  new base T (17/11);  total (101/82)
Fisher's exact test for biased strand distribution p-value = 5.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.99e-01

Only 100 of 198 total aligned reads displayed.

TCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTAC  >  NC_000913/3831658‑3831856
                                                                                                   |                                                                                                   
tCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATgc                                                                                                     <  1:593576/100‑1 (MQ=255)
tCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATgc                                                                                                     >  2:1776116/1‑100 (MQ=255)
 cccGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATgcc                                                                                                    >  2:3428588/1‑100 (MQ=255)
  ttccgatctcTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATgccg                                                                                                   <  2:750526/91‑1 (MQ=255)
     tGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATgccgccg                                                                                                >  1:1069405/1‑100 (MQ=255)
      ggCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATgccgccgc                                                                                               >  2:677418/1‑100 (MQ=255)
          gCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCcgcg                                                                                           >  1:1617045/1‑100 (MQ=255)
           cTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCcgcg                                                                                           <  2:3628768/99‑1 (MQ=255)
                cgcgTTCGATGGTTACGGTGCCAACCGGACGCTGGTTACCGATCAGACATGGCGCGTCTGGAGTGGCTCGGCGAATCCGAATGCCGCCTCCGCGACATTc                                                                                     <  2:2778354/100‑1 (MQ=255)
                 gcgcTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCa                                                                                    >  1:2460302/1‑100 (MQ=255)
                  cgcTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGTCGCGGCGGCGGCATTCag                                                                                   <  1:1390106/100‑2 (MQ=255)
                  cgcTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCa                                                                                    >  2:2662056/1‑99 (MQ=255)
                   gcTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGTCGCGGCGGCGGCATTCag                                                                                   <  2:1181106/99‑2 (MQ=39)
                   gcTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAAt                                                                                  >  1:2890050/1‑100 (MQ=255)
                   gcTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAAt                                                                                  <  2:1784697/100‑1 (MQ=255)
                    cTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGTCGCGGCGGCGGCATTCagat                                                                                 <  2:1134528/100‑4 (MQ=38)
                    cTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGTCGCGGCGGCGGCATTCagat                                                                                 <  2:130293/100‑4 (MQ=38)
                     tCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCg                                                                                <  2:1841772/100‑1 (MQ=255)
                       gATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCgtg                                                                              <  1:3171124/100‑1 (MQ=255)
                       gATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCgtg                                                                              >  2:3225755/1‑100 (MQ=255)
                        aTGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTAGAGTGGCTCGGCGAATCTGAATGTCGCGGCGGCGGCATTCagattcgc                                                                             <  1:1090826/100‑8 (MQ=37)
                         tGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGtt                                                                            >  1:2238533/1‑100 (MQ=255)
                         tGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGtt                                                                            <  1:1718562/100‑1 (MQ=255)
                          ggTTACGGTGCCAACCTGACGCTGGTTACCTATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGtt                                                                            <  2:3191799/99‑1 (MQ=255)
                          ggTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTg                                                                           <  1:2876528/100‑1 (MQ=255)
                          ggTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTg                                                                           <  2:2170221/100‑1 (MQ=255)
                          ggTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTg                                                                           <  2:425475/100‑1 (MQ=255)
                           gTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGa                                                                          >  1:1949141/1‑100 (MQ=255)
                             tACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATg                                                                        >  1:2168080/1‑100 (MQ=255)
                             tACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATg                                                                        >  2:1728860/1‑100 (MQ=255)
                               cGGTGCCAACCTTACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGtt                                                                      >  2:1286366/1‑100 (MQ=255)
                               cGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGtt                                                                      <  1:3570273/100‑1 (MQ=255)
                               cGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGtt                                                                      >  1:3650224/1‑100 (MQ=255)
                                 gTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGACTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTtc                                                                     >  2:1886232/1‑99 (MQ=255)
                                         cTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGt                                                            >  1:1796898/1‑100 (MQ=255)
                                           gACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTGCGTCGACCGTCg                                                          >  1:1442281/1‑100 (MQ=255)
                                            aCGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGGTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGc                                                         >  2:1538498/1‑100 (MQ=255)
                                             cGCTGGTTACCTATCATAAATGGGCCGTCTGTAGTGGCTCTTCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTCATCTTCGTCGCCCGTCg                                                          <  2:1625978/98‑1 (MQ=255)
                                             cGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGc                                                         <  2:2336319/99‑1 (MQ=255)
                                             cGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCt                                                        >  1:2539078/1‑100 (MQ=255)
                                              gCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTg                                                       <  2:2975960/100‑1 (MQ=255)
                                               cTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTgg                                                      <  2:2819226/100‑1 (MQ=255)
                                               cTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTgg                                                      >  1:1168880/1‑100 (MQ=255)
                                               cTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTgg                                                      <  2:2434621/100‑1 (MQ=255)
                                                tGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGTCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGgc                                                     <  2:1968544/100‑1 (MQ=255)
                                                    tACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCTGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTg                                                 >  1:3102747/1‑100 (MQ=255)
                                                    tACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTg                                                 <  1:3216411/100‑1 (MQ=255)
                                                    tACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTg                                                 <  2:3609775/100‑1 (MQ=255)
                                                      ccGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGAc                                               <  1:822554/100‑1 (MQ=255)
                                                      ccGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGAc                                               >  2:2730034/1‑100 (MQ=255)
                                                       cGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGAcc                                              >  1:1082426/1‑100 (MQ=255)
                                                       cGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGAcc                                              >  1:1246075/1‑100 (MQ=255)
                                                       cGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGAc                                               <  2:1794885/99‑1 (MQ=255)
                                                        naTCATAAATGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGc                                            >  2:954518/2‑100 (MQ=255)
                                                        gATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGGCCg                                             >  2:1369912/1‑100 (MQ=255)
                                                        gATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGAcc                                              >  2:1539503/1‑99 (MQ=255)
                                                        gATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCg                                             <  1:673802/100‑1 (MQ=255)
                                                        gATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCg                                             <  2:2921023/100‑1 (MQ=255)
                                                        gATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCg                                             >  2:3271888/1‑100 (MQ=255)
                                                        gATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCg                                             <  2:1069396/100‑1 (MQ=255)
                                                          tCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCTGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGcc                                           <  2:3102723/100‑1 (MQ=255)
                                                           cATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCa                                          >  2:1214751/1‑100 (MQ=255)
                                                           cATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCa                                          <  2:1674241/100‑1 (MQ=255)
                                                             tAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAAc                                        >  1:1999691/1‑100 (MQ=255)
                                                             tAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAAc                                        >  1:3741449/1‑100 (MQ=255)
                                                             tAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAAc                                        <  1:1728870/100‑1 (MQ=255)
                                                             tAAATGGGGCGTCTGGAGTGGCTCGGCGAAGCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAAc                                        >  2:2258303/1‑100 (MQ=255)
                                                             naaaTGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATGCAAGCGTGTTTATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAAc                                        >  2:1511825/2‑100 (MQ=255)
                                                              aaaTGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAAc                                        <  1:1564195/99‑1 (MQ=255)
                                                               aaTGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAAcag                                      <  2:3116078/100‑1 (MQ=255)
                                                                aTGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCGGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGGTCAGCGCCagcgac                                     <  2:3073461/100‑6 (MQ=255)
                                                                aTGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAAcagc                                     >  1:1543512/1‑100 (MQ=255)
                                                                 tGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAAcagca                                    >  2:855864/1‑100 (MQ=255)
                                                                   gggCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACagcagc                                  <  2:1481709/100‑1 (MQ=255)
                                                                       gTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACgg                              <  1:914202/100‑1 (MQ=255)
                                                                       gTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACgg                              >  2:816845/1‑100 (MQ=255)
                                                                          tGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGa                           <  2:707623/100‑1 (MQ=255)
                                                                          tGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGa                           <  1:2758288/100‑1 (MQ=255)
                                                                           ggAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGGCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTg                            >  2:1242384/1‑98 (MQ=255)
                                                                           ggAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGaa                          >  1:2555454/1‑100 (MQ=255)
                                                                           ggAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGaa                          >  2:1019416/1‑100 (MQ=255)
                                                                           ggAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGaa                          >  1:1980522/1‑100 (MQ=255)
                                                                                ggCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCa                     >  2:2300891/1‑100 (MQ=255)
                                                                                ggCTCGACGCATCTGAATGCCGCCGCCGCGACATTCAACCGAGTTGATGTTCGTCGCCCGTCGCTGGCGCTGAACGCCAACAGCAGCACGGTGCata                        >  2:1836767/1‑97 (MQ=255)
                                                                                  cTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGc                   >  1:752699/1‑100 (MQ=255)
                                                                                   tCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCg                  <  2:1343332/100‑1 (MQ=255)
                                                                                     aacacgattGAATGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGaa                >  1:976703/9‑100 (MQ=25)
                                                                                         cgattGAATGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGa            >  2:3496602/5‑100 (MQ=38)
                                                                                         cgattGAATGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGa            >  1:2916141/5‑100 (MQ=38)
                                                                                         aaTCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGa            >  1:1337962/1‑100 (MQ=255)
                                                                                          gattGAATGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAg           >  2:905494/4‑100 (MQ=38)
                                                                                          aTCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAg           <  1:2337679/100‑1 (MQ=255)
                                                                                           attGAATGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGt          >  1:518856/3‑100 (MQ=38)
                                                                                            cTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTg         <  2:3324220/100‑1 (MQ=255)
                                                                                             tGAATGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGc        >  1:3342327/1‑100 (MQ=255)
                                                                                             tGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCGACAGCAGCACGGTGAATATCAGCGAACTGAGTGc        >  1:2039223/1‑100 (MQ=255)
                                                                                              gAATGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCa       >  2:2048715/1‑100 (MQ=255)
                                                                                                 tGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCAttt    >  2:3149930/1‑100 (MQ=255)
                                                                                                   ccgccgccgCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTAc  <  1:2373803/100‑1 (MQ=255)
                                                                                                   ccgccgccgCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTAc  >  1:40248/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTAC  >  NC_000913/3831658‑3831856

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: