Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,498,894 C→G 16.6% intergenic (+19/‑44) ydcH → / → rimL DUF465 family protein/ribosomal‑protein‑L7/L12‑serine acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,498,8940CG16.6% 387.4 / 60.4 145intergenic (+19/‑44)ydcH/rimLDUF465 family protein/ribosomal‑protein‑L7/L12‑serine acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (65/54);  new base G (9/15);  total (76/69)
Fisher's exact test for biased strand distribution p-value = 1.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 155 total aligned reads displayed.

TCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGAAAA  >  NC_000913/1498796‑1498993
                                                                                                  |                                                                                                   
    cATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAggc                                                                                                >  1:3430524/1‑100 (MQ=255)
     aTGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAggcg                                                                                               <  1:3513664/100‑1 (MQ=255)
     aTGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAggcg                                                                                               <  1:3419084/100‑1 (MQ=255)
     aTGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAggcg                                                                                               <  1:3059250/100‑1 (MQ=255)
     aTGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAggcg                                                                                               >  2:3419064/1‑100 (MQ=255)
      tGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAggcgg                                                                                              <  1:1565319/100‑1 (MQ=255)
       gAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAggcggc                                                                                             >  2:1762245/1‑100 (MQ=255)
        aaaaaaCAAAAGCTACAGTTAAACGATGAGCTGCTCAAAATCCTGCAGCAGGAGAGCGTCAGCGAGGTGTAAACCTTCCTAAGCCGCCAGCGAggggg                                                                                              >  2:1243783/1‑98 (MQ=255)
        aaaaaaCAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCTGGCTGGCGGCt                                                                                            <  2:420175/100‑1 (MQ=255)
        aaaaaaCAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCt                                                                                            >  1:138338/1‑100 (MQ=255)
         aaaaaCAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCtt                                                                                           >  2:2230410/1‑100 (MQ=255)
          aaaaCAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCTGGCTGGCGGCtta                                                                                          <  1:294589/100‑2 (MQ=255)
          aaaaCAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCTGGCTGGCGGCtta                                                                                          <  1:3156357/100‑2 (MQ=255)
          aaaaCAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCttt                                                                                          >  2:3552469/1‑100 (MQ=255)
            aaCAAAAGTTACAGTTAAAAGATGAGATGTTCAAGAGCTTGCATTAGGAGAGCTTCAAAGAGGTGTAAACTTTCCTAAGCCGCCTGGCTGGCGGCttag                                                                                         <  2:3025003/99‑3 (MQ=255)
            aaCAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCttttt                                                                                        >  2:300068/1‑100 (MQ=255)
              cAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCTGGCTGGCGGCttaggaa                                                                                      <  1:1586694/100‑6 (MQ=38)
              cAAAAGCTACAGTGAAAAGATGAGGTGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCTGGCTGGCGGCttaggaa                                                                                      <  2:3254074/100‑6 (MQ=38)
               aaaaGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTaa                                                                                     <  1:5297/100‑1 (MQ=255)
               aaaaGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTaa                                                                                     <  1:2610218/100‑1 (MQ=255)
                aaaGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCTGGCTGGCGGCttaggaaag                                                                                    <  2:2689137/100‑8 (MQ=38)
                   gCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTaacaac                                                                                 >  1:813722/1‑100 (MQ=255)
                    cTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACt                                                                                >  1:964049/1‑100 (MQ=255)
                     tACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTg                                                                               <  1:294263/100‑1 (MQ=255)
                     tACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTg                                                                               >  1:1794652/1‑100 (MQ=255)
                     tACAGCTAAAAGATGAGATGCTCAAAATCGTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTg                                                                               <  2:2484138/100‑1 (MQ=255)
                      aCAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCTGGCGGCTTTTTTAACAACt                                                                                >  2:2951737/1‑98 (MQ=255)
                       cAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCa                                                                             >  2:4016573/1‑100 (MQ=255)
                         gTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATg                                                                           >  2:1260368/1‑100 (MQ=255)
                           tAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGa                                                                         >  1:1052554/1‑100 (MQ=255)
                            aaaaGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGa                                                                         <  1:3870540/99‑1 (MQ=255)
                            aaaaGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGAt                                                                        >  1:945715/1‑100 (MQ=255)
                              cttccgatctgCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTg                                                                      <  2:4034852/91‑1 (MQ=255)
                                 atgagatgCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACt                                                                   <  1:1659469/100‑1 (MQ=255)
                                  tgagatgCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTg                                                                  >  2:2353378/1‑100 (MQ=255)
                                   gagatgCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTgg                                                                 <  1:1981357/100‑1 (MQ=255)
                                   gagatgCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTgg                                                                 >  1:374946/1‑100 (MQ=255)
                                    agatgCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGTGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGga                                                                <  2:625062/100‑1 (MQ=255)
                                     gatgCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGgag                                                               <  1:582563/100‑1 (MQ=255)
                                     gatgCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGgag                                                               <  2:1179963/100‑1 (MQ=255)
                                       tgCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGAt                                                             >  2:636430/1‑100 (MQ=255)
                                        gCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGagatcggaag                                                            >  1:4034878/1‑91 (MQ=255)
                                         cTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATaa                                                           >  2:2676807/1‑100 (MQ=255)
                                          tCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAg                                                          >  1:443329/1‑100 (MQ=255)
                                          tCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAg                                                          >  2:2998334/1‑100 (MQ=255)
                                           cAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGa                                                         <  1:2196136/100‑1 (MQ=255)
                                            aaaaTCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGAt                                                        >  1:2640074/1‑100 (MQ=255)
                                             aaaTCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATg                                                       <  2:491274/100‑1 (MQ=255)
                                                tCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACt                                                    <  1:2998353/100‑1 (MQ=255)
                                                 ccTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTg                                                   <  2:3366633/100‑1 (MQ=255)
                                                    gcagcaGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGaaa                                                >  1:3026853/1‑100 (MQ=255)
                                                     cagcagGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACGACTGCATGGATTGACTGGAGATAAGATGACTGAAAc                                               >  2:2280859/1‑100 (MQ=255)
                                                     cagcagGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAAc                                               <  1:382458/100‑1 (MQ=255)
                                                        cagGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGAt                                            >  1:1691678/1‑100 (MQ=255)
                                                           gagaGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATaaa                                         <  2:1961291/100‑1 (MQ=255)
                                                             gagCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAg                                       <  1:2149609/100‑1 (MQ=255)
                                                             gagCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAg                                       <  2:135048/100‑1 (MQ=255)
                                                             gagCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAg                                       <  2:2807804/100‑1 (MQ=255)
                                                               gCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTa                                     >  1:3619423/1‑100 (MQ=255)
                                                               gCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTa                                     >  1:1544757/1‑100 (MQ=255)
                                                                 gTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAg                                   >  1:3852887/1‑100 (MQ=255)
                                                                  tCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGc                                  <  2:1048865/100‑1 (MQ=255)
                                                                   cAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCg                                 >  2:2756250/1‑100 (MQ=255)
                                                                   cAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCg                                 <  2:2186998/100‑1 (MQ=255)
                                                                     aaGAGGAGTAAACTTTCATAAGCCGCCAGACAGGCGAATTTTTTAACAACTGCCTGGGAAGACCGGAGATAAGAAGATCGAAACGATAAAAGTAAGCGaa                               >  2:3500951/1‑100 (MQ=255)
                                                                       gagGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTGTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATc                             <  2:2910392/100‑1 (MQ=255)
                                                                         ggTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCAc                           >  1:917641/1‑100 (MQ=255)
                                                                         ggTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCAc                           >  2:4090248/1‑100 (MQ=255)
                                                                          gtgtAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCAc                           >  2:410628/1‑99 (MQ=255)
                                                                          gtgtAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACt                          >  1:3193433/1‑100 (MQ=255)
                                                                            gtAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTg                        >  2:2572693/1‑100 (MQ=255)
                                                                              aaaCTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGaa                      >  2:809/1‑100 (MQ=255)
                                                                              aaaCTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGaa                      >  2:110158/1‑100 (MQ=255)
                                                                                aCTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAAtt                    <  1:2280872/100‑1 (MQ=255)
                                                                                 gttaaaaaaGCCGCCTGGCTGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTa                   >  2:2011026/8‑100 (MQ=38)
                                                                                 gttaaaaaaGCCGCCTGGCTGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTa                   >  1:1958711/8‑100 (MQ=38)
                                                                                  ttaaaaaaGCCGCCTGGCTGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTAc                  >  1:1397170/7‑100 (MQ=38)
                                                                                   ttCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACa                 <  1:34692/100‑1 (MQ=255)
                                                                                   ttCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACa                 <  2:2907923/100‑1 (MQ=255)
                                                                                    tCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACAt                <  1:335391/100‑1 (MQ=255)
                                                                                    aaaaaaGCCGCCTGGCTGGCGGCTTTTTTAACAACTGCATGGCTTGACTGGCGATAAGATGACTGAACCGATAAAAGTACGCGAATCACTTGAATTACAt                >  2:1558396/5‑100 (MQ=255)
                                                                                    aaaaaaGCCGCCTGGCTGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGACACGATAAAAGTAAGCGAATCACTTGAATTACAt                >  2:1027348/5‑100 (MQ=255)
                                                                                    aaaaaaGCCGCCTGGCTGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACAt                >  2:3923237/5‑100 (MQ=255)
                                                                                     ccTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATg               >  1:3209861/1‑100 (MQ=255)
                                                                                      cTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGc              <  1:3552492/100‑1 (MQ=255)
                                                                                      cTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGc              >  1:4121326/1‑100 (MQ=255)
                                                                                      cTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGc              <  2:1717055/100‑1 (MQ=255)
                                                                                      cTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGc              >  1:1172119/1‑100 (MQ=255)
                                                                                      aaaaGCCGCCTGGCTGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCAATTGAATTACATGc              >  1:3447387/3‑100 (MQ=255)
                                                                                       aaaGCCGCCTGGCTGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCt             >  2:3254140/2‑100 (MQ=255)
                                                                                         aGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGt           >  1:226541/1‑100 (MQ=255)
                                                                                           ccgccAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTg         >  2:3171482/1‑100 (MQ=255)
                                                                                            cgccAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGc        >  2:1553154/1‑100 (MQ=255)
                                                                                            cgccAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGc        >  1:2178504/1‑100 (MQ=255)
                                                                                            cgccAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGc        <  1:1835204/100‑1 (MQ=255)
                                                                                             gccagccaGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCa       >  1:834768/1‑100 (MQ=255)
                                                                                             gccagccaGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCa       >  1:1767717/1‑100 (MQ=255)
                                                                                                agccagGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGaa    >  1:1519964/1‑100 (MQ=255)
                                                                                                  ccagGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAATATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGaaaa  >  1:3687904/1‑100 (MQ=255)
                                                                                                  ccagGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGaaaa  <  2:2578247/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
TCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGAAAA  >  NC_000913/1498796‑1498993

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: