Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,831,764 C→G 16.0% A436A (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,7640CG16.0% 588.1 / 90.2 215A436A (GCC→GCGyicHputative inner membrane‑anchored periplasmic AsmA family protein
Reads supporting (aligned to +/- strand):  ref base C (102/76);  new base G (18/16);  total (123/92)
Fisher's exact test for biased strand distribution p-value = 7.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 234 total aligned reads displayed.

GCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAA  >  NC_000913/3831665‑3831863
                                                                                                   |                                                                                                   
gCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATgccgccgcc                                                                                                     <  1:251777/100‑1 (MQ=255)
gCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATgccgccgcc                                                                                                     >  2:1923495/1‑100 (MQ=255)
 cAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGccgccgccg                                                                                                    >  1:3270345/1‑100 (MQ=255)
    cTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGa                                                                                                 >  1:1151091/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)
       aCCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCCATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACa                                                                                               >  1:345557/1‑99 (MQ=255)
        ccGCGCTCGATGGTTACGGTGCCAACCTGACTCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGATCGGCGAATCTGAATGCCGCCGCCGCGACAt                                                                                              >  1:1552969/1‑99 (MQ=255)
        ccGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACAtt                                                                                             >  1:2233638/1‑100 (MQ=255)
          gcgcTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCa                                                                                           <  2:1608580/100‑1 (MQ=255)
          gcgcTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCa                                                                                           >  1:4356820/1‑100 (MQ=255)
           cgcTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGTCGCGGCGGCGGCATTCa                                                                                           <  1:702854/99‑1 (MQ=255)
           cgcTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGCCTGGAGTGGCTCGCCGAATCTGAATGTCGCGGCGGCGGCATCCag                                                                                          <  1:3310465/100‑2 (MQ=255)
            gcTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGTCGCGGCGGCGGCATTCaga                                                                                         <  1:1423007/100‑3 (MQ=38)
             cTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGTCGCGGCGGCGGCATTCagat                                                                                        <  1:56149/100‑4 (MQ=38)
              tCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCTTCTGGAGTGGCTCGGCGAATCTGAATGTCGCGGCGGCGGCATTCagat                                                                                        <  1:4084690/99‑4 (MQ=38)
               cGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGTCGCGGCGGCGGCATTCagattc                                                                                      <  2:1342957/100‑6 (MQ=37)
                gATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCgtg                                                                                     <  2:3050342/100‑1 (MQ=255)
                 aTGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCgtgt                                                                                    >  1:4148352/1‑100 (MQ=255)
                   ggTTACGGTGCCAATCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTg                                                                                  <  1:3009610/100‑1 (MQ=255)
                   ggTTACGGTGCCAACCTGACGCTTGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGTCGCGGCGGCGGCATTCagattcgc                                                                                    <  2:2570230/98‑8 (MQ=25)
                   ggTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTg                                                                                  >  2:4310271/1‑100 (MQ=255)
                    gTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGa                                                                                 >  2:2504410/1‑100 (MQ=255)
                    gTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGa                                                                                 >  1:1341440/1‑100 (MQ=255)
                    gTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGa                                                                                 <  1:2339124/100‑1 (MQ=255)
                      tACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGCATGCCGCCGCCGCGACATTCAATCGTGTTGATg                                                                               >  2:668067/1‑100 (MQ=255)
                         ggTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTtc                                                                            <  2:2235751/100‑1 (MQ=255)
                             ccAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTtcgtcg                                                                        >  2:2144224/1‑100 (MQ=255)
                               aaCCTGAGGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGcc                                                                      >  2:2317482/1‑100 (MQ=255)
                                 ccTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCg                                                                    <  2:1332060/100‑1 (MQ=255)
                                 ccTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCg                                                                    >  1:3392472/1‑100 (MQ=255)
                                    gACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCg                                                                 <  2:495150/100‑1 (MQ=255)
                                        cTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTgg                                                             >  2:2196457/1‑100 (MQ=255)
                                        cTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTgg                                                             <  2:1874493/100‑1 (MQ=255)
                                          ggTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACCTTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGgc                                                            >  1:2194731/1‑99 (MQ=255)
                                           gTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGgcgc                                                          <  1:2287623/100‑1 (MQ=255)
                                               ccGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGAc                                                      >  2:2291976/1‑100 (MQ=255)
                                               ccGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGAc                                                      >  2:2387534/1‑100 (MQ=255)
                                               ccGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGAc                                                      <  2:3949609/100‑1 (MQ=255)
                                                cGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGAcc                                                     >  1:2119617/1‑100 (MQ=255)
                                                cGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGAcc                                                     >  2:3379355/1‑100 (MQ=255)
                                                cGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGAcc                                                     <  1:453232/100‑1 (MQ=255)
                                                cGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGAcc                                                     >  2:2140428/1‑100 (MQ=255)
                                                 gATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCg                                                    >  1:3576229/1‑100 (MQ=255)
                                                 gATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCg                                                    >  2:724074/1‑100 (MQ=255)
                                                 gATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCg                                                    >  2:4215082/1‑100 (MQ=255)
                                                 gATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACAGTCAATCGTGCTGATGTTCGCCGCCCGTCGCTTGCGCTGAcc                                                     >  1:2952496/1‑99 (MQ=255)
                                                   tCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGcc                                                  <  2:3920930/100‑1 (MQ=255)
                                                   tCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGcc                                                  <  1:106195/100‑1 (MQ=255)
                                                    cATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCa                                                 <  2:3411711/100‑1 (MQ=255)
                                                     aTAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCaa                                                >  2:1755561/1‑100 (MQ=255)
                                                        aaTGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAAcag                                             <  2:820500/100‑1 (MQ=255)
                                                          tGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAAcagca                                           >  1:3940192/1‑100 (MQ=255)
                                                          tGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAAcagca                                           <  2:2810595/100‑1 (MQ=255)
                                                             ggCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAgcagca                                        >  2:375990/1‑100 (MQ=255)
                                                                gTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCACCAGTAGCACgg                                     <  1:2318150/100‑1 (MQ=255)
                                                                gTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACgg                                     >  2:940352/1‑100 (MQ=255)
                                                                 tCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGt                                    >  2:2736483/1‑100 (MQ=255)
                                                                  cTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTg                                   <  2:1727054/100‑1 (MQ=255)
                                                                  cTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTg                                   <  1:2025346/100‑1 (MQ=255)
                                                                   tGGAGTGGCTCGGCGAATCTGCATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACgg                                     >  1:1749652/1‑97 (MQ=255)
                                                                      aGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAata                               <  1:1854114/100‑1 (MQ=255)
                                                                        tGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATc                             <  1:1286879/100‑1 (MQ=255)
                                                                          gCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAg                           <  2:4108217/100‑1 (MQ=255)
                                                                          gCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAg                           <  2:4394473/100‑1 (MQ=255)
                                                                          gCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAg                           >  2:4345403/1‑100 (MQ=255)
                                                                           cTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAg                           <  1:1200969/99‑1 (MQ=255)
                                                                           cTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGGCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGc                          <  2:956140/100‑1 (MQ=255)
                                                                            tCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCg                         >  2:1080357/1‑100 (MQ=255)
                                                                            tCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCg                         <  2:2619470/100‑1 (MQ=255)
                                                                               gCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAAc                      <  1:2640209/100‑1 (MQ=255)
                                                                               gCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAAc                      >  2:2279984/1‑100 (MQ=255)
                                                                                cGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACt                     >  1:2856268/1‑100 (MQ=255)
                                                                                 gAATCTGAATGCCTCCGCCGCGACATTCAATCGTGTTGAGGTTCGTCGCCCGTCTTTCGCGCCGACCGCCAACAGCAGCACGGTTAATATCAGCGAACTg                    <  1:993279/100‑1 (MQ=255)
                                                                                  aaTCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGa                   >  2:561856/1‑100 (MQ=255)
                                                                                  aaTCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGa                   >  2:4070032/1‑100 (MQ=255)
                                                                                   tattGAATGTCGCGGCGGAGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGGGAATATCAGCGAACTGAg                  >  1:1500942/4‑100 (MQ=255)
                                                                                   aTCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAg                  <  1:1668926/100‑1 (MQ=255)
                                                                                   aTCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCCCCCGTCGCTGGCGCTGACCGCCACAAGCAGCACGGTGAATACCAGCGAACTg                    >  2:4123807/1‑98 (MQ=255)
                                                                                    attGAATGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAg                  >  1:1827840/3‑99 (MQ=38)
                                                                                    attGAATGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGt                 >  2:3535721/3‑100 (MQ=38)
                                                                                    attGAATGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGt                 >  1:892588/3‑100 (MQ=38)
                                                                                     cTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTg                <  2:3191967/100‑1 (MQ=255)
                                                                                      tGAATGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGc               >  2:2035913/1‑100 (MQ=255)
                                                                                      tGAATGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGc               >  1:3810298/1‑100 (MQ=255)
                                                                                      tGAATGCCGCCGGCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGc               >  2:2891130/1‑100 (MQ=255)
                                                                                      tGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAgaga               >  2:3255220/1‑97 (MQ=255)
                                                                                      tGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGc               >  2:1050233/1‑100 (MQ=255)
                                                                                      tGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCCGTGCCCGCCTACAGCAGCACGGTGAATATCAGAGCACTAAGGGc               >  1:4245274/1‑100 (MQ=255)
                                                                                      tGAATGCCGCCGCCGCGACATACAATCGTGTTGATGTTCGTAGACCATCGCTGGCGCTGCCCGCCCACAGCAGCATGGTCAATATCAGCGAACTGAGTga               >  1:1516453/1‑99 (MQ=255)
                                                                                       gAATGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCa              >  1:492705/1‑100 (MQ=255)
                                                                                       gAATGCCGCCGCCGCGACATGCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCa              >  1:2698733/1‑100 (MQ=255)
                                                                                        aaTGTCGCGGCGGCGGCATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCAt             >  2:1393402/1‑100 (MQ=255)
                                                                                         aTGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCAtt            <  2:3512217/100‑1 (MQ=255)
                                                                                            ccgccgccgCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTAc         <  2:2700899/100‑1 (MQ=255)
                                                                                             cgccgccgcGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACt        <  1:324957/100‑1 (MQ=255)
                                                                                             cgccgccgcGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGGGAATATCAGCGAACTGAGTGCATTTACt        >  2:1981599/1‑100 (MQ=255)
                                                                                              gccgccgcGACATTCAATCGTGTTGGTGTTCGTCGCCCGTCGCTGGCGCTGTCCGCCAATAGGAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTg       <  2:3015396/100‑1 (MQ=255)
                                                                                               ccgccgcGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTgn      <  1:728327/100‑2 (MQ=255)
                                                                                                   |                                                                                                   
GCAGCTCACCGCGCTCGATGGTTACGGTGCCAACCTGACGCTGGTTACCGATCATAAATGGGGCGTCTGGAGTGGCTCGGCGAATCTGAATGCCGCCGCCGCGACATTCAATCGTGTTGATGTTCGTCGCCCGTCGCTGGCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAA  >  NC_000913/3831665‑3831863

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: