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.9% 632.3 / ‑3.5 177G14G (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/0);  new base G (81/94);  total (83/94)
Fisher's exact test for biased strand distribution p-value = 2.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

Only 100 of 181 total aligned reads displayed.

TCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCACCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACAGCTTTTCCGTAGC  >  NC_000913/3806929‑3807123
                                                                                                 |                                                                                                 
tCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCccgcc                                                                                                 <  2:1441570/100‑1 (MQ=255)
tCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTCTTGATGCAATGCGCATAAAGTCACGTTGAAGCccgcc                                                                                                 <  2:2232536/100‑1 (MQ=255)
 cGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCccgcca                                                                                                >  2:3350737/1‑100 (MQ=255)
  gggCAATCGCCTTCCCACGACTGTGTATATACCCGCACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATCCGCATAACGTCACGTTTAAGCCCGCCaa                                                                                               >  1:1864046/1‑100 (MQ=255)
   ggCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCaaa                                                                                              >  1:21085/1‑100 (MQ=255)
    gCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAAt                                                                                             >  2:3394032/1‑100 (MQ=255)
     cAATCGCCTTCCCACGACTGTGTATATACCCGAACATTGTGGCACCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATg                                                                                            <  2:2526854/100‑1 (MQ=255)
     cAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATg                                                                                            <  2:2793059/100‑1 (MQ=255)
        tCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGaa                                                                                         <  2:3686470/100‑1 (MQ=255)
         cGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGaaa                                                                                        <  1:3773860/100‑1 (MQ=255)
                 cACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTtat                                                                                >  2:2945750/1‑100 (MQ=255)
                  aCGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTtata                                                                               >  1:3525477/1‑100 (MQ=255)
                   cGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTtatat                                                                              >  2:2162591/1‑100 (MQ=255)
                   cGACTGTATATATACCCGAACATGGTGACCCCGTGCGGCAACTGTTGATGCAATGCGCATATAGTCACGTTGAAGCCCGCCAAATGGAAAATATTtatat                                                                              <  2:3427946/100‑1 (MQ=255)
                     aCTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATaa                                                                            >  2:2064036/1‑100 (MQ=255)
                     aCTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATaa                                                                            <  1:3044768/100‑1 (MQ=255)
                       tgtgTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAAc                                                                          >  1:314999/1‑100 (MQ=255)
                         tgtATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACaa                                                                        >  1:2654409/1‑100 (MQ=255)
                         tgtATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATAGTTATATAAACaa                                                                        >  1:778479/1‑100 (MQ=255)
                            atataCCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATACAATGCTCATAAATTCACGTTGAAGCCAGCCAAATGGAAAATATTTATATAAACAAAAc                                                                     <  2:392746/100‑1 (MQ=255)
                             tataCCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACg                                                                    >  2:3083261/1‑100 (MQ=255)
                             tataCCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACg                                                                    <  1:313799/100‑1 (MQ=255)
                             tataCCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACg                                                                    >  2:1528192/1‑100 (MQ=255)
                               taCCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGcc                                                                  >  1:3185365/1‑100 (MQ=255)
                               taCCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGcc                                                                  >  1:3811874/1‑100 (MQ=255)
                               taCCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGcc                                                                  <  1:3419454/100‑1 (MQ=255)
                               taCCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGcc                                                                  <  2:1230506/100‑1 (MQ=255)
                                aCCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATTTAAACAAAACGCCa                                                                 >  2:978332/1‑100 (MQ=255)
                                aCCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCACATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCa                                                                 >  1:3277659/1‑100 (MQ=255)
                                 cccGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCAc                                                                >  1:1081471/1‑100 (MQ=255)
                                 cccGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCAc                                                                >  2:26279/1‑100 (MQ=255)
                                 cccGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCAc                                                                >  2:3686808/1‑100 (MQ=255)
                                 cccGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCAc                                                                >  1:1898450/1‑100 (MQ=255)
                                 cccGAACATGGTGGCCCCGTGCGGCAACTGTCGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCAc                                                                <  2:2815274/100‑1 (MQ=255)
                                  ccGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACg                                                               <  2:3525439/100‑1 (MQ=255)
                                   cGAACATGGTTGCCCCGTGCGGCAACTGTTGATTCAATGCGAATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGa                                                              <  2:670671/100‑1 (MQ=255)
                                     aaCATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATc                                                            <  2:1218140/100‑1 (MQ=255)
                                      aCATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCagtcatt                                                           >  2:2481847/1‑94 (MQ=255)
                                      aCATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCa                                                           <  2:1424631/100‑1 (MQ=255)
                                          ggtggCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAAt                                                       >  2:2638876/1‑100 (MQ=255)
                                           gtggCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAAtt                                                      <  1:2446234/100‑1 (MQ=255)
                                            tggCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATtg                                                     <  1:679324/100‑1 (MQ=255)
                                            tggCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATtg                                                     <  1:718251/100‑1 (MQ=255)
                                              gCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATtgtg                                                   <  2:1383907/100‑1 (MQ=255)
                                              gCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATtgtg                                                   <  2:1590917/100‑1 (MQ=255)
                                              gCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATtgtg                                                   >  1:2085227/1‑100 (MQ=255)
                                              gCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATtgtg                                                   >  1:3641419/1‑100 (MQ=255)
                                              gCCCCGTGCGGCAACTGCTGATGCAATGCACATAAAGACACGTTTAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACACCACAATCATAATtg                                                     <  2:1307081/98‑1 (MQ=255)
                                               ccccGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGc                                                  <  2:4039930/100‑1 (MQ=255)
                                               ccccGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGc                                                  <  2:1930075/100‑1 (MQ=255)
                                               ccccGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGc                                                  <  1:123960/100‑1 (MQ=255)
                                                   gTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTc                                              >  2:3044944/1‑100 (MQ=255)
                                                    tGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTcc                                             <  1:1943477/100‑1 (MQ=255)
                                                    tGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTcc                                             <  1:3551514/100‑1 (MQ=255)
                                                     gCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCt                                            >  1:1993571/1‑100 (MQ=255)
                                                      cGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTg                                           >  2:3578237/1‑100 (MQ=255)
                                                       ggCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGt                                          >  2:1261374/1‑100 (MQ=255)
                                                       ggCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGt                                          >  2:1078682/1‑100 (MQ=255)
                                                        gCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTc                                         <  2:2412466/100‑1 (MQ=255)
                                                         cAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAACCGCCACGATCATAATTGTGCATTCCTGTc                                         <  2:3247979/99‑1 (MQ=255)
                                                         cAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCa                                        >  1:1366489/1‑100 (MQ=255)
                                                          aaCTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTc                                         <  2:2683752/98‑1 (MQ=255)
                                                           aCTGTTGATGCAATGCGGATAAAGTCACTTTGAAGCCCGCCAAATTGAAAAGATTTAGATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGc                                      <  2:809514/100‑1 (MQ=255)
                                                           aCTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGc                                      <  1:512535/100‑1 (MQ=255)
                                                           aCTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGc                                      <  1:546883/100‑1 (MQ=255)
                                                           aCTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGc                                      <  2:3874203/100‑1 (MQ=255)
                                                              gTTGATGAAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAACCGCCACGATCATAATTGTGCATTCCTGTCAGCTGa                                   <  2:3139667/100‑1 (MQ=255)
                                                                tGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCTGCTGAAg                                 <  1:2199633/100‑1 (MQ=255)
                                                                tGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAg                                 <  1:3525282/100‑1 (MQ=255)
                                                                tGATGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAg                                 >  1:1798419/1‑100 (MQ=255)
                                                                 taggCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAgg                                <  1:3265308/97‑1 (MQ=255)
                                                                  aTGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAggg                               <  2:314994/100‑1 (MQ=255)
                                                                   tGCAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAgggg                              <  2:186363/100‑1 (MQ=255)
                                                                     cAATGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATCGTGCATTCCTGTCAGCTGAAGGGGCa                            >  2:1993655/1‑100 (MQ=255)
                                                                       aTGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCAtc                          <  1:2162611/100‑1 (MQ=255)
                                                                        tGCGCATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTATGCATTCCTGTCAGCTGAAGGGGCAtct                         >  2:176613/1‑100 (MQ=255)
                                                                         gcgcATAAAGTCACGTTGAAGCCCACCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCAtatt                        >  1:2463728/1‑100 (MQ=255)
                                                                          cgcATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCAtcttc                       >  2:4177972/1‑100 (MQ=255)
                                                                            caaaaaGTCACATTAAAGCCCGCCAAATAGACAATATTTATATAAACCAAACGCCACGATCATAATTTTGCATACCTCTCAGCTAAAGGGGCAtcttct                      <  2:3630563/96‑1 (MQ=255)
                                                                            cATAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTg                     >  2:1839157/1‑100 (MQ=255)
                                                                              tAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGc                   <  1:3241518/100‑1 (MQ=255)
                                                                              tAAAGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGc                   <  1:2990233/100‑1 (MQ=255)
                                                                               aaaGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCa                  <  1:1319146/100‑1 (MQ=255)
                                                                               aaaGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCa                  <  2:1371563/100‑1 (MQ=255)
                                                                                 aGTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAAc                <  2:2023790/100‑1 (MQ=255)
                                                                                  gacACGTTGAAGCACGCCAAATGGAAAATATGTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACa               >  1:2815705/3‑100 (MQ=255)
                                                                                  gTCACGTTGAAGCCCGCCAAATGGAAAATATTTATGTAAACAAACCGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAAc                <  1:1474823/99‑1 (MQ=255)
                                                                                  gTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACa               >  1:306453/1‑100 (MQ=255)
                                                                                  gTCACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACa               <  1:3847656/100‑1 (MQ=255)
                                                                                  gTCAAGTTGAAGCCCGCAAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACa               >  2:3061540/1‑100 (MQ=255)
                                                                                    cACGTTGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACAGc             <  2:3811832/100‑1 (MQ=255)
                                                                                        ttGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACAGCtttt         >  2:2540816/1‑100 (MQ=255)
                                                                                         tGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACAGCTTTTc        <  1:2540843/100‑1 (MQ=255)
                                                                                         tGAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACAGCTTTTc        <  1:4103031/100‑1 (MQ=255)
                                                                                          gAAGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACAGCttt          >  2:297792/1‑97 (MQ=255)
                                                                                           aaGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACAGCTTTTCCg      <  2:2358094/100‑1 (MQ=255)
                                                                                            aGCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACAGCTTTTCCGt     <  1:3394070/100‑1 (MQ=255)
                                                                                             gCCCGCCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACAGCTTTTCCGTa    >  1:4013151/1‑100 (MQ=255)
                                                                                              cccACCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACAGCTTTTCCGTAg   >  2:704411/1‑100 (MQ=255)
                                                                                               ccaccaAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACAGCTTTTCCGTAgc  >  1:228578/1‑100 (MQ=255)
                                                                                                 |                                                                                                 
TCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATAAAGTCACGTTGAAGCCCACCAAATGGAAAATATTTATATAAACAAAACGCCACGATCATAATTGTGCATTCCTGTCAGCTGAAGGGGCATCTTCTGGCAACAGCTTTTCCGTAGC  >  NC_000913/3806929‑3807123

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