Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,831,757 C→T 15.9% 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.9% 528.4 / 71.1 217A434V (GCC→GTC) yicHputative inner membrane‑anchored periplasmic AsmA family protein
Reads supporting (aligned to +/- strand):  ref base C (105/77);  new base T (19/16);  total (124/93)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 235 total aligned reads displayed.

TCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTAC  >  NC_000913/3831658‑3831856
                                                                                                   |                                                                                                   
tCCCGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATgc                                                                                                     >  1:1246730/1‑100 (MQ=255)
  ccGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATgccg                                                                                                   <  2:2946889/100‑1 (MQ=255)
  ccGTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATgccg                                                                                                   <  1:2872208/100‑1 (MQ=255)
    gTGGCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATgccgcc                                                                                                 <  2:1471902/100‑1 (MQ=255)
      ggCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATgccgccgc                                                                                               >  1:3415187/1‑100 (MQ=255)
       gCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAACTGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATgccgccgcc                                                                                              >  1:1110406/1‑100 (MQ=255)
       gCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATgccgccgcc                                                                                              >  2:1923495/1‑100 (MQ=255)
       gCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATgccgccgcc                                                                                              >  1:3050324/1‑100 (MQ=255)
       gCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATgccgccgcc                                                                                              <  1:251777/100‑1 (MQ=255)
         aGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGCATCTGAATGCcgccgccgc                                                                                            >  1:2848417/1‑100 (MQ=255)
           cTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGa                                                                                          >  2:2878974/1‑100 (MQ=255)
            tCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGAc                                                                                         <  1:372650/100‑1 (MQ=255)
             cACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACa                                                                                        >  1:1011319/1‑100 (MQ=255)
             cACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACa                                                                                        <  1:4080203/100‑1 (MQ=255)
                cgcgCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTc                                                                                     <  2:352512/100‑1 (MQ=255)
                 gcgcTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGTCGCGGCGGCGGCATTCa                                                                                    <  1:4340185/100‑1 (MQ=255)
                  cgcTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGTCGCGGCGGCGGCATTCag                                                                                   <  1:1761506/100‑2 (MQ=255)
                  cgcTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCaa                                                                                   <  1:832675/100‑1 (MQ=255)
                  cgcTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGCCTGGAGTGGCTCGCCGAATCTGAATGTCGCGGCGGCGGCATCCag                                                                                   <  1:3310465/100‑2 (MQ=255)
                  cgcTCGATGGTTACGCTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGTCGCGGCGGCGGCATTCag                                                                                   <  2:4343599/100‑2 (MQ=255)
                   gcTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAAt                                                                                  <  1:897145/100‑1 (MQ=255)
                   gcTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAAt                                                                                  >  1:495149/1‑100 (MQ=255)
                    cTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATc                                                                                 >  2:1023106/1‑100 (MQ=255)
                     tCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCg                                                                                <  2:3400667/100‑1 (MQ=255)
                     tCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCg                                                                                >  1:3577260/1‑100 (MQ=255)
                       gATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGTCGCGGCGGCGGCATTCagattcg                                                                              <  2:1388499/100‑7 (MQ=37)
                        aTGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCgtgt                                                                             <  2:347702/100‑1 (MQ=255)
                         tGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGtt                                                                            <  1:573086/100‑1 (MQ=255)
                          ggTTACGGTGCCAACCTGACGCTTGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGTCGCGGCGGCGGCATTCagattcgc                                                                             <  2:2570230/98‑8 (MQ=25)
                          ggTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTg                                                                           <  2:4413879/100‑1 (MQ=255)
                           gTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGa                                                                          <  2:4387007/100‑1 (MQ=255)
                           gTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGa                                                                          >  1:1341440/1‑100 (MQ=255)
                             tACGTTGCAAACCTGACGCTGGTGACCCATCATAAAGGGGGCATCAGGCGTGGCCCGGCGCCTCTTACTACCGCCGCAGCGCCATCCAATc                                                                                 >  1:2245671/1‑91 (MQ=255)
                             tACGGTGCGGACCTGGCCCTGATTAGCGATCATAAATGGGGCGTCTGGCGTGGCTCGGCAAATCTGAATGCCGCAGCCTCGACCTTGAATCGTCttta                                                                          >  1:1166796/1‑96 (MQ=255)
                             tACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGCATGCCGCCGCCGCGACATTCAATCGTGTTGATg                                                                        >  2:668067/1‑100 (MQ=255)
                              aCGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGt                                                                       <  2:1064939/100‑1 (MQ=255)
                               cGGTGCCAACCTGACGCTGGTTACCGATCATAAATTGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACAGTCAATCGTGTTGATGt                                                                       >  2:2196598/1‑99 (MQ=255)
                               cGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGtt                                                                      >  1:2982316/1‑100 (MQ=255)
                                      aaCCTGAGGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGcc                                                               >  2:2317482/1‑100 (MQ=255)
                                      aaCCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTGGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGcc                                                               >  1:2292022/1‑100 (MQ=255)
                                        ccTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCg                                                             <  2:1332060/100‑1 (MQ=255)
                                          tGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTc                                                           >  2:2184415/1‑100 (MQ=255)
                                             cGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCt                                                        >  2:2287250/1‑100 (MQ=255)
                                             cGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCt                                                        >  2:3055119/1‑100 (MQ=255)
                                               cTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTgc                                                      <  2:1834151/100‑2 (MQ=255)
                                                 ggTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGgcg                                                    >  1:1733502/1‑100 (MQ=255)
                                                    tACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTg                                                 <  1:4122132/100‑1 (MQ=255)
                                                      ccGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGAc                                               >  2:3575065/1‑100 (MQ=255)
                                                      ccGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGAc                                               >  2:2387534/1‑100 (MQ=255)
                                                      ccGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGAc                                               <  2:3949609/100‑1 (MQ=255)
                                                      ccGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGAc                                               <  1:1503117/100‑1 (MQ=255)
                                                        gATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCg                                             >  2:4215082/1‑100 (MQ=255)
                                                        gATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCg                                             >  2:724074/1‑100 (MQ=255)
                                                        gATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGGTGATGTTCGTCGCCCGTCGCTGGCGCTGACCg                                             >  1:1655290/1‑100 (MQ=255)
                                                          tCATAAATGGGGCGTCTGGAGTGGCTCGTCGAATCTGAATGCCGCCGCCGCGATATGCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGa                                                <  1:2407955/95‑1 (MQ=255)
                                                          tCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGcc                                           <  2:3920930/100‑1 (MQ=255)
                                                          tCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGcc                                           <  2:2279643/100‑1 (MQ=255)
                                                           cATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCCCTGACCGcc                                           >  1:2596262/1‑99 (MQ=255)
                                                            aTAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCaa                                         >  2:155354/1‑100 (MQ=255)
                                                            aTAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCaa                                         >  2:1755561/1‑100 (MQ=255)
                                                              aaaTGGGGGGGCTGGAGTGGCTCGGCGAATCTGAATGCCGCGGCCGCGACATTCAATCGTGTGGATGTTCGTCGCCCGTCGCTGGCGGTCAGCGCCagc                                        <  1:1522878/99‑3 (MQ=255)
                                                               aaTGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAAcag                                      <  2:820500/100‑1 (MQ=255)
                                                               aaTGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAAcag                                      >  2:2639184/1‑100 (MQ=255)
                                                                 tGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAAcagca                                    <  2:2810595/100‑1 (MQ=255)
                                                                    ggCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAgcagca                                 >  2:375990/1‑100 (MQ=255)
                                                                       gTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACgg                              <  2:543295/100‑1 (MQ=255)
                                                                       gTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACgg                              >  2:1990571/1‑100 (MQ=255)
                                                                       gTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACgg                              >  2:940352/1‑100 (MQ=255)
                                                                        tCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGt                             >  1:3367805/1‑100 (MQ=255)
                                                                        tCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGt                             >  2:2736483/1‑100 (MQ=255)
                                                                         cTGGAGTGGCTCGGCGAATCTGAATTCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTg                            <  1:3902062/100‑1 (MQ=255)
                                                                         cTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTg                            <  2:3347853/100‑1 (MQ=255)
                                                                         cTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTg                            <  2:1727054/100‑1 (MQ=255)
                                                                         cTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCATCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTg                            >  1:465137/1‑100 (MQ=255)
                                                                          tGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGa                           >  1:4394441/1‑100 (MQ=255)
                                                                               tGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATc                      <  1:1286879/100‑1 (MQ=255)
                                                                                ggCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCa                     >  1:1937844/1‑100 (MQ=255)
                                                                                 gCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAg                    <  2:4108217/100‑1 (MQ=255)
                                                                                 gCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAg                    >  2:2136217/1‑100 (MQ=255)
                                                                                  cTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGGCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGc                   <  2:956140/100‑1 (MQ=255)
                                                                                   tCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCg                  <  2:2619470/100‑1 (MQ=255)
                                                                                    cGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACAGCCAACAGCAGCACGGTGAATATCAGCg                  >  1:355616/1‑99 (MQ=255)
                                                                                      gCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAAc               <  1:4180205/100‑1 (MQ=255)
                                                                                      acacgattGAATGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAAc               >  2:2463629/8‑100 (MQ=37)
                                                                                        gAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTg             >  1:1431168/1‑100 (MQ=255)
                                                                                        gAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTg             >  2:1146932/1‑100 (MQ=255)
                                                                                         cgattGAATGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGa            >  1:703631/5‑100 (MQ=38)
                                                                                         aaTCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGa            >  2:4070032/1‑100 (MQ=255)
                                                                                          tcttGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAg           <  1:3255202/97‑1 (MQ=255)
                                                                                          aTCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAg           >  1:2195947/1‑100 (MQ=255)
                                                                                          aTCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAg           <  1:2140421/100‑1 (MQ=255)
                                                                                          aTCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCCCCCGTCGCTGGCGCTGACCGCCACAAGCAGCACGGTGAATACCAGCGAACTg             >  2:4123807/1‑98 (MQ=255)
                                                                                           attGAATGTCGCGGCGGCGGCATTCAATCGTGTTGATTGTCGTCGCGCGTCGCTGGCGCTGACCGACAACAGCAGAACGGGGAATATCAGCGAACt              >  1:2324139/3‑96 (MQ=21)
                                                                                           attGAATGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGt          >  2:626080/3‑100 (MQ=38)
                                                                                           attGAATGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGt          >  2:2142095/3‑100 (MQ=38)
                                                                                            cTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTg         <  2:3191967/100‑1 (MQ=255)
                                                                                            cTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGAAGCACGGTGAATATCAGCGAACTGAGTg         >  1:3951393/1‑100 (MQ=255)
                                                                                             tGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGc        >  1:1846790/1‑100 (MQ=255)
                                                                                              gAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCa       >  1:2046994/1‑100 (MQ=255)
                                                                                              gAATGCCGCCGCCGCGACATGCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCa       >  1:2698733/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: