Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,806,430 C→A 9.1% 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.1% 193.6 / 3.5 65R213L (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 (32/26);  new base A (2/5);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 2.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGACCGGAGAAGAAATGCACATTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGG  >  NC_000913/3806286‑3806572
                                                                                                                                                |                                                                                                                                              
cGACCGGAGAAGAAATGCACATTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGa                                                                                                                                              <  1:135898/147‑1 (MQ=255)
cGACCGGAGAAGAAATGCACATTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGa                                                                                                                                              <  1:14945/147‑1 (MQ=255)
cGACCGGAGAAGAAATGCACATTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGa                                                                                                                                              <  1:174079/147‑1 (MQ=255)
cGACCGGAGAAGAAATGCACATTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGAt                                                                                                                                             <  7:124018/148‑1 (MQ=255)
       agaagaAATGCCCATTGCTCCGCCCGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGGTTATCCAGCCCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTTAGCGATCTACAcc                                                                                                                                      >  7:76955/1‑148 (MQ=255)
        gaagaaATGCACATTGCTCCGCACGCAGAGTTTTTCTGCCAGCGCTTAAAATTTTAGCTGCTTATACTGAACAACAACAAATAAAAGCGTATTGTGAAGTAATGATTCCGCTAACGATGACAAAGCTTCAATTGAGCGATCTACACCt                                                                                                                                     >  6:83147/1‑148 (MQ=255)
          agaaATGCACATTGCTCCTCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCttt                                                                                                                                   <  7:294035/148‑1 (MQ=255)
          agaaATGCACATTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACATCAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCtt                                                                                                                                    >  8:94493/1‑147 (MQ=255)
                     ttGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAAt                                                                                                                        >  2:154278/1‑148 (MQ=255)
                      tGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATc                                                                                                                       >  8:177624/1‑148 (MQ=255)
                          ccGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCATCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGAt                                                                                                                   >  8:302886/1‑148 (MQ=255)
                          ccGCACGCCGACTTTTTCTGCCAGCGCTTCAAATTTTCGCGTCTTATCCTGACAAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATAGAGCCATCTACACCTTTACGCCCAAAAACTGAt                                                                                                                   >  4:197310/1‑148 (MQ=255)
                             cACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCTGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGTTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCttt                                                                                                                                   <  5:78255/129‑1 (MQ=255)
                             cACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACACCAAATAAAAGCGTATTGGGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCa                                                                                                                <  3:144562/148‑1 (MQ=255)
                             cACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTCAGCGATCTACACCttt                                                                                                                                   >  6:78255/1‑129 (MQ=255)
                              aCGCCGAGTTTTTCTGCCAGAGCTGCATAATTTCGCTGCTTATCCAGCCCATCAACAAATAAAAGCGTATTGTGACGTAATGATTTTTTTTTTTTATTGATAGGGTGAATTGAGAGATCTACACCTTTACGGCCAAAATCTGATCCaa                                                                                                               <  6:99781/148‑1 (MQ=255)
                              aCGCCGAGTTGTTCTGCCACCGCTTCAAATTTTCGCGGCTTATCCTCACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTTTTTTTTTTTTTTGATAGGTTGAATTGAGAGATCTACACCTTTACGGCCAAAATCTGATCCaa                                                                                                               <  8:259604/148‑1 (MQ=255)
                               cGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAAc                                                                                                              <  7:177624/148‑1 (MQ=255)
                                   gattttttCTGCCAGCGCTTCATATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGc                                                                                                          <  3:197310/145‑1 (MQ=255)
                                       ttttCTGCCAGCGCTACAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGAGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCTATCAACACCATTACGGCCAAAATCTGATCCAAACTGCAGt                                                                                                       >  2:155144/1‑147 (MQ=255)
                                              ccAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGt                                                                                                <  3:105295/147‑1 (MQ=255)
                                                         aaTTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCtt                                                                                    >  4:246277/1‑148 (MQ=255)
                                                          aTTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCttt                                                                                   <  1:99406/148‑1 (MQ=255)
                                                           ttttCGAGGCTTATCCTGACCAACACAAAATAAAAGCCTATTGTGACTTAATGATTTAGTTTTTTAATTGATAGGTTTAATTGAGAGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCtttt                                                                                  <  8:76955/148‑1 (MQ=255)
                                                                      tATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAATGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCAtttt                                                                        >  2:160486/1‑147 (MQ=255)
                                                                       aTCCTGACCTACAACAAATAAAAGCGTATTGTGACGTAATGATTACGGTAACGATGCCAAAGCTTCAATTGAGAGATCTACAACTTTACGGACAAAATCTGATCAAACCTTCAGTAATAAGTTTTGTATCTATTTTATGACATTTTTc                                                                      >  3:13157/1‑148 (MQ=255)
                                                                       aTCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCTGTAACGATGCAAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGAAATAAGTTTTGTTGCTCTTTTATGCCAttttt                                                                       <  1:155144/147‑1 (MQ=255)
                                                                           tGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGc                                                                  >  1:290364/1‑148 (MQ=255)
                                                                             aCCAACAACAAATAAAAGCGTATTGTGACGTAATCATTCCGGTAACGAAGCCAAAGCTTCAATTCAGCGATCTACACCTTTACGGCCAAAATATGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCg                                                                 >  2:42694/1‑147 (MQ=255)
                                                                                aacatcaaATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATaa                                                             <  7:94493/148‑1 (MQ=255)
                                                                                aacaacaaATAAAAGCTTATTGTGACGTAATGATTCATGTATAGATGACAAAGCTTCAATTGAGCTATATACACCTTTAATGACAAAATCTGATCCAACCAGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATaa                                                             <  5:83147/148‑1 (MQ=255)
                                                                                    acaaATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGATAGATCTACACCTTTAAGGCAAAAATCTGATACAACCTGCATTAATAATTTTTGTTGCTCTTTTATGCCATTTATCTGGCGATAAAtt                                                          >  5:260330/1‑147 (MQ=255)
                                                                                      aaaTAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCa                                                       <  3:246277/148‑1 (MQ=255)
                                                                                                        aCGTAATGATTCCGGTAATGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGc                                     <  1:160486/148‑1 (MQ=255)
                                                                                                        aCGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGc                                     >  4:137034/1‑148 (MQ=255)
                                                                                                        aCGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGc                                     >  2:257539/1‑148 (MQ=255)
                                                                                                          gTAATGATTACGTTAATGATTTCATAGGTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTAGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCAc                                                      <  8:236053/129‑1 (MQ=255)
                                                                                                           tAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTAGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACg                                                     >  7:236053/1‑129 (MQ=255)
                                                                                                            aaTGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGt                                                <  1:311138/133‑1 (MQ=255)
                                                                                                            aaTGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGt                                                >  2:311138/1‑133 (MQ=255)
                                                                                                            aaTGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCAc                                 >  4:194857/1‑148 (MQ=255)
                                                                                                            aaTGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCAc                                 >  1:244541/1‑148 (MQ=255)
                                                                                                            aaTGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCAc                                 >  8:319849/1‑148 (MQ=255)
                                                                                                                 ttCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTAt                            <  1:154278/148‑1 (MQ=255)
                                                                                                                     ggTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTg                                                                    >  8:42042/1‑104 (MQ=255)
                                                                                                                     ggTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTg                                                                    <  7:42042/104‑1 (MQ=255)
                                                                                                                      gTAACGATGCCAAAGCTTCAATTGAGCGATCTACCCCTTTACGGCCATAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGCATTTGCTCACTGTATTTTCt                       >  5:161464/1‑148 (MQ=255)
                                                                                                                      gTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCtttt                                                                                  >  5:9101/1‑89 (MQ=255)
                                                                                                                        aaCAATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGt                     >  6:186586/1‑148 (MQ=255)
                                                                                                                         aCGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGt                     <  3:194857/147‑1 (MQ=255)
                                                                                                                            aTGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGttt                                              <  1:105103/119‑1 (MQ=255)
                                                                                                                            aTGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGttt                                              >  4:278367/1‑119 (MQ=255)
                                                                                                                            aTGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGttt                                              <  3:278367/119‑1 (MQ=255)
                                                                                                                            aTGCCAAAGCTGCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGttt                                              >  2:105103/1‑119 (MQ=255)
                                                                                                                             ttccAAAGTTTCAATTGAGAGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCgg                 <  6:161464/145‑1 (MQ=255)
                                                                                                                               ccAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCCACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTACTCACTGTATTTTCTGTCCGGAt               <  7:302886/147‑1 (MQ=255)
                                                                                                                                  aaGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATaaa            >  8:291863/1‑147 (MQ=255)
                                                                                                                                   aGCTTCAATTGAGCGATCTACACCTTTACTGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATa          >  2:2047/1‑148 (MQ=255)
                                                                                                                                   aGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATa          >  6:105339/1‑148 (MQ=255)
                                                                                                                                   aGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATACAAAATGCAGTAATAAGTTTTTTTGCTATTTTATTCCATTTTTCTGGCGATAAATATCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATa          >  3:40610/1‑148 (MQ=255)
                                                                                                                                    gCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATAc         >  8:93499/1‑148 (MQ=255)
                                                                                                                                    gCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATAc         <  5:186586/148‑1 (MQ=255)
                                                                                                                                    cgggtatttGAGAGATATACACCTTTACTGAAATAATATGATCCAAACTGCAGTAATAAGTTTTGTTGATCTTTTATGCAATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATa          <  6:21666/140‑1 (MQ=255)
                                                                                                                                         aaTTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCggg     <  1:2047/147‑1 (MQ=255)
                                                                                                                                            tGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTAACGGCTGTTTGGGATTTGCTCACTGTATTTTCTCTCCGGATAAAAACCGggagg  >  8:118631/1‑147 (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: ‑