Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,498,892 A→T 16.6% intergenic (+17/‑46) 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,8920AT16.6% 363.8 / 60.9 145intergenic (+17/‑46)ydcH/rimLDUF465 family protein/ribosomal‑protein‑L7/L12‑serine acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (67/54);  new base T (9/15);  total (76/69)
Fisher's exact test for biased strand distribution p-value = 1.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.46e-01

Only 100 of 157 total aligned reads displayed.

TGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGAA  >  NC_000913/1498793‑1498991
                                                                                                   |                                                                                                   
tGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCgcca                                                                                                     >  1:3020195/1‑100 (MQ=255)
 ggTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCgccag                                                                                                    >  1:1239147/1‑100 (MQ=255)
 ggTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCgccag                                                                                                    >  1:3253755/1‑100 (MQ=255)
 ggTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGAAGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCgccag                                                                                                    >  1:2907939/1‑100 (MQ=255)
   tCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCgccagcc                                                                                                  <  2:3020176/100‑1 (MQ=255)
    ccGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCgccagcca                                                                                                 >  2:1590804/1‑100 (MQ=255)
      gCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAgg                                                                                               >  2:2704675/1‑100 (MQ=255)
       cATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAggc                                                                                              >  1:3430524/1‑100 (MQ=255)
        aTGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAggcg                                                                                             >  2:3419064/1‑100 (MQ=255)
        aTGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAggcg                                                                                             <  1:3059250/100‑1 (MQ=255)
        aTGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAggcg                                                                                             <  1:3513664/100‑1 (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:305117/100‑2 (MQ=255)
             aaaaCAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCTGGCTGGCGGCtta                                                                                        <  1:294589/100‑2 (MQ=255)
             aaaaCAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCttt                                                                                        >  2:3552469/1‑100 (MQ=255)
              aaaCAAAAGCTACAGTTAAAAGATGAGGTGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCTGGCTGGCGGCttag                                                                                       <  2:1389058/100‑3 (MQ=255)
              aaaCAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCTGGCTGGCGGCttag                                                                                       <  1:3724763/100‑3 (MQ=255)
              aaaCAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCTGGCTGGCGGCttag                                                                                       <  1:3157226/100‑3 (MQ=255)
              aaaCAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCtttt                                                                                       <  1:1435586/100‑1 (MQ=255)
              aaaCAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAAATTTCCTCAGCCGACAGACAggcggatctt                                                                                       >  1:3866882/1‑95 (MQ=255)
               aaCAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCttttt                                                                                      >  2:300068/1‑100 (MQ=255)
               aaCAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCttttt                                                                                      <  1:917696/100‑1 (MQ=255)
                aaaaaaGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCtttttt                                                                                     >  1:3824022/3‑100 (MQ=255)
                aCAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCTGGCTTGCGGCttag                                                                                       <  2:1052719/98‑3 (MQ=255)
                 cAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTa                                                                                    >  2:319330/1‑100 (MQ=255)
                  aaaaGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTaa                                                                                   <  1:2610218/100‑1 (MQ=255)
                  aaaaGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTaa                                                                                   >  1:1638659/1‑100 (MQ=255)
                   aaaGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCTGGCTGGCGGCttaggaaag                                                                                  <  2:2689137/100‑8 (MQ=38)
                   aaaGCTACAGTTAAAAGAACAGATCATCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCTGGCTGGCGGCttaggaa                                                                                    <  2:2973811/98‑6 (MQ=37)
                       cTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACt                                                                              >  1:964049/1‑100 (MQ=255)
                        tACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTg                                                                             >  1:1794652/1‑100 (MQ=255)
                        tACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTg                                                                             <  1:294263/100‑1 (MQ=255)
                        tACAGCTAAAAGATGAGATGCTCAAAATCGTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTg                                                                             <  2:2484138/100‑1 (MQ=255)
                         aCAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCTGGCGGCTTTTTTAACAACt                                                                              >  2:2951737/1‑98 (MQ=255)
                            gTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATg                                                                         >  2:1260368/1‑100 (MQ=255)
                            gTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATg                                                                         >  1:4045854/1‑100 (MQ=255)
                             ttAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATgg                                                                        <  1:3401417/100‑1 (MQ=255)
                             ttAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATgg                                                                        >  1:1904914/1‑100 (MQ=255)
                              tAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGa                                                                       >  1:1052554/1‑100 (MQ=255)
                                 cttccgatctgCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTg                                                                    <  2:4034852/91‑1 (MQ=255)
                                     tgagatgCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTg                                                                >  2:2353378/1‑100 (MQ=255)
                                       agatgCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGTGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGga                                                              <  2:625062/100‑1 (MQ=255)
                                       agatgCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGga                                                              >  1:2464087/1‑100 (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)
                                             tCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAg                                                        >  1:443329/1‑100 (MQ=255)
                                                aaaTCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATg                                                     <  2:491274/100‑1 (MQ=255)
                                                  agccTGCGGCAGGAGAGCGTCAAAGAGGTTTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGa                                                       <  1:1787641/94‑1 (MQ=255)
                                                   tCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACt                                                  <  1:2998353/100‑1 (MQ=255)
                                                    ccTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTg                                                 <  2:3366633/100‑1 (MQ=255)
                                                      tGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGaa                                               <  1:518826/100‑1 (MQ=255)
                                                       gcagcaGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGaaa                                              >  2:1015250/1‑100 (MQ=255)
                                                        cagcagGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACGACTGCATGGATTGACTGGAGATAAGATGACTGAAAc                                             >  2:2280859/1‑100 (MQ=255)
                                                        cagcagGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAAc                                             <  1:382458/100‑1 (MQ=255)
                                                         agcagGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACg                                            <  1:3435136/100‑1 (MQ=255)
                                                          gcagGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGa                                           >  1:2187011/1‑100 (MQ=255)
                                                           cagGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGAt                                          >  1:1691678/1‑100 (MQ=255)
                                                              gagaGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATaaa                                       <  2:1961291/100‑1 (MQ=255)
                                                                gagCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATaaatg                                     <  1:721593/100‑3 (MQ=255)
                                                                gagCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAg                                     <  2:2807804/100‑1 (MQ=255)
                                                                gagCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAg                                     <  2:135048/100‑1 (MQ=255)
                                                                gagCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAg                                     <  1:1991757/100‑1 (MQ=255)
                                                                  gCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTa                                   >  1:3619423/1‑100 (MQ=255)
                                                                    gTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAg                                 <  1:1708177/100‑1 (MQ=255)
                                                                     tCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGc                                <  2:1048865/100‑1 (MQ=255)
                                                                      cAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCg                               <  2:2186998/100‑1 (MQ=255)
                                                                        aaGAGGAGTAAACTTTCATAAGCCGCCAGACAGGCGAATTTTTTAACAACTGCCTGGGAAGACCGGAGATAAGAAGATCGAAACGATAAAAGTAAGCGaa                             >  2:3500951/1‑100 (MQ=255)
                                                                          gagGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTGTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATc                           <  2:2910392/100‑1 (MQ=255)
                                                                            ggTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCAc                         >  2:4090248/1‑100 (MQ=255)
                                                                            ggTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCAc                         >  1:917641/1‑100 (MQ=255)
                                                                             gtgtAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCAc                         >  2:410628/1‑99 (MQ=255)
                                                                              tgtAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACtt                       >  1:2275431/1‑100 (MQ=255)
                                                                               gtAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTg                      >  2:2572693/1‑100 (MQ=255)
                                                                                 aaaCTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGaa                    >  2:809/1‑100 (MQ=255)
                                                                                    gttaaaaaaGCCGCCTGGCTGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTa                 >  1:1958711/8‑100 (MQ=38)
                                                                                    gttaaaaaaGCCGCCTGGCTGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTa                 >  2:2011026/8‑100 (MQ=38)
                                                                                      ttCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACa               <  1:34692/100‑1 (MQ=255)
                                                                                       tCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACAt              <  1:335391/100‑1 (MQ=255)
                                                                                       aaaaaaGCCGCCTGGCTGGCGGCTTTTTTAACAACTGCATGGCTTGACTGGCGATAAGATGACTGAACCGATAAAAGTACGCGAATCACTTGAATTACAt              >  2:1558396/5‑100 (MQ=255)
                                                                                       aaaaaaGCCGCCTGGCTGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACAt              >  2:3923237/5‑100 (MQ=255)
                                                                                        ccTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATg             >  1:3209861/1‑100 (MQ=255)
                                                                                        aaaaaGCCGCCTGGCTGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATg             >  2:4014860/4‑100 (MQ=255)
                                                                                         cTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGc            >  1:1172119/1‑100 (MQ=255)
                                                                                         cTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGc            <  1:3552492/100‑1 (MQ=255)
                                                                                         cTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGAAGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGc            >  2:2894468/1‑100 (MQ=255)
                                                                                          aaaGCCGCCTGGCTGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCt           >  2:3254140/2‑100 (MQ=255)
                                                                                            aGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGt         >  1:226541/1‑100 (MQ=255)
                                                                                              ccgccAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGGATTACATGCTGTTg       >  1:1558313/1‑100 (MQ=255)
                                                                                              ccgccAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTg       >  2:3171482/1‑100 (MQ=255)
                                                                                               cgccAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGc      >  1:2178504/1‑100 (MQ=255)
                                                                                                gccagccaGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCa     >  1:834768/1‑100 (MQ=255)
                                                                                                 ccagccagGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAg    >  2:2591952/1‑100 (MQ=255)
                                                                                                   agccagGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGaa  >  1:1519964/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TGGTCCGCATGAAAAAACAAAAGCTACAGTTAAAAGATGAGATGCTCAAAATCCTGCAGCAGGAGAGCGTCAAAGAGGTGTAAACTTTCCTAAGCCGCCAGCCAGGCGGCTTTTTTAACAACTGCATGGATTGACTGGAGATAAGATGACTGAAACGATAAAAGTAAGCGAATCACTTGAATTACATGCTGTTGCAGAA  >  NC_000913/1498793‑1498991

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: