Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,807,026 A→G 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,0260AG98.3% 417.5 / ‑2.6 119G14G (GGT→GGCwaaGUDP‑glucose:(heptosyl)lipopolysaccharide alpha‑1,3‑glucosyltransferase; lipopolysaccharide core biosynthesis protein; lipopolysaccharide glucosyltransferase I
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (77/40);  total (77/42)
Fisher's exact test for biased strand distribution p-value = 1.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

Only 100 of 122 total aligned reads displayed.

TCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCACCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACAGCTTTTCCGTA  >  NC_000913/3806929‑3807121
                                                                                                 |                                                                                               
tCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCccgcc                                                                                               >  2:776123/1‑100 (MQ=255)
 cGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCccgcca                                                                                              >  2:1932972/1‑100 (MQ=255)
 cGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCccgcca                                                                                              <  2:700670/100‑1 (MQ=255)
  gggCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCaa                                                                                             >  2:860454/1‑100 (MQ=255)
      aaTCGCCTTCCCACGACTGTGTATATACCCGAAGATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGTATAAAGTCACGTTGAAGCCCGCCAAATgg                                                                                         >  2:1851374/1‑100 (MQ=255)
      aaTCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATgg                                                                                         >  1:1993756/1‑100 (MQ=255)
       aTCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGa                                                                                        >  1:886132/1‑100 (MQ=255)
         cGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTCTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGaaa                                                                                      <  2:1486587/100‑1 (MQ=255)
          gCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGaaaa                                                                                     >  2:535721/1‑100 (MQ=255)
          gCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGaaaa                                                                                     >  2:1967693/1‑100 (MQ=255)
              tCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATAtt                                                                                 >  2:2564694/1‑100 (MQ=255)
               cccACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATAttt                                                                                >  1:1744497/1‑100 (MQ=255)
                ccACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTta                                                                               >  1:1546310/1‑100 (MQ=255)
                 cACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTtat                                                                              >  1:1757140/1‑100 (MQ=255)
                 cACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTtat                                                                              >  1:1407313/1‑100 (MQ=255)
                  aCGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACCGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTtata                                                                             <  2:1979506/100‑1 (MQ=255)
                          gtATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACaaa                                                                     >  1:2048605/1‑100 (MQ=255)
                           tatataCCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACaaaa                                                                    >  2:43176/1‑100 (MQ=255)
                            atataCCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAAc                                                                   >  2:2621830/1‑100 (MQ=255)
                            atataCCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAAc                                                                   >  2:2381086/1‑100 (MQ=255)
                            atataCCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAAc                                                                   >  1:54263/1‑100 (MQ=255)
                            atataCCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAAc                                                                   >  2:2633000/1‑100 (MQ=255)
                            atataCCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAAc                                                                   >  2:2966128/1‑100 (MQ=255)
                            atataCCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAAc                                                                   >  1:708893/1‑100 (MQ=255)
                             tataCCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACg                                                                  <  1:2915954/100‑1 (MQ=255)
                             tataCCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACg                                                                  >  2:411860/1‑100 (MQ=255)
                              ataCCCGAACATGGTGTCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATCTTTATATAAACAAAACGc                                                                 >  2:1702652/1‑100 (MQ=255)
                              ataCCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGc                                                                 >  2:668269/1‑100 (MQ=255)
                               taCCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGcc                                                                >  2:665057/1‑100 (MQ=255)
                               taCCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGcc                                                                <  2:250465/100‑1 (MQ=255)
                                aCCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGcc                                                                <  2:1288148/99‑1 (MQ=255)
                                 cccGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCAc                                                              >  1:1683168/1‑100 (MQ=255)
                                 cccGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCAc                                                              >  2:845086/1‑100 (MQ=255)
                                   cGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGa                                                            >  1:1317931/1‑100 (MQ=255)
                                   cGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGa                                                            >  1:2776480/1‑100 (MQ=255)
                                     aaCATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCACATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATc                                                          >  1:782455/1‑100 (MQ=255)
                                       cATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCAt                                                        <  1:1070751/100‑1 (MQ=255)
                                          ggtggCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAAt                                                     >  1:2702753/1‑100 (MQ=255)
                                           gtggCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAAtt                                                    >  1:34846/1‑100 (MQ=255)
                                              gCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATtgtg                                                 >  2:637740/1‑100 (MQ=255)
                                               ccccGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGc                                                >  2:110131/1‑100 (MQ=255)
                                                  cGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCACCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCAtt                                             <  1:2985511/100‑1 (MQ=255)
                                                  cGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCACCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCAtt                                             <  1:230698/100‑1 (MQ=255)
                                                    tGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTcc                                           <  1:2363583/100‑1 (MQ=255)
                                                    tGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTcc                                           >  1:1304003/1‑100 (MQ=255)
                                                      cGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTg                                         <  2:2863231/100‑1 (MQ=255)
                                                       ggCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGt                                        >  1:3057493/1‑100 (MQ=255)
                                                       ggCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGt                                        >  2:295227/1‑100 (MQ=255)
                                                        gCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTc                                       <  2:1744478/100‑1 (MQ=255)
                                                        gCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTc                                       >  2:1247080/1‑100 (MQ=255)
                                                        gCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTc                                       >  1:1041760/1‑100 (MQ=255)
                                                        gCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTc                                       <  1:2621853/100‑1 (MQ=255)
                                                        gCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTc                                       <  1:10122/100‑1 (MQ=255)
                                                         cAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCa                                      <  1:686354/100‑1 (MQ=255)
                                                         cAACTGTTGATACAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCa                                      >  2:1349862/1‑100 (MQ=255)
                                                          aaCTGTTGATGCAATGCGCATAAAGTAACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAg                                     <  1:1223753/100‑1 (MQ=255)
                                                           aCTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGc                                    <  2:1182154/100‑1 (MQ=255)
                                                           aCTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGc                                    >  1:161173/1‑100 (MQ=255)
                                                           aCTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGc                                    >  1:2159560/1‑100 (MQ=255)
                                                            cTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCt                                   >  2:1780548/1‑100 (MQ=255)
                                                              gTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGa                                 >  1:981926/1‑100 (MQ=255)
                                                              gTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGa                                 <  2:54263/100‑1 (MQ=255)
                                                              gTTGATGCAATGCGCATAAAGTAACGTTGAAGCCCGCCAAATGGAGAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTg                                  >  2:373020/1‑99 (MQ=255)
                                                                 gATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAgg                              >  2:2740899/1‑100 (MQ=255)
                                                                   tGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAgggg                            >  2:1833292/1‑100 (MQ=255)
                                                                   tGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAgggg                            >  1:638916/1‑100 (MQ=255)
                                                                   tGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAgggg                            <  2:2237484/100‑1 (MQ=255)
                                                                     cAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACACAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCa                          >  2:1413520/1‑100 (MQ=255)
                                                                     cAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCa                          <  1:1932994/100‑1 (MQ=255)
                                                                     cAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCa                          >  1:434774/1‑100 (MQ=255)
                                                                     cAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCa                          <  1:3083429/100‑1 (MQ=255)
                                                                          cgcATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCAtcttc                     >  2:1943601/1‑100 (MQ=255)
                                                                           gcATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCAtcttct                    >  2:2413841/1‑100 (MQ=255)
                                                                           gcATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCAtcttct                    >  1:1722381/1‑100 (MQ=255)
                                                                            cATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTg                   >  1:1762258/1‑100 (MQ=255)
                                                                            cATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTg                   <  1:2375067/100‑1 (MQ=255)
                                                                            cATAAAGTCACGATGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTg                   >  1:2026487/1‑100 (MQ=255)
                                                                              tAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGc                 <  2:161172/100‑1 (MQ=255)
                                                                              tAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGc                 <  1:637742/100‑1 (MQ=255)
                                                                              tAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGc                 <  2:2974663/100‑1 (MQ=255)
                                                                              naaaGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGc                 >  2:1572189/2‑100 (MQ=255)
                                                                               aaaGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCa                >  2:2218424/1‑100 (MQ=255)
                                                                                aaGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCaa               >  1:1708942/1‑100 (MQ=255)
                                                                                 aGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAAc              <  2:1735894/100‑1 (MQ=255)
                                                                                 aGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAAc              <  1:2381109/100‑1 (MQ=255)
                                                                                  gTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACa             >  1:2230900/1‑100 (MQ=255)
                                                                                  gTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACa             >  2:1648293/1‑100 (MQ=255)
                                                                                  gTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACa             >  1:1433683/1‑100 (MQ=255)
                                                                                   tCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAATCGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACAg            <  2:2861640/100‑1 (MQ=255)
                                                                                    cACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACAGc           >  1:2188511/1‑100 (MQ=255)
                                                                                    cACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACAGc           <  1:411862/100‑1 (MQ=255)
                                                                                     aCGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACAGCt          <  1:591962/100‑1 (MQ=255)
                                                                                      cGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACAGCtt         >  2:3065232/1‑100 (MQ=255)
                                                                                        ttGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACAGCtttt       <  2:1110282/100‑1 (MQ=255)
                                                                                          gAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACAGCTTTTcc     >  1:2583693/1‑100 (MQ=255)
                                                                                          gAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAACTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACAGCTTTTcc     <  1:482654/100‑1 (MQ=255)
                                                                                           aaGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACAGCTTTTCCg    <  1:341050/100‑1 (MQ=255)
                                                                                            aGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACAGCTTTTCCGt   >  1:2617986/1‑100 (MQ=255)
                                                                                             gCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACAGCTTTTCCGTa  >  2:1389026/1‑100 (MQ=255)
                                                                                             gCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACAGCTTTTCCGTa  >  1:1210821/1‑100 (MQ=255)
                                                                                                 |                                                                                               
TCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCACCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACAGCTTTTCCGTA  >  NC_000913/3806929‑3807121

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