Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,806,430 C→A 10.7% 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,4300CA10.7% 179.4 / 5.3 67R213L (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 (40/19);  new base A (2/6);  total (42/25)
Fisher's exact test for biased strand distribution p-value = 4.47e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGCGACCGGAGAAGAAATGCACATTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAG  >  NC_000913/3806283‑3806575
                                                                                                                                                   |                                                                                                                                                 
ttGCGACCGGAGAAGAAATGCACATTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTTTCCTGACCAACACCAAATAAAAGCGTATTGTCACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGc                                                                                                                                                   >  3:79939/1‑148 (MQ=255)
   cGACCGGAGAAGAAATGCACATTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGAt                                                                                                                                                <  6:178302/148‑1 (MQ=255)
   cGACCGGAGAAGAAATGCACATTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATc                                                                                                                                               <  6:240259/149‑1 (MQ=255)
       cGGAGAAGAAATGCACATTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTaca                                                                                                                                           <  1:221460/149‑1 (MQ=255)
                     acaTTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTTTTTTTTTTTATTGATAGGTTGAATTGAGAGATCTACACCTTTACGGCCaaa                                                                                                                             <  2:233879/149‑1 (MQ=255)
                     acaTTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCaa                                                                                                                              >  2:145910/1‑148 (MQ=255)
                       aTTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCaaaa                                                                                                                            >  4:265114/1‑148 (MQ=255)
                       aTTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAAt                                                                                                                           >  7:136701/1‑149 (MQ=255)
                       aTTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATAATGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTACTGTAACGATCCCAAAGCATCAATTGAGCGATCTACACCTTTACCGCCaaaa                                                                                                                            >  3:40068/1‑148 (MQ=255)
                        ttGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATc                                                                                                                          <  1:145910/149‑1 (MQ=255)
                         tGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCt                                                                                                                         >  8:368692/1‑149 (MQ=255)
                              cGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGAt                                                                                                                      <  7:388124/147‑1 (MQ=255)
                                cACGCCGACTTTTTCAGCCAGCGCTTCAACCTTTCGCGGCTTATCCTGCCCAACAACAAATAAAAGCGTCTTGTGACGTAATTATTTTTTTTTTTTTTTGATAGGTTGAATTGAGAGATCTACACCTTTACGGCCAAAATCTGATCCaa                                                                                                                  <  4:5366/149‑1 (MQ=255)
                                 actcagAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGCCTTATCCAGACAAACAAAAAATAAAAGCGAATTGTGACTTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAAc                                                                                                                 >  2:151422/6‑149 (MQ=255)
                                 aCGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATTATTCCGGTAACGATGACAAAGATTCAATTGAGATATCTACACCTTTACTGCCAAAATCTGATCCAAc                                                                                                                 >  1:388617/1‑149 (MQ=255)
                                     cGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAATGCGTATTGTGACGTAATGCTTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTg                                                                                                              >  3:75347/1‑148 (MQ=255)
                                     cGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCAGATCCAACCTg                                                                                                              >  2:349659/1‑148 (MQ=255)
                                            ttCTGCCAGCTCTTCAAATTTTCGCGGATTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTTTTTTTTTTTATTGATAGGGTGAATTGAGAGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGtaat                                                                                                       <  2:290511/148‑1 (MQ=255)
                                             tCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGtaata                                                                                                      <  7:161896/148‑1 (MQ=255)
                                                 ccAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGttt                                                                                                 >  3:243993/1‑149 (MQ=255)
                                                        ttCAAATTTTCGCGGCTTATCCGGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGc                                                                                           >  8:279053/1‑148 (MQ=255)
                                                         tCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGATTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTGTTGTTGctc                                                                                         >  2:271412/1‑149 (MQ=255)
                                                           aaaTTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCtt                                                                                       <  1:349659/149‑1 (MQ=255)
                                                                   gcgGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGACAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCa                                                                               >  3:68620/1‑149 (MQ=255)
                                                                     ggCTTATCCTGACCAACAACACATAACAGCGTATTGTGACGTAATGATTTTTTTTTTTTATTGATAGGGTGAATTGAGAGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCAtt                                                                             <  4:75347/149‑1 (MQ=255)
                                                                           tCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCAGTTTATTATTTGATAGCTTGAATTGAGAGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTg                                                                       <  4:79939/149‑1 (MQ=255)
                                                                            ccTGACCAACAACAAATAAAAGCTTATTGTGACGTAAAGATTCCGGTAACGATGCCAAAGCTTCAAATGAGAGATCTACACCTTTACTCCCAAAATATCATCCAACCTGAAGTAATATGATTTGTTGATATTTTATGCCATTTTTCTg                                                                       >  7:368570/1‑148 (MQ=255)
                                                                            ccTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTTACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTTCTCTTTTATGCCATTTTTCTgt                                                                      >  7:373841/1‑148 (MQ=255)
                                                                              tGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTACGGTAAAGATGCCAAAGCTTCAATTGAGCTATCTACACCTTTACGGCAAAAATCTGATAAAAACTTCATTAATAAGTTTTGTTGATCTTTTATGAAATTTTTCTTGc                                                                     >  1:90182/1‑148 (MQ=255)
                                                                                   aacaacaaATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACATTTACGGCCAAAATATGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGACATTTTTCTGGCGATaaa                                                               >  3:117363/1‑149 (MQ=255)
                                                                                     caacaaATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAAtt                                                             >  2:255650/1‑149 (MQ=255)
                                                                                     caacaaATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTAGTTGCTCTTTTATGCCATTTTTCTGGCGATAAAtt                                                             >  6:314814/1‑149 (MQ=255)
                                                                                      aacaaATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAAtt                                                             >  2:166526/1‑148 (MQ=255)
                                                                                       acaaATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTCGCGATAAATTTc                                                           <  7:193096/149‑1 (MQ=255)
                                                                                       acaaATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGACAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAAttt                                                            >  7:344340/1‑148 (MQ=255)
                                                                                        cataTAAAAGCGTATTGTGAAGTAATGATTCAGTTTTTGTTGTGATAGCTTCAATTGAGAGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCa                                                          <  8:344340/146‑1 (MQ=255)
                                                                                        caaATAAAAGCGTATTGTGACGTAATGATTCAGGTAACGATGCCAAAGCTTCAAGTGAGCGATCTACACCTTTACTGACAAAATCTGATCCAACCTGCAGTAATAAGTTTTTTTGATCTTTTATGCCATTTTTCTGGCGATAAATTTCa                                                          >  3:328124/1‑149 (MQ=255)
                                                                                               aaGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTACAATTGAGCGATCTACACCATTACGGCCAAAATCTGAACCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCAATTTTTCTGGCGATAAATTTCACGGCTGt                                                   >  8:224659/1‑149 (MQ=255)
                                                                                                   gTATTGTGACGTAATGATTACGGTAACGATGCAAAAGCTTCTATTGATCTATCTACACCTTTACGTCCTAAATCTGATCCAACCTGCAGTAATTAGTTTTGTTTCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTATTTgg                                               <  1:151422/149‑1 (MQ=255)
                                                                                                    tATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTggg                                              <  7:368692/149‑1 (MQ=255)
                                                                                                         tgACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTACTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTg                                         >  8:177249/1‑149 (MQ=255)
                                                                                                         tgACGTAATGAATCCGGTAACGAAGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGACATTTTTCTGGCGATAAATTTAACGGCTGTTTGGGATTTg                                         >  4:314506/1‑149 (MQ=255)
                                                                                                           aCGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCt                                       >  4:340228/1‑149 (MQ=255)
                                                                                                             gTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCa                                     >  3:258991/1‑149 (MQ=255)
                                                                                                             gTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCa                                     <  1:255650/149‑1 (MQ=255)
                                                                                                              tAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCAc                                    >  6:255810/1‑149 (MQ=255)
                                                                                                               aaTGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGttt                                                 >  6:262588/1‑135 (MQ=255)
                                                                                                               aaTGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGttt                                                 <  5:262588/135‑1 (MQ=255)
                                                                                                                aTGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATTAATTTCACGGCTGTTTGGGATTTGCTCACTg                                  >  2:27897/1‑149 (MQ=255)
                                                                                                                 tGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTg                                  >  4:161106/1‑148 (MQ=255)
                                                                                                                  taTTCCGGTAACGATGCCAAAGCATCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGCTCCAACCTCCAGTAAAAAGGTTTGTAGCTCGTTAAAGCCATTTTTCTGGCGATAAATATCACGCCTGTTTGCGATTTGCTCACAGTa                                >  2:65904/2‑149 (MQ=255)
                                                                                                                   aTTCCGGTAACGATGCCAAAGCTTCGATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTa                                <  7:86518/148‑1 (MQ=255)
                                                                                                                     tCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTCCAGTAATAAGTATTATTTCTCTTTTATTCAATTTTTATGGAGATAAATTTAACTGCTGTTTTGTATTTGCTCACTGTAttt                             >  6:201235/1‑149 (MQ=255)
                                                                                                                       cGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTc                           >  5:205791/1‑149 (MQ=255)
                                                                                                                         gTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATaaa                                                               >  7:33207/1‑111 (MQ=255)
                                                                                                                            aCGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCa                                     >  5:90884/1‑134 (MQ=255)
                                                                                                                               aTGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGTCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGa                   <  3:265114/149‑1 (MQ=255)
                                                                                                                               aTGACAAAGCTTAAATTGAGCGATCTACAGATTTACGGCCAAAATCTGATCCAACCTGAATTAATAAGTTTTGTTGCTCTTTTATGACATTTTTCTCGCGATAAATTTCACTGCTGTTTGGGATTTGCTCACTTTATTTTCTGTCCGGa                   <  5:201235/149‑1 (MQ=255)
                                                                                                                                   aaaaGCTTCAATTGAGCTATATACACCTTTAAGTCCAAAATATTATCAAACCTGAAGTAATAAGTTTTGTTGATCTTTTATTCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATaaa               <  7:389945/148‑1 (MQ=255)
                                                                                                                                    aaaGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATaaa               >  4:160409/1‑148 (MQ=255)
                                                                                                                                    aaaGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCAAGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATaaa               >  6:62994/1‑148 (MQ=255)
                                                                                                                                      aGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATAc            <  4:258991/149‑1 (MQ=255)
                                                                                                                                       gCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCCTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATAcc           >  5:357124/1‑149 (MQ=255)
                                                                                                                                         ttCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCgg         >  6:9593/1‑149 (MQ=255)
                                                                                                                                         ttCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAAATTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCgg         >  6:97417/1‑149 (MQ=255)
                                                                                                                                           cacttGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTGTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCGGTTTGGGAGTTGCTCACTGTATTTTCTGTCCGGATAAATACCgg         <  5:32439/144‑1 (MQ=255)
                                                                                                                                            aaTTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGggag      >  1:392650/1‑149 (MQ=255)
                                                                                                                                                 aGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGgaag  <  5:255810/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TTGCGACCGGAGAAGAAATGCACATTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAG  >  NC_000913/3806283‑3806575

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: ‑