Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,806,430 C→A 9.3% 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.3% 213.6 / 3.6 74R213L (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 (44/22);  new base A (2/6);  total (46/28)
Fisher's exact test for biased strand distribution p-value = 4.70e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGACCGGAGAAGAAATGCACATTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGA  >  NC_000913/3806286‑3806570
                                                                                                                                                |                                                                                                                                            
cGACCGGAGAAGAAATGCACATTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGGTGCCAAAGCTTCAATTGAGCGAt                                                                                                                                           <  3:165211/148‑1 (MQ=255)
cGACCGGAGAAGAAATGCACATTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGAt                                                                                                                                           <  6:155684/148‑1 (MQ=255)
     ggAGAAGAAATGCACATTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTacac                                                                                                                                     >  7:316309/1‑149 (MQ=255)
        gaagaaATGCACATTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCCCCTTATCCTCACCAACAACAAATAAAAACCAATTGAGACGTAAACATACCCGAAACCATCCCAAACCATCAAATCAGCAATCTACAcca                                                                                                                                   >  6:278786/1‑147 (MQ=255)
            aaaTGCACATTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACg                                                                                                                              >  2:151371/1‑149 (MQ=255)
             aaTTCACAATGCTCCGCACGCCGAGTATTTCAGCCAGCGCTTCAAATTTTCGCGGCATATCCTGACCACCAACAAATAAAAGCGTATTGTGACGTAATGCTTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACgg                                                                                                                             >  8:291204/1‑149 (MQ=255)
               tGCACATTGCTCAGCACGCCGAGTTTTTCTGCCAGCGCTTCATATTTTCGCGTCTTATCCTGACCAACAACAAATAAAAGCGTAATGTGACGTAATGATTCCGGTAACGATGCCATCGGTTCAATTGAGAGATCTACACATTTACGGc                                                                                                                            >  5:271725/1‑148 (MQ=255)
                 cacaTTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATc                                                                                                                                          <  4:61038/132‑1 (MQ=255)
                  acaTTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCaa                                                                                                                         >  2:332901/1‑148 (MQ=255)
                    aTTGCTCCGCCCGCCGAGGTTTTCTGCCAGCGCGTCAAATTTTCGCCGCATATCCTGCCCAACAACAAATAAAAGCGTTTTGTGACGTAATGATTTTTTTTTTTTTTTGATAGGTTGAATTGAGAGATCTACACCTTTACGGCCAAAAt                                                                                                                      <  8:106209/149‑1 (MQ=255)
                        cTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGAGCTACACCTTTACGGCCAAAATCTg                                                                                                                   <  5:278786/148‑1 (MQ=255)
                           ctcACGCCGATTAATTCTGCGAGCGATTCAAATTTTCGCGGCTTATCCTGACCACCACCATATAAAAGGGTATTGTGACGTAATGATTCAGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGa                                                                                                                  <  3:27820/144‑1 (MQ=255)
                            gCACGCCGAGTTTTTCTGCCAGAGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTTATcc                                                                                                               >  6:210100/1‑148 (MQ=255)
                             cACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCaa                                                                                                             >  8:164092/1‑149 (MQ=255)
                                gCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTTTTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCt                                                                                                          >  5:93583/1‑149 (MQ=255)
                                gCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCt                                                                                                          >  4:199426/1‑149 (MQ=255)
                                 ccGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTTAGCGATCTACACCTTTACGGCCAAAATCTGATCCAAACTg                                                                                                         >  7:278838/1‑149 (MQ=255)
                                         ttCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGtaat                                                                                                  <  3:199426/148‑1 (MQ=255)
                                          tCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGtaata                                                                                                 <  7:164092/148‑1 (MQ=255)
                                               cAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGtttt                                                                                           >  6:140536/1‑149 (MQ=255)
                                                aGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTttg                                                                                          <  3:64334/149‑1 (MQ=255)
                                                     ttCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGct                                                                                     >  2:266813/1‑149 (MQ=255)
                                                     ttCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGct                                                                                     >  7:303083/1‑149 (MQ=255)
                                                     ttCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTAGTTGct                                                                                     >  2:284312/1‑149 (MQ=255)
                                                        aaaTTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCtt                                                                                  >  2:185599/1‑149 (MQ=255)
                                                         aaTTTTCGCGGCTTATCCTGACAAACAACAAATAAAATCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGAGATCTACACCTATACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCttt                                                                                 >  5:47660/1‑149 (MQ=255)
                                                           ttttCGCGGCTTATCCTGACCAACAAGAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTttgtt                                                                                        <  5:275041/140‑1 (MQ=255)
                                                           ttttCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCTGTAACGATGCCAAAGATTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGAAACAACCTCCAGTAATAAGTTttgat                                                                                        >  6:275041/1‑138 (MQ=255)
                                                             ttCGCTGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTAt                                                                              <  3:268315/148‑1 (MQ=255)
                                                                        tCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTATg                                                                  >  3:250335/1‑149 (MQ=255)
                                                                          cttaCCAACAACAAATAAAAGCGTATTGTGACGTAATGATTTTTTTTTTTTTTTGATAGGGTGAATTGAGAGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGAAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGc                                                                <  8:278838/146‑1 (MQ=255)
                                                                          cTGACCAACAACAACTAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGACAAAGCTTAAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAAAAAGTTTTGTTGCTCTTTTATGCCATTTTTCTgg                                                                 <  5:61059/148‑1 (MQ=255)
                                                                          cTGACCAACAACAAATAAACGCGTATTGTGCCGTAATGATTTATTTTTATAATTGATAGGTTGAATTGAGAGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTgg                                                                 <  8:316309/148‑1 (MQ=255)
                                                                           tGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGt                                                                                                     >  2:330290/1‑111 (MQ=255)
                                                                           tGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGt                                                                                                     <  1:330290/111‑1 (MQ=255)
                                                                                 acaacaaATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGa                                                                                                                  >  8:265837/1‑92 (MQ=255)
                                                                                  caacaaATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAAtt                                                        >  6:284249/1‑149 (MQ=255)
                                                                                   aacaaATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAAttt                                                       >  7:264929/1‑149 (MQ=255)
                                                                                      aaaTAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCAc                                                    >  2:20065/1‑149 (MQ=255)
                                                                                          aaaaGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATaa                                                           <  5:237671/138‑1 (MQ=255)
                                                                                          aaaaGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTCCAGAAAAAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATaa                                                           >  6:237671/1‑138 (MQ=255)
                                                                                          aaaaGCGTATTGTGACGTAATGATTCCAGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCt                                                >  8:334427/1‑149 (MQ=255)
                                                                                              gCGTACTGTGCCGAAATGATTCCGGTAACGATGCCAAAGCCTCAATTGAGCGATCTACGCCTTTACGGCCAAAATCTCATCCAACCTGCAGTACTAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGttt                                            <  7:291204/149‑1 (MQ=255)
                                                                                                   ttGTGACTTAATGATTCATTTAATGATTCGATAGCTTGAATTGAGAGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGACATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGAt                                       <  2:112612/149‑1 (MQ=255)
                                                                                                      tgACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACATGCAGTAATAAGTTATGTTGCTCTTTTATGCCATTTTTCTGGCCATAAATTTCACGGCTGTTTGGGAATTg                                    >  5:45933/1‑149 (MQ=255)
                                                                                                        aCGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGc                                   <  1:332901/148‑1 (MQ=255)
                                                                                                          gTAATGATTCCGGTAACGATGGCATAGCTTCAATTGAGAGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCa                                <  8:303083/149‑1 (MQ=255)
                                                                                                           tAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCt                                                >  2:258936/1‑132 (MQ=255)
                                                                                                           tAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCt                                                <  1:258936/132‑1 (MQ=255)
                                                                                                           tAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCAc                               >  2:292417/1‑149 (MQ=255)
                                                                                                            aaTGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACt                              >  4:16862/1‑149 (MQ=255)
                                                                                                            aaTGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACt                              >  4:12276/1‑149 (MQ=255)
                                                                                                              tGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTg                             >  1:169622/1‑148 (MQ=255)
                                                                                                                aTTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGAtt                                      >  2:104289/1‑137 (MQ=255)
                                                                                                                aTTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGAtt                                      <  1:104289/137‑1 (MQ=255)
                                                                                                                aTTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTa                           <  1:151371/148‑1 (MQ=255)
                                                                                                                    cGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTAACGGCTGTATGGGATTTGCTCACTGTATTTTc                      >  8:306879/1‑149 (MQ=255)
                                                                                                                       tAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGt                   <  3:12276/149‑1 (MQ=255)
                                                                                                                       tAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGt                   <  3:16862/149‑1 (MQ=255)
                                                                                                                           gATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCa                                                                          >  2:293350/1‑90 (MQ=255)
                                                                                                                             tGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTc                  <  6:205169/144‑1 (MQ=255)
                                                                                                                             tGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTc                  >  5:205169/1‑144 (MQ=255)
                                                                                                                               ccAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATa            >  2:137138/1‑148 (MQ=255)
                                                                                                                                cAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATaaa          >  3:175692/1‑149 (MQ=255)
                                                                                                                                cAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATaaa          >  8:267314/1‑149 (MQ=255)
                                                                                                                                 aaaGCTTCAATTGCGCGATCTACACCTTTACGGCCAAAATCTGATCCACCCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATaaa          >  5:133470/1‑148 (MQ=255)
                                                                                                                                 aaaGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCTATAAATTTCACTGCTGTATGGGATTTGCTAACTGTATTTTCTTTCCTGATaaa          >  8:137186/1‑148 (MQ=255)
                                                                                                                                 aaaGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATaaa          >  1:50565/1‑148 (MQ=255)
                                                                                                                                 aaaGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATaaa          >  5:199371/1‑148 (MQ=255)
                                                                                                                                      ttCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCgg    >  2:182737/1‑149 (MQ=255)
                                                                                                                                      ttCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCgg    >  6:309406/1‑149 (MQ=255)
                                                                                                                                      ttCAATTGAGAGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCgg    <  6:93583/149‑1 (MQ=255)
                                                                                                                                       tCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCgg    >  3:272473/1‑148 (MQ=255)
                                                                                                                                         aaTTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGgga  <  5:284249/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                            
CGACCGGAGAAGAAATGCACATTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGA  >  NC_000913/3806286‑3806570

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: