Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,806,430 C→A 9.2% R213L (CGC→CTC)  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,806,4300CA 9.2% 203.0 / 4.1 68R213L (CGC→CTC) waaGUDP‑glucose:(heptosyl)lipopolysaccharide alpha‑1,3‑glucosyltransferase; lipopolysaccharide core biosynthesis protein; lipopolysaccharide glucosyltransferase I
Reads supporting (aligned to +/- strand):  ref base C (47/14);  new base A (2/5);  total (49/19)
Fisher's exact test for biased strand distribution p-value = 1.55e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGACCGGAGAAGAAATGCACATTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGG  >  NC_000913/3806286‑3806572
                                                                                                                                                |                                                                                                                                              
cGACCGGAGAAGAAATGCACATTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATc                                                                                                                                            >  4:237743/1‑149 (MQ=255)
      gagaAGAAATGCACATTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGc                                                                                                                                                >  8:196250/1‑139 (MQ=255)
      gagaAGAAATGCACATTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGc                                                                                                                                                >  8:191589/1‑139 (MQ=255)
      gagaAGAAATGCACATTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGc                                                                                                                                                <  7:196250/139‑1 (MQ=255)
      gagaAGAAATGCACATTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGc                                                                                                                                                <  7:191589/139‑1 (MQ=255)
          agaaATGCACATTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTa                                                                                                                                  <  7:58061/149‑1 (MQ=255)
          agaaATGCACATTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGACTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTa                                                                                                                                  >  2:195603/1‑149 (MQ=255)
             aaTGCAGATTGCTCCGCACTCCGACTTTTTCTGCCTGCGCTTCAAAATTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACgg                                                                                                                               <  3:249389/149‑1 (MQ=255)
                  acaTTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCaa                                                                                                                           >  5:228337/1‑148 (MQ=255)
                  acaTAGCTCCGCACGCCGAGTTTTTCTGCAAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTTTTTTTTTTTTTTGATAGGTTGAATTGAGAGATCTACACCTTTACGGCCaaa                                                                                                                          <  2:258322/149‑1 (MQ=255)
                    aTTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGTCCAAAAt                                                                                                                        >  5:197457/1‑149 (MQ=255)
                    aTTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAAGTGAGCGATCTACACCTTTACGGCCAAAAt                                                                                                                        >  3:201534/1‑149 (MQ=255)
                     ttGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTATCGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCaaaa                                                                                                                         <  1:181015/147‑1 (MQ=255)
                     ttGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTATCGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCaaaa                                                                                                                         >  2:181015/1‑147 (MQ=255)
                         tCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACTGCAAAAATCTGAt                                                                                                                   >  1:169219/1‑149 (MQ=255)
                          ccGCACGCCGAGTCTTTCTGACAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGTTAACGATGCCAAAGCTTCAATTGACCGATCTACACCTTTACGGCCAAAATCTGATc                                                                                                                  <  5:273661/149‑1 (MQ=255)
                           cGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATTATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATc                                                                                                                  >  6:204423/1‑148 (MQ=255)
                            gCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATcc                                                                                                                 >  8:173496/1‑148 (MQ=255)
                            gCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCa                                                                                                                <  3:237743/149‑1 (MQ=255)
                              aCGCCGATTTCTCCTGCCAGCGCAGCATATTTTCGCGGCTTATCCTGACCAACAACTAATAAAAGCGTTTTGTCACGTAATGATTTTTTTTTTTTTTTGATAGGTGGAATTGAGAGATCTACACCTTTAGGGCCAAAATCTGATCCAAc                                                                                                              <  6:19414/149‑1 (MQ=255)
                                              ccAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTCCACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGttt                                                                                              >  5:172802/1‑149 (MQ=255)
                                               cAGCGCTTCAAATTTTCGCGACTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGtttt                                                                                             <  1:195603/149‑1 (MQ=255)
                                                        aaaTTTTGGCGGCTTATCCTGACCAACACCAAATAATTGCGTATTGTGACGTAATGAATCCGGTTACCATGCCAAAGCTTCTCTTGAGCGATCTAGACCTTTACGGCCAAAATCTGATCCAACCTGCAGtaataa                                                                                                  <  5:190232/135‑1 (MQ=255)
                                                        aaaTTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGAGACGAAATCATTCCGGTAACGATGCCAAATCTTCAATTGAGCGATCTAAAACTTTACGGCCAAAATCTGAAACAACCTGCAGtaataa                                                                                                  >  6:190232/1‑135 (MQ=255)
                                                          aTTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCtttt                                                                                  <  3:269840/149‑1 (MQ=255)
                                                           ttttCGCGGCTTATCCTGCCCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACCCCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGATGCTCTTTTa                                                                                 >  1:206771/1‑149 (MQ=255)
                                                                gcgGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCa                                                                            >  8:112591/1‑149 (MQ=255)
                                                                   gCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCAttt                                                                         <  7:173496/149‑1 (MQ=255)
                                                                          cTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGc                                                                  >  8:283169/1‑149 (MQ=255)
                                                                            gACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGa                                                                >  8:329067/1‑149 (MQ=255)
                                                                             aCCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGa                                                                >  7:315787/1‑148 (MQ=255)
                                                                              ccAACAACAAATAAAAGCGTATTGTGACGAAATGATTTTTTTTTTTTTTTGATAGGTTGAATAGAGAGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATa                                                              <  8:319570/149‑1 (MQ=255)
                                                                                 acaacaaATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAAt                                                           <  7:112591/149‑1 (MQ=255)
                                                                                  caacaaATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAAtt                                                          >  8:144868/1‑149 (MQ=255)
                                                                                   aacaaATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGATATCTACACATTTACGGCCAAAATCTGATCCACCCTGCAGTAATAAGTTTTGTTGCTATTTTATGCCATTTATCTGGCGATAAAttt                                                         >  3:239068/1‑149 (MQ=255)
                                                                                    acaaATAAAAGAGTATTGTGACGTAATGATTCAGGTAACGATGCCAAAGCTTCAATTGAGAGATATACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTATTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAAttt                                                         >  4:216220/1‑148 (MQ=255)
                                                                                              gCGGATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGCGCGCTCTACACCTTTACGGCCAAAATCTGATCCACCCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGttt                                              >  7:100985/1‑149 (MQ=255)
                                                                                                   ttGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGAt                                         >  2:235009/1‑149 (MQ=255)
                                                                                                     gtgACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGAttt                                       >  8:285180/1‑149 (MQ=255)
                                                                                                      tgACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTg                                      >  8:310516/1‑149 (MQ=255)
                                                                                                        aCGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTgg                                            >  2:300598/1‑141 (MQ=255)
                                                                                                        aCGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTgg                                            <  1:300598/141‑1 (MQ=255)
                                                                                                        aCGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCt                                    >  5:169213/1‑149 (MQ=255)
                                                                                                          gTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCa                                  >  1:94028/1‑149 (MQ=255)
                                                                                                           tAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCAc                                 >  6:220503/1‑149 (MQ=255)
                                                                                                            aaTGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACt                                >  1:251103/1‑149 (MQ=255)
                                                                                                            aaTGATTCCGGTAACCAAGCCAAAGCTACAATTGAGCGATATACACCATTACTGCCAAAATCACATCCAACCTCCAGAAATAACTTTTGTTGCTCTTTTATGCAATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACt                                >  2:189996/1‑149 (MQ=255)
                                                                                                               gATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTa                             >  2:58134/1‑149 (MQ=255)
                                                                                                                aTTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTAt                            >  8:296765/1‑149 (MQ=255)
                                                                                                                 ttCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTAtt                           >  5:309643/1‑149 (MQ=255)
                                                                                                                 ttCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTCATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGCCACATTGGCATTTGCTCACTGTAtt                           >  4:11933/1‑149 (MQ=255)
                                                                                                                          ttttttgaaaGCGTGAATTGAGAGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTcc                   <  6:228337/141‑1 (MQ=255)
                                                                                                                          cGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCg                  >  5:268256/1‑149 (MQ=255)
                                                                                                                           gATGCCAAAGCTTCAATTGAGCGATCTACACCTTTGCGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCg                  >  3:290817/1‑148 (MQ=255)
                                                                                                                           gATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCg                  >  8:128403/1‑148 (MQ=255)
                                                                                                                              ttgaaaGCTTGAATTGAGAGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGAt               <  4:201534/145‑1 (MQ=255)
                                                                                                                              gCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGAt               >  5:57704/1‑148 (MQ=255)
                                                                                                                                 aaaGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGttt                                                                                              <  7:114773/66‑1 (MQ=255)
                                                                                                                                 aaaGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGttt                                                                                              >  8:114773/1‑66 (MQ=255)
                                                                                                                                 aaaGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTTTTGCTATTTTATGAAATTTTTCTGGCGATAAATTTCACGGCTGTTTGGTATTTGCTCACTGTATTTTCTGTCCGGATaaa            >  6:63277/1‑148 (MQ=255)
                                                                                                                                 aaaGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATaaa            >  2:293705/1‑148 (MQ=255)
                                                                                                                                 aaaGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATaaa            >  1:330702/1‑148 (MQ=255)
                                                                                                                                   aGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATAc         >  7:133359/1‑149 (MQ=255)
                                                                                                                                      ttCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCgg      >  1:118271/1‑149 (MQ=255)
                                                                                                                                      ttCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCgg      >  4:26131/1‑149 (MQ=255)
                                                                                                                                        cAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCgg      <  7:283169/147‑1 (MQ=255)
                                                                                                                                        cAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCgg      <  1:58134/147‑1 (MQ=255)
                                                                                                                                         aaTTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGggag   >  8:189452/1‑149 (MQ=255)
                                                                                                                                         aaTTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGggag   >  1:38565/1‑149 (MQ=255)
                                                                                                                                          aTTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGggag   >  1:289619/1‑148 (MQ=255)
                                                                                                                                          aTTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTATGTACGGATAAATACCGggag   >  7:269831/1‑148 (MQ=255)
                                                                                                                                             gAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGggagg  >  5:126204/1‑146 (MQ=255)
                                                                                                                                                |                                                                                                                                              
CGACCGGAGAAGAAATGCACATTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGG  >  NC_000913/3806286‑3806572

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