Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,021,640 +C coding (137/366 nt) fliO → flagellar biosynthesis protein FliO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,021,6361.C98.6% 860.6 / 3.2 210coding (133/366 nt)fliOflagellar biosynthesis protein FliO
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base C (97/110);  total (99/111)
Fisher's exact test for biased strand distribution p-value = 6.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 242 total aligned reads displayed.

GCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑G‑‑CCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCTTCACTGGGCGCGCGTGAAAGGGTTGTGGTGGTCGATGTGGAAGATGCACGGC  >  NC_000913/2021536‑2021735
                                                                                                       |                                                                                                    
  cGGTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑G‑‑cc                                                                                                   >  2:3068429/1‑101 (MQ=255)
   ggTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑G‑‑ccc                                                                                                  <  1:3068161/101‑1 (MQ=255)
   ggTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑G‑‑c                                                                                                    <  1:3068080/99‑1 (MQ=255)
   ggTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGAGATTT‑‑G‑‑ccc                                                                                                  >  1:3068432/1‑101 (MQ=255)
    gTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑G‑‑cccc                                                                                                 >  2:3068435/1‑101 (MQ=255)
    gTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑G‑‑cccc                                                                                                 >  1:3068436/1‑101 (MQ=255)
    gTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑G‑‑cccc                                                                                                 <  2:3067918/101‑1 (MQ=255)
       tCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑G‑‑cccccta                                                                                              <  2:3067781/101‑4 (MQ=255)
       tCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑G‑‑cccccta                                                                                              <  1:3068063/101‑4 (MQ=255)
       tCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑G‑‑cccccct                                                                                              <  1:3067962/101‑4 (MQ=255)
        ctgctgCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑G‑‑ccccctaa                                                                                             >  2:3068439/1‑97 (MQ=255)
         tgctgcGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑G‑‑ccccctaaa                                                                                            >  1:3068442/1‑96 (MQ=255)
             gcgcCACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCTAAACGCa                                                                                         >  2:3068448/1‑101 (MQ=255)
              cgcCACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCCTAAACGCa                                                                                         >  2:3068450/1‑101 (MQ=255)
               gcCACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCTAAACGCACt                                                                                       <  2:3067857/101‑1 (MQ=255)
                ccACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCCTAAACGCACt                                                                                       <  1:3068012/101‑1 (MQ=255)
                ccACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCTAAACGCACTg                                                                                      <  2:3067990/101‑1 (MQ=255)
                ccACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCTAAACGCACTg                                                                                      <  1:3067969/101‑1 (MQ=255)
                 cACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCCTAAACGCACTg                                                                                      >  1:3068452/1‑101 (MQ=255)
                 cACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCTAAACGCACTgg                                                                                     <  2:3068140/101‑1 (MQ=255)
                   ctgctgCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCCTAAACGCACTGGc                                                                                    <  2:3068123/101‑1 (MQ=255)
                    tgctgcAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCCTAAACGCACTGGCg                                                                                   >  1:3068456/1‑101 (MQ=255)
                     gctgcAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCTAAACGCACTGGCGtt                                                                                 >  2:3068457/1‑101 (MQ=255)
                      ctgcAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCCTAAACGCACTGGCGtt                                                                                 >  2:3068459/1‑101 (MQ=255)
                      ctgcAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCTAAACGCACTGGCGTTa                                                                                <  2:3068050/101‑1 (MQ=255)
                       tgcAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCCTAAACGCACTGGCGTTa                                                                                <  2:3068204/101‑1 (MQ=255)
                       tgcAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCCTAAACGCACTGGCGTTa                                                                                <  2:3068045/101‑1 (MQ=255)
                        gcAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCCTAAACGCACTGGCGTTaa                                                                               <  1:3067920/101‑1 (MQ=255)
                         cAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCTAAACGCACTGGCGTTAACg                                                                             >  2:3068462/1‑101 (MQ=255)
                          aGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCCTAAACGCACTGGCGTTAACg                                                                             >  1:3068464/1‑101 (MQ=255)
                            gTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCCTAAACGCACTGGCGTTAACGGt                                                                           <  2:3068115/101‑1 (MQ=255)
                             tGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCTAAACGCACTGGCGTTAACGGTCt                                                                         >  1:3068468/1‑101 (MQ=255)
                              gAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCCTAAACGCACTGGCGTTAACGGTCt                                                                         <  2:3068073/101‑1 (MQ=255)
                               aGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCCTAAACGCACTGGCGTTAACGGTCTg                                                                        >  2:3068471/1‑101 (MQ=255)
                               aGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCTAAACGCACTGGCGTTAACGGTCTGa                                                                       >  1:3068470/1‑101 (MQ=255)
                                 cggcgCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCCTAAACGCACTGGCGTTAACGGTCTGaa                                                                      <  2:3068402/101‑1 (MQ=255)
                                 cggcgCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCCTAAACGCACTGGCGTTAACGGTCTGaa                                                                      >  1:3068474/1‑101 (MQ=255)
                                 cggcgCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCCTAAACGCACTGGCGTTAACGGTCTGaa                                                                      >  2:3068475/1‑101 (MQ=255)
                                   gcgcACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCCTAAACGCACTGGCGTTAACGGTCTGa                                                                       <  1:3068133/98‑1 (MQ=255)
                                   gcgcACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAAt                                                                   <  2:3067996/101‑1 (MQ=255)
                                    cgcACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAAtt                                                                  <  2:3068090/101‑1 (MQ=255)
                                       aCTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAg                                                                >  1:3068480/1‑101 (MQ=255)
                                        cTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAgc                                                               >  2:3068481/1‑101 (MQ=255)
                                         tCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAgcgc                                                             <  1:3067724/101‑1 (MQ=255)
                                           aTCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGcc                                                            >  2:3068484/1‑101 (MQ=255)
                                            tCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCa                                                           <  2:3068310/101‑1 (MQ=255)
                                             cGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAg                                                          <  2:3068066/101‑1 (MQ=255)
                                             cGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCTAAACGCAc                                                                                        >  1:3068486/1‑70 (MQ=255)
                                              gCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTg                                                        <  1:3068212/101‑1 (MQ=255)
                                              gCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTg                                                        <  1:3068053/101‑1 (MQ=255)
                                                 attattGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCtt                                                     >  1:3068489/1‑101 (MQ=255)
                                                  ttattGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCtt                                                     >  1:3068495/1‑101 (MQ=255)
                                                  ttattGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCtt                                                     <  1:3068046/101‑1 (MQ=255)
                                                  ttattGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCTTc                                                    >  1:3068492/1‑101 (MQ=255)
                                                  ttattGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCTTc                                                    >  2:3068490/1‑101 (MQ=255)
                                                   tattGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCTTc                                                    <  2:3068216/101‑1 (MQ=255)
                                                   tattGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCTTc                                                    <  1:3067982/101‑1 (MQ=255)
                                                     ttGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCTTCAc                                                  <  1:3067766/101‑1 (MQ=255)
                                                       gCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCTTCACTgg                                               <  2:3068173/101‑1 (MQ=255)
                                                        cccTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCTTCACTgg                                               <  2:3068096/101‑1 (MQ=255)
                                                        cccTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCTTCACTg                                                >  1:3068498/1‑99 (MQ=255)
                                                         ccTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCTTCACTGGgc                                             <  1:3067932/101‑1 (MQ=255)
                                                          cTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCTTCACTGGgc                                             >  1:3068501/1‑101 (MQ=255)
                                                          cTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCTTCACTGGgcg                                            <  1:3068030/101‑1 (MQ=255)
                                                             aTCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCTTCACTGGgcgcg                                          >  1:3068504/1‑101 (MQ=255)
                                                               ccTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGAGCCAGTGCTTCACTGGGCGCGCGt                                       >  2:3068505/1‑101 (MQ=255)
                                                                cTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAgc                                                               >  1:3068509/1‑77 (MQ=255)
                                                                cTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAgagc                                                             >  2:3068508/1‑78 (MQ=255)
                                                                cTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCTTCACTGGGCGCGCGTg                                      <  1:3068172/101‑1 (MQ=255)
                                                                 tCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCTTCACTGGGCGCGCGTg                                      <  2:3067899/101‑1 (MQ=255)
                                                                  cGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCTTCACTGGGCGCGCGTGaa                                    <  2:3068187/101‑1 (MQ=255)
                                                                   gctgctGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCTTCACTGGGCGCGCGTGaa                                    >  2:3068514/1‑101 (MQ=255)
                                                                   gctgctGCCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCTTCACTGGGCGCGCGTGaaa                                   <  2:3068105/101‑1 (MQ=255)
                                                                     tgctgcCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCTTCACTGGGCGCGCGTGAAAg                                  <  1:3068120/101‑1 (MQ=255)
                                                                     tgctgcCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCTTCACTGGGCGCGCGTGAAAg                                  >  2:3068516/1‑101 (MQ=255)
                                                                     tgctgcCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCTTCACTGGGCGCGCGTGAAAgg                                 <  2:3068002/101‑1 (MQ=255)
                                                                      gctgcCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCTTCACTGGGCGCGCGTGAAAggg                                <  2:3068144/101‑1 (MQ=255)
                                                                       ctgcCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCTTCACTGGGCGCGCGTGAAAggg                                >  2:3068519/1‑101 (MQ=255)
                                                                        tgcCTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCCTAAACGCACTgg                                                                                     <  1:3067910/47‑1 (MQ=37)
                                                                          ccTGGCTGGTAAAACGGTTGGGATTT‑‑GCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCTTCACTGGGCGCGCGTGAAAGGGTtgt                            >  2:3068521/1‑101 (MQ=255)
                                                                           ctggctggTAAAACGGTTGGGATTT‑‑GCCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCTTCACTGGGCGCGCGTGAAAGGGTtgt                            <  2:3068136/101‑1 (MQ=255)
                                                                           ctggctggTAAAACGGTTGGGATTT‑‑GCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCTTCACTGGGCGCGCGTGAAAGGGTtgtg                           <  2:3068021/101‑1 (MQ=255)
                                                                            tggctggTAAAACGGTTGGGATTT‑‑GCCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCTTCACTGGGCGCGCGTGAAAGGGTtgtg                           <  1:3068109/101‑1 (MQ=255)
                                                                             ggctggTAAAACGGTTGGGATTT‑‑GCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCTTCACTGGGCGCGCGTGAAAGGGTTgtggt                         >  2:3068523/1‑101 (MQ=255)
                                                                              gctggTAAAACGGTTGGGATTT‑‑GCCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCTTCACTGGGCGCGCGTGAAAGGGTTgtggt                         <  1:3068308/101‑1 (MQ=255)
                                                                               ctggTAAAACGGTTGGGATTT‑‑GCCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCTTCACTGGGCGCGCGTGAAAGGGTTgtggtg                        <  1:3067843/101‑1 (MQ=255)
                                                                               ctggTAAAACGGTTGGGATTT‑‑GCCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCTTCACTGGGCGCGCGTGAAAGGGTTgtggtg                        >  1:3068525/1‑101 (MQ=255)
                                                                                tggTAAAACGGTTGGGATTT‑‑GCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCTTCACTGGGCGCGCGTGAAAGGGTTGTggtggt                      <  1:3068037/101‑1 (MQ=255)
                                                                                  gTAAAACGGTTGGGATTT‑‑GCCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCTTCACTGGGCGCGCGTGAAAGGGTTGTGGTGGTc                     >  2:3068528/1‑101 (MQ=255)
                                                                                   tAAAACGGTTGGGATTT‑‑GCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCTTCACTGGGCGCGCGTGAAAGGGTTGTGGTGGTCGa                   <  2:3067875/101‑1 (MQ=255)
                                                                                   tAAAACGGTTGGGATTT‑‑GCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCTTCACTGGGCGCGCGTGAAAGGGTTGTGGTGGTCGa                   >  1:3068530/1‑101 (MQ=255)
                                                                                    aaaaCGGTTGGGATTT‑‑GCCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCTTCACTGGGCGCGCGTGAAAGGGTTGTGGTGGTCGa                   <  1:3068282/101‑1 (MQ=255)
                                                                                    aaaaCGGTTGGGATTT‑‑GCCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGAGCCAGTGCTTCACTGGGCGCGCGTGAAAGGGTTGTGGTGGTCGa                   >  1:3068534/1‑101 (MQ=255)
                                                                                    aaaaCGGTTGGGATTT‑‑GCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCTTCACTGGGCGCGCGTGAAAGGGTTGTGGTGGTCGAt                  >  2:3068532/1‑101 (MQ=255)
                                                                                       aCGGTTGGGATTT‑‑GCCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCTTCACTGGGCGCGCGTGAAAGGGTTGTGGTGGTCGAtgt                >  2:3068537/1‑101 (MQ=255)
                                                                                          gTTGGGATTTGCG‑‑CCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCTTCACTGGGCGCGCGTGAAAGGGTTGTGGTGGTCGAtgtgca             <  1:3068365/101‑3 (MQ=255)
                                                                                           ttGGGATTT‑‑GCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCTTCACTGGGCGCGCGTGAAAGGGTTGTGGTGGTCGATGTGGAAg           <  2:3068011/101‑1 (MQ=255)
                                                                                             gggATTT‑‑GCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCTTCACTGGGCGCGCGTGAAAGGGTTGTGGTGGTCGATGTGGAAGAt         >  1:3068541/1‑101 (MQ=255)
                                                                                             gggATTT‑‑GCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCTTCACTGGGCGCGCGTGAAAGGGTTGTGGTGGTCGATGTGGAAGAt         >  2:3068540/1‑101 (MQ=255)
                                                                                               gATTT‑‑GCCCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCTTCACTGGGCGCGCGTGAAAGGGTtg                             >  1:3068544/1‑79 (MQ=255)
                                                                                                       |                                                                                                    
GCCGGTTTCTGCTGCGCCACTGCTGCAGGTGAGCGGCGCACTCATCGCCATTATTGCCCTGATCCTCGCTGCTGCCTGGCTGGTAAAACGGTTGGGATTT‑‑G‑‑CCCCTAAACGCACTGGCGTTAACGGTCTGAAAATTAGCGCCAGTGCTTCACTGGGCGCGCGTGAAAGGGTTGTGGTGGTCGATGTGGAAGATGCACGGC  >  NC_000913/2021536‑2021735

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: