Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,831,761 C→G 15.5% A435A (GCC→GCG 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,7610CG15.5% 496.9 / 77.3 182A435A (GCC→GCGyicHputative inner membrane‑anchored periplasmic AsmA family protein
Reads supporting (aligned to +/- strand):  ref base C (82/70);  new base G (17/11);  total (100/82)
Fisher's exact test for biased strand distribution p-value = 5.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

Only 100 of 196 total aligned reads displayed.

TGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGA  >  NC_000913/3831663‑3831859
                                                                                                  |                                                                                                  
tGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATgccgccg                                                                                                   >  1:1069405/1‑100 (MQ=255)
   cAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGccgccgccg                                                                                                >  1:2346100/1‑100 (MQ=255)
       tCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGAc                                                                                            <  2:674495/100‑1 (MQ=255)
       tCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGAc                                                                                            <  1:1561824/100‑1 (MQ=255)
         aCCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACAt                                                                                          <  1:2300275/100‑1 (MQ=255)
           cgcgTTCGATGGTTACGGTGCCAACCGGACGCTGGTTACCGATCAGACATGGCGCGTCTGGAGTGGCTCGGCGAATCCGAATGCCGCCTCCGCGACATTc                                                                                        <  2:2778354/100‑1 (MQ=255)
              gcTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGTCGCGGCGGCGGCATTCag                                                                                      <  2:1181106/99‑2 (MQ=39)
               cTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGTCGCGGCGGCGGCATTCagat                                                                                    <  2:1134528/100‑4 (MQ=38)
               cTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGTCGCGGCGGCGGCATTCagat                                                                                    <  1:565327/100‑4 (MQ=38)
               cTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGTCGCGGCGGCGGCATTCagat                                                                                    <  2:130293/100‑4 (MQ=38)
               cTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGTCGCGGCGGCGGCATTCagat                                                                                    <  1:1432675/100‑4 (MQ=38)
               cTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATc                                                                                    >  1:867744/1‑100 (MQ=255)
                 cGATGGTTACGGTGCTAACCTGACGCTGTTTATCGATCATAAATGGGGCGTCTGGAGTCGCTCGGCGAATCTGAATGTCGCGGCGGCGGCATTCagattc                                                                                  <  2:1854125/100‑6 (MQ=25)
                 cGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGTCGCGGCGGCGGCATTCagattc                                                                                  <  1:2824019/100‑6 (MQ=37)
                  gATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCgtg                                                                                 >  2:3225755/1‑100 (MQ=255)
                   aTGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTAGAGTGGCTCGGCGAATCTGAATGTCGCGGCGGCGGCATTCagattcgc                                                                                <  1:1090826/100‑8 (MQ=37)
                     ggTTACGGTGCCAACCTGACGCTGGTTACCTATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGtt                                                                               <  2:3191799/99‑1 (MQ=255)
                     ggTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTg                                                                              <  2:425475/100‑1 (MQ=255)
                      gTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGa                                                                             >  1:3224209/1‑100 (MQ=255)
                       ttACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGAt                                                                            >  2:2025291/1‑100 (MQ=255)
                        tACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATg                                                                           >  1:2168080/1‑100 (MQ=255)
                        tACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATg                                                                           >  2:1728860/1‑100 (MQ=255)
                        tACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATg                                                                           <  2:3621810/100‑1 (MQ=255)
                         aCGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGt                                                                          >  1:202788/1‑100 (MQ=255)
                          cGGTGCCAACCTTACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGtt                                                                         >  2:1286366/1‑100 (MQ=255)
                          cGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGCTGtt                                                                         >  1:1871569/1‑100 (MQ=255)
                              gCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTtcgtc                                                                     <  2:3481023/100‑1 (MQ=255)
                               ccAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTtcgtcg                                                                    <  1:3560927/100‑1 (MQ=255)
                               ccAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTtcgtcg                                                                    >  1:3682729/1‑100 (MQ=255)
                                cAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGc                                                                   >  2:1917056/1‑100 (MQ=255)
                                 aaCCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGcc                                                                  <  1:2910928/100‑1 (MQ=255)
                                    cTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCg                                                                <  1:1596009/99‑1 (MQ=255)
                                       aCGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGGTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGc                                                            >  2:1538498/1‑100 (MQ=255)
                                       aCGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCGGCCGCCGCCACATTCAATCGTGTTGATGTTACTCGCCCGTCg                                                             >  2:3083498/1‑99 (MQ=255)
                                        cGCTGGTTACCTATCATAAATGGGCCGTCTGTAGTGGCTCTTCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTCATCTTCGTCGCCCGTCg                                                             <  2:1625978/98‑1 (MQ=255)
                                        cGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGc                                                            <  2:2336319/99‑1 (MQ=255)
                                        cGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCt                                                           >  1:2539078/1‑100 (MQ=255)
                                          cTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTgg                                                         <  2:2434621/100‑1 (MQ=255)
                                          cTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTgg                                                         >  1:1168880/1‑100 (MQ=255)
                                           tGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGgc                                                        <  1:2395811/100‑1 (MQ=255)
                                             gTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGgcgc                                                      <  1:1956941/100‑1 (MQ=255)
                                                 ccGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGAc                                                  >  2:2730034/1‑100 (MQ=255)
                                                 ccGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCAACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGAc                                                  >  2:742462/1‑100 (MQ=255)
                                                  cGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGAcc                                                 >  1:1082426/1‑100 (MQ=255)
                                                  cGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGAcc                                                 >  1:1246075/1‑100 (MQ=255)
                                                   naTCATAAATGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGc                                               >  2:954518/2‑100 (MQ=255)
                                                   gATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGGCCg                                                >  2:1369912/1‑100 (MQ=255)
                                                   gATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGAcc                                                 >  2:1539503/1‑99 (MQ=255)
                                                   gATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCg                                                <  2:1069396/100‑1 (MQ=255)
                                                   gATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCg                                                >  2:3271888/1‑100 (MQ=255)
                                                   gATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCg                                                <  1:673802/100‑1 (MQ=255)
                                                      cATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCa                                             >  2:1214751/1‑100 (MQ=255)
                                                      cATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCa                                             <  2:1674241/100‑1 (MQ=255)
                                                        tAAATGGGGCGTCTGGAGTGGCTCGGCGAAGCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAAc                                           >  2:2258303/1‑100 (MQ=255)
                                                        naaaTGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATGCAAGCGTGTTTATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAAc                                           >  2:1511825/2‑100 (MQ=255)
                                                         aaaTGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAAca                                          >  2:1838486/1‑100 (MQ=255)
                                                           aTGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAAcagc                                        >  1:1719405/1‑100 (MQ=255)
                                                            tGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAAcagca                                       <  1:2210177/100‑1 (MQ=255)
                                                            tGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAAcagca                                       >  2:855864/1‑100 (MQ=255)
                                                              gggCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACagcagc                                     <  2:1481709/100‑1 (MQ=255)
                                                              gggCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACagcagc                                     <  2:2956908/100‑1 (MQ=255)
                                                               ggCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTAGTCGCCCGTCGCTGGCGCTGACCGCCAACAgcagca                                    >  1:3643133/1‑100 (MQ=255)
                                                                  gTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACgg                                 >  2:816845/1‑100 (MQ=255)
                                                                   tCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGt                                >  1:3104465/1‑100 (MQ=255)
                                                                    cTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTg                               <  1:784547/100‑1 (MQ=255)
                                                                      ggAGTGGCTCGGCGACTCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGaa                             <  1:3763815/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)
                                                                       gAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAat                            >  1:2933089/1‑100 (MQ=255)
                                                                       gAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAat                            <  2:2850817/100‑1 (MQ=255)
                                                                        aGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAata                           <  2:3104441/100‑1 (MQ=255)
                                                                        aGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAata                           <  1:1532935/100‑1 (MQ=255)
                                                                           ggCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCa                        >  2:2300891/1‑100 (MQ=255)
                                                                           ggCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCa                        >  2:3570546/1‑100 (MQ=255)
                                                                           ggCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCa                        >  1:2631590/1‑100 (MQ=255)
                                                                            gCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAg                       >  1:2887753/1‑100 (MQ=255)
                                                                             cTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGc                      <  1:1761223/100‑1 (MQ=255)
                                                                              tCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCg                     <  2:2989452/100‑1 (MQ=255)
                                                                              tCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCg                     <  2:1343332/100‑1 (MQ=255)
                                                                               caacacgattGAATGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGa                    >  1:1369744/10‑100 (MQ=37)
                                                                                aacacgattGAATGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGaa                   >  1:976703/9‑100 (MQ=25)
                                                                                   gAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAAc                  >  2:1793974/1‑98 (MQ=255)
                                                                                   gAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTg                >  2:2643161/1‑100 (MQ=255)
                                                                                     gattGAATGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAg              >  2:905494/4‑100 (MQ=38)
                                                                                     aTCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAg              <  2:1796888/100‑1 (MQ=255)
                                                                                     aTCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAg              <  2:2168066/100‑1 (MQ=255)
                                                                                       cTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTg            >  1:1982296/1‑100 (MQ=255)
                                                                                       cTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTg            <  2:3324220/100‑1 (MQ=255)
                                                                                        tGAATGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGc           >  2:68024/1‑100 (MQ=255)
                                                                                         gAATGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCa          >  1:3177634/1‑100 (MQ=255)
                                                                                         gAATGTCGCGGCGGCGGCATTCAATCGTGTTGACGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCa          >  2:2684160/1‑100 (MQ=255)
                                                                                          aaTGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCAt         >  1:496084/1‑100 (MQ=255)
                                                                                           aTGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCAtt        >  1:1702099/1‑100 (MQ=255)
                                                                                            tGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCAttt       >  2:3149930/1‑100 (MQ=255)
                                                                                             gccgccgccGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTa      >  2:178348/1‑100 (MQ=255)
                                                                                               cgccgccgcGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACt    <  2:3652170/100‑1 (MQ=255)
                                                                                               cgccgccgcGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACt    <  1:2125425/100‑1 (MQ=255)
                                                                                               cgccgccgcGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACt    >  1:2474747/1‑100 (MQ=255)
                                                                                                 ccgccgcGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGa  <  2:2206781/100‑1 (MQ=255)
                                                                                                  |                                                                                                  
TGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGA  >  NC_000913/3831663‑3831859

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: