Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,807,026 A→G 100% G14G (GGT→GGC waaG ← UDP‑glucose:(heptosyl)lipopolysaccharide alpha‑1,3‑glucosyltransferase; lipopolysaccharide core biosynthesis protein; lipopolysaccharide glucosyltransferase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,807,0260AG97.7% 590.7 / ‑0.9 178G14G (GGT→GGCwaaGUDP‑glucose:(heptosyl)lipopolysaccharide alpha‑1,3‑glucosyltransferase; lipopolysaccharide core biosynthesis protein; lipopolysaccharide glucosyltransferase I
Reads supporting (aligned to +/- strand):  ref base A (2/3);  new base G (113/60);  total (115/63)
Fisher's exact test for biased strand distribution p-value = 3.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 181 total aligned reads displayed.

TCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCACCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACAGCTTTTCCGTAGCGG  >  NC_000913/3806929‑3807125
                                                                                                 |                                                                                                   
tCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCccgcc                                                                                                   >  1:1562724/1‑100 (MQ=255)
   ggCAATCGCCTTCCCACGACTTTGTATATACCCGAACATGGTGGCCCCGTGCTGCAACTGTTGATGCAATGCACATAAAGACACGTTGAAGCCCGCCaaa                                                                                                >  1:3735900/1‑100 (MQ=255)
   ggCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTTCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTTAAGCCCGCCaaa                                                                                                >  1:3750789/1‑100 (MQ=255)
   ggCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCaaa                                                                                                >  2:2979073/1‑100 (MQ=255)
   ggCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCaaa                                                                                                >  2:3645066/1‑100 (MQ=255)
   gaaaaTCGCCTTCCCACAACTGTAAATATACCCGAAAATAATAGCCTAGAACAACAACTGTTGATGCAATAAGCAAAAAGTCAAGTTGAAAACCACTaaa                                                                                                >  2:2161875/4‑100 (MQ=255)
    tccgatctcTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAAt                                                                                               <  1:2904969/92‑1 (MQ=255)
    gCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAAt                                                                                               >  1:3358671/1‑100 (MQ=255)
     cAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAAt                                                                                               <  1:2141514/99‑1 (MQ=255)
     cAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATg                                                                                              >  1:2818407/1‑100 (MQ=255)
      aaTCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGTCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATgg                                                                                             <  1:2489615/100‑1 (MQ=255)
      aaTCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATgg                                                                                             >  2:1277365/1‑100 (MQ=255)
          gCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATACGCAAAAAGTCACGTTGAAGCCCGCCAAATGGaaaa                                                                                         >  1:3219407/1‑100 (MQ=255)
           ccTTCCCACGACTGTGTATATAACCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAat                                                                                        <  1:2084138/100‑1 (MQ=255)
            cTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAAtagatcgga                                                                                       >  2:2904943/1‑92 (MQ=255)
             ttCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAatat                                                                                      >  2:1153901/1‑100 (MQ=255)
               cccACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATAttt                                                                                    >  1:2906594/1‑100 (MQ=255)
               cccACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATAttt                                                                                    >  2:1626326/1‑100 (MQ=255)
               cccACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATAttt                                                                                    >  2:759541/1‑100 (MQ=255)
               cccACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATAttt                                                                                    >  1:1337631/1‑100 (MQ=255)
                ccACTACTGTGTATATACCCGAACATGGTGGCCCCTTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTta                                                                                   <  1:319584/100‑1 (MQ=255)
                ccACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTta                                                                                   <  1:1802789/100‑1 (MQ=255)
                ccACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTta                                                                                   <  1:3898267/100‑1 (MQ=255)
                 cACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGTCAACTGTTGATACAATGCGCATAAAATCACGTTGAAGCCCGCCAAATGGAAAATAttaa                                                                                   >  1:1152618/1‑97 (MQ=255)
                  aCGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTtata                                                                                 <  2:1981033/100‑1 (MQ=255)
                       tgtgTATGTACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCTCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAAc                                                                            <  2:469430/100‑1 (MQ=255)
                       tgtgTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTTATGCAATGCGCATAAAGTCACGTTGAAGCCCACCAAATGGAAAATATTTATATAAAc                                                                            <  2:599506/100‑1 (MQ=255)
                        gtgtATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACa                                                                           >  1:2267966/1‑100 (MQ=255)
                        gtgtATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCACGCCAAATGGAAAATATTTATATAAACa                                                                           >  1:657302/1‑100 (MQ=255)
                         tgtATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACaa                                                                          >  1:3604235/1‑100 (MQ=255)
                         tgtATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACaa                                                                          >  2:1927593/1‑100 (MQ=255)
                           tttataCCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACaaaa                                                                        >  1:2651702/3‑100 (MQ=255)
                           tatataCCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACaaaa                                                                        >  1:2203732/1‑100 (MQ=255)
                             tataCCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACaaaa                                                                        >  1:255790/1‑98 (MQ=255)
                             tataCCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAAc                                                                       >  1:275556/1‑99 (MQ=255)
                             tataCCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACg                                                                      >  1:1715069/1‑100 (MQ=255)
                               taCCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGcc                                                                    >  1:711560/1‑100 (MQ=255)
                               taCCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGcc                                                                    >  1:1088434/1‑100 (MQ=255)
                                 cccGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCAc                                                                  >  1:3475774/1‑100 (MQ=255)
                                 cccGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCAc                                                                  >  1:565001/1‑100 (MQ=255)
                                 cccGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCAc                                                                  >  2:2867609/1‑100 (MQ=255)
                                 cccGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCAc                                                                  >  2:3179198/1‑100 (MQ=255)
                                   cGAACATGGTGTCTCCGTGCGGCAACTGTTTTTGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACg                                                                 <  1:577307/99‑1 (MQ=255)
                                    gAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGa                                                                >  1:1834817/1‑99 (MQ=255)
                                     aaCATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATc                                                              <  2:2216791/100‑1 (MQ=255)
                                      aCATGGTGGCAACGTGCGGCAAATGTTGATGCAATGCGCATAAAGTCTCGTTGAAGCCCGCCAAAGGGAAACTATTTATATAAACAAAACGCCAc                                                                  >  2:740331/1‑95 (MQ=255)
                                         tggtggCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATaa                                                          >  1:2671999/1‑100 (MQ=255)
                                         tggtggCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAACCCCGCCAAATGGAAAAAATATATATAAACAAAACGCCACAATCaaaa                                                          >  1:1580789/1‑100 (MQ=255)
                                          ggtggCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAAt                                                         <  1:3791544/100‑1 (MQ=255)
                                          ggtggCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAAt                                                         <  2:350307/100‑1 (MQ=255)
                                           gtggCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAAtt                                                        <  1:889207/100‑1 (MQ=255)
                                            ggcccagatGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATtg                                                       <  2:203851/92‑1 (MQ=255)
                                             ggCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATtgt                                                      <  2:3297348/100‑1 (MQ=255)
                                               ccccGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGc                                                    >  2:2942861/1‑100 (MQ=255)
                                                cccGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCa                                                   <  2:2465363/100‑1 (MQ=255)
                                                    aacGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCAACGATAATAAttataca                                                   >  2:227652/3‑91 (MQ=255)
                                                      cGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTg                                             >  2:3352705/1‑100 (MQ=255)
                                                       ggCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTg                                             <  2:1332975/99‑1 (MQ=255)
                                                       ggCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGt                                            >  1:2134795/1‑100 (MQ=255)
                                                        gCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTc                                           <  2:3810345/100‑1 (MQ=255)
                                                         cAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCa                                          <  1:3514316/100‑1 (MQ=255)
                                                         cAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCAAGATCATAATTGTACAttaatatca                                          >  1:2417872/1‑93 (MQ=255)
                                                         cAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTCATATAAACAACACGCCACGATCATAATTATGCATTCCTGtta                                          <  1:2643665/100‑3 (MQ=255)
                                                           aCTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGc                                        <  1:402664/100‑1 (MQ=255)
                                                           aCTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGc                                        <  1:2591182/100‑1 (MQ=255)
                                                            gatctGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCt                                       <  1:61556/96‑1 (MQ=255)
                                                            cTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCt                                       >  2:1599379/1‑100 (MQ=255)
                                                             tGTTGATGCAATGCGCATAAAGTCACGTTGAAGCACGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGc                                        <  2:2671975/98‑1 (MQ=255)
                                                             tGTTGATGCAATGCGCATAAAGTAACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATAATAATTGTGCATTCCTGTCAGCTg                                      >  2:2632048/1‑100 (MQ=255)
                                                              gTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGa                                     >  1:2222554/1‑100 (MQ=255)
                                                               ttGATGCAATGCGCATAAAGTCACGTTGAAGCCCACCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGaa                                    <  1:81744/100‑1 (MQ=255)
                                                                 gATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCtagatc                                  >  2:61553/1‑95 (MQ=255)
                                                                 gATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAgg                                  >  2:1763812/1‑100 (MQ=255)
                                                                 gATGCAATGCGCATAAAGTCACGTTGAAGACCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAgg                                  <  2:849848/100‑1 (MQ=255)
                                                                  aTGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAggg                                 >  1:2929746/1‑100 (MQ=255)
                                                                  aTGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAggg                                 <  2:1000221/100‑1 (MQ=255)
                                                                   tGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAgggg                                >  1:149798/1‑100 (MQ=255)
                                                                   tGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAgggg                                >  1:1973052/1‑100 (MQ=255)
                                                                    gCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGc                               >  2:3911380/1‑100 (MQ=255)
                                                                      aaTGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCAt                             >  1:2201804/1‑100 (MQ=255)
                                                                       aTGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCAtc                            >  2:414617/1‑100 (MQ=255)
                                                                           gcATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCAtcttct                        >  1:1063310/1‑100 (MQ=255)
                                                                            cATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTg                       <  2:1333349/100‑1 (MQ=255)
                                                                             aTAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTAAGCTGAAGGGGCATCTTCTgg                      >  1:521323/1‑100 (MQ=255)
                                                                              tAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGc                     >  1:2449070/1‑100 (MQ=255)
                                                                              tAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGACATCTTCTGGc                     >  2:1464644/1‑100 (MQ=255)
                                                                               aaaGTCACGTTGAAGCCTGCAAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCa                    <  2:2534046/100‑1 (MQ=255)
                                                                                  gTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACa                 >  1:2584735/1‑100 (MQ=255)
                                                                                  gTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACa                 >  2:2355610/1‑100 (MQ=255)
                                                                                  gTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAACTGAAGGGGAATCTTCTGGCAACa                 >  2:2777109/1‑100 (MQ=255)
                                                                                  gTCACGTTGAAGCCCGCCAAATGAAAAACATTTATATAAACAAACCCCCACGATCATAATTGTGCTTTCCTGTCATCTGAAGGTGCATCTTCTGGCAAc                  <  2:2702771/99‑1 (MQ=255)
                                                                                  gTCACGTTGAAGCCCGCCAAATAGAAAATATTTATACAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACa                 <  2:3642648/100‑1 (MQ=255)
                                                                                    cACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACAg                >  2:1921868/1‑99 (MQ=255)
                                                                                    cACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACAGc               <  1:2026288/100‑1 (MQ=255)
                                                                                      cGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGCTCATAATTGTGCATTCCTGTCAGATGAAGGGGCATCTTCTGGCAACAGCtt             >  1:229992/1‑100 (MQ=255)
                                                                                      cGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAAAGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACAGCtt             >  2:2998834/1‑100 (MQ=255)
                                                                                       gTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACAGCttt            >  2:547347/1‑100 (MQ=255)
                                                                                           aaGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACAGCTTTTCCg        <  2:2094557/100‑1 (MQ=255)
                                                                                           aaGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACAGCTTTTCCg        <  2:653351/100‑1 (MQ=255)
                                                                                           aaGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACAGCTTTTCCg        >  1:2968783/1‑100 (MQ=255)
                                                                                                 |                                                                                                   
TCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCACCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACAGCTTTTCCGTAGCGG  >  NC_000913/3806929‑3807125

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