Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,806,430 C→A 9.8% R213L (CGC→CTC)  waaG ← lipopolysaccharide glucosyltransferase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,806,4300CA 9.8% 246.4 / 3.8 81R213L (CGC→CTC) waaGlipopolysaccharide glucosyltransferase I
Reads supporting (aligned to +/- strand):  ref base C (48/24);  new base A (2/7);  total (50/31)
Fisher's exact test for biased strand distribution p-value = 2.35e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTGCGACCGGAGAAGAAATGCACATTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAA  >  NC_000913/3806282‑3806577
                                                                                                                                                    |                                                                                                                                                   
atTGCGACCGGAGAAGAAATGCACATTGCTCCGCACGCCGAGTATTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCATGACCAACAACAAAAAAAAGAGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCg                                                                                                                                                    >  2:1277316/1‑150 (MQ=255)
 ttGCGACCGGAGAAGAAATGCACATTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGa                                                                                                                                                   <  1:21840/150‑1 (MQ=255)
   gCGACCGGAGAAGAAATGCACATTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATc                                                                                                                                                 <  1:608452/150‑1 (MQ=255)
     gACCGGAGAAGAAATGCACATTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTa                                                                                                                                               >  1:1246043/1‑150 (MQ=255)
      aCCGGAGAAGAAATGCACATTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTac                                                                                                                                              <  2:1309888/150‑1 (MQ=255)
      aCCGGAGAAGAAATGCACATTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTac                                                                                                                                              <  1:620043/150‑1 (MQ=255)
       ccGGAGAAGAAATGCACATTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTaca                                                                                                                                             <  2:1076631/150‑1 (MQ=255)
         ggAGAAGAAATGCACATTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACAcc                                                                                                                                           >  2:1022414/1‑150 (MQ=255)
            gaagaaATGCACATTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCttt                                                                                                                                        >  2:838463/1‑150 (MQ=255)
              agaaATGCACATTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTAc                                                                                                                                      >  1:1495411/1‑150 (MQ=255)
              agaaATGCACATTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTAc                                                                                                                                      <  1:838463/150‑1 (MQ=255)
                aaaTGCACATTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACgg                                                                                                                                    <  1:896898/150‑1 (MQ=255)
                    gCACATTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCaa                                                                                                                                >  1:1002540/1‑150 (MQ=255)
                     cacaTTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCaaa                                                                                                                               <  1:125341/150‑1 (MQ=255)
                         ttGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCAGGTAACGATGCCAAAGCTTCAATTGAGAGATCTACACCTTTACGGCCAAAATCt                                                                                                                           >  1:1122275/1‑150 (MQ=255)
                             tCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATc                                                                                                                       <  1:352757/150‑1 (MQ=255)
                                gCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCaa                                                                                                                    >  1:246983/1‑150 (MQ=255)
                                  aCGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATCAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAAcc                                                                                                                  <  1:1181013/150‑1 (MQ=255)
                                    gCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTAAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCGAACCTg                                                                                                                >  1:797584/1‑150 (MQ=255)
                                     ccGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGc                                                                                                               <  1:870389/150‑1 (MQ=255)
                                      cGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTTTGACGTAATGATTACGGTAACGATGCAAAAGCTTCAATTGAGAGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCa                                                                                                              >  1:1006304/1‑150 (MQ=255)
                                        aGTTTTTCTGACAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGt                                                                                                            <  1:952199/150‑1 (MQ=255)
                                             ttCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTTTTTTTTTTTATTGATAGGGTGAATTGAGAGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGtaataa                                                                                                       <  2:1431127/150‑1 (MQ=255)
                                             ttCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGtaataa                                                                                                       >  1:1226585/1‑150 (MQ=255)
                                                 gCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGttt                                                                                                   >  2:160622/1‑150 (MQ=255)
                                                   cAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTttg                                                                                                 >  1:999837/1‑150 (MQ=255)
                                                    aGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTttgt                                                                                                >  2:646625/1‑150 (MQ=255)
                                                     tcgcTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTTTTTTTTTTTATTGATAGGGTGAATTGAGAGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTttgtt                                                                                               <  2:630379/149‑1 (MQ=255)
                                                     ccgcTTCAAATTTTCGCCGCTTATCCTGACCAACAACAAATAAAATCGTATTGTGACGTAATGATTTTTTTTTTTTATTGATAGGGTGAATAGAGAGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTttgtt                                                                                               <  2:1419850/149‑1 (MQ=255)
                                                       gcTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGc                                                                                             >  1:807390/1‑150 (MQ=255)
                                                         ttCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTACAATTCAGCGATCTACACCTTTACGGAAAAAAACTGATCAAACCTGCAGTAATAAGATTTTTTGcta                                                                                           >  1:1162491/1‑149 (MQ=255)
                                                          tCAAATTTTCGCGGCTTATCCTGACCACCAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATATGATCCAACCTGCAGTAATAAGTTTTGTTGctct                                                                                          >  1:1504107/1‑150 (MQ=255)
                                                          tCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGctct                                                                                          >  1:192686/1‑150 (MQ=255)
                                                              aTTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTa                                                                                      >  1:622299/1‑150 (MQ=255)
                                                                      ggCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTTTTTTTTTTTTATGATAGGGTGAATTGAGAGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTGAGGCCAtgg                                                                              <  2:98298/150‑3 (MQ=255)
                                                                            tCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTgg                                                                        >  1:256427/1‑150 (MQ=255)
                                                                            tCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTgg                                                                        >  2:777329/1‑150 (MQ=255)
                                                                             ccTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGc                                                                       >  1:1001266/1‑150 (MQ=255)
                                                                                  ccAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATaa                                                                  >  1:621702/1‑150 (MQ=255)
                                                                                        acaaATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCa                                                            >  1:476628/1‑150 (MQ=255)
                                                                                           aaTAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGAt                                                                    <  1:750594/139‑1 (MQ=255)
                                                                                           aaTAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGAt                                                                    >  2:750594/1‑139 (MQ=255)
                                                                                           aaTAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAtttaacccg                                                         >  2:393140/1‑144 (MQ=255)
                                                                                           aaTAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTAAACCTTTACGGACAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTATTTTATGCCATTTTTCTGTCGATAAATTTCACgg                                                         >  1:1291486/1‑150 (MQ=255)
                                                                                             tAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCt                                                       >  2:470209/1‑150 (MQ=255)
                                                                                             tAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCt                                                       >  2:1368862/1‑150 (MQ=255)
                                                                                                  gCGTATTGTGACGTAATGATTATTTTAATTATTTGAAAGGTTGAATTGAGAGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTg                                                  <  2:256427/150‑1 (MQ=255)
                                                                                                   cGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTgg                                                 >  2:78165/1‑150 (MQ=255)
                                                                                                    gTATTGTGACGTAATGATTCCGTTAAATATGTGATAGCTTGAATTGAGAGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTggg                                                <  2:192686/150‑1 (MQ=255)
                                                                                                    gTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCGCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTggg                                                <  1:24821/150‑1 (MQ=255)
                                                                                                      aTTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGTAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGAt                                              >  1:397802/1‑150 (MQ=255)
                                                                                                           gACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCt                                         >  2:103681/1‑150 (MQ=255)
                                                                                                            aCGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTc                                        >  1:953084/1‑150 (MQ=255)
                                                                                                               tAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACt                                     <  1:470209/150‑1 (MQ=255)
                                                                                                                 aTGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTATGGCGATAAATTTCACGGCTGTTTGGGATTTGCTAACTGt                                   >  1:524107/1‑150 (MQ=255)
                                                                                                                   gATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTAt                                 <  1:103681/150‑1 (MQ=255)
                                                                                                                   gATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTAt                                 >  1:90796/1‑150 (MQ=255)
                                                                                                                    aTTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTAtt                                <  1:1277316/150‑1 (MQ=255)
                                                                                                                     ttCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTAttt                               <  1:69331/150‑1 (MQ=255)
                                                                                                                       ccGTTAACTATGGCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTc                             <  2:476628/150‑1 (MQ=255)
                                                                                                                         ggTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGGGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTg                           >  1:866804/1‑150 (MQ=255)
                                                                                                                         ggTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTg                           >  2:1260462/1‑150 (MQ=255)
                                                                                                                         ggTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTg                           >  2:1266638/1‑150 (MQ=255)
                                                                                                                            aaCGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCAAGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTcc                        >  1:523981/1‑150 (MQ=255)
                                                                                                                             aCGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCAAGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCg                       >  1:1092546/1‑150 (MQ=255)
                                                                                                                             aCGATGACAAAGCTACAAATGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCg                       >  2:1238398/1‑150 (MQ=255)
                                                                                                                              cGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCgg                      >  1:1080385/1‑150 (MQ=255)
                                                                                                                              cGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCgg                      >  1:792013/1‑150 (MQ=255)
                                                                                                                               gATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGa                     <  1:646625/150‑1 (MQ=255)
                                                                                                                                  gCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATaa                  <  1:160622/150‑1 (MQ=255)
                                                                                                                                    cAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAAt                >  1:1462550/1‑150 (MQ=255)
                                                                                                                                    cAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAAt                >  1:672166/1‑150 (MQ=255)
                                                                                                                                        gCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCg            >  2:1011245/1‑150 (MQ=255)
                                                                                                                                         cTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCgg                      <  2:1470520/139‑1 (MQ=255)
                                                                                                                                         cTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCg                       >  1:848911/1‑138 (MQ=255)
                                                                                                                                         cTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCg                       <  2:848911/138‑1 (MQ=255)
                                                                                                                                         cTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGATCACTGTATTTTCTGTCCgg                      >  1:1470520/1‑139 (MQ=255)
                                                                                                                                         cTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGAAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTCTTTGGGATTTGCTCACTGTATTTTCTGTCCCGATAAATACCgg           >  2:264914/1‑150 (MQ=255)
                                                                                                                                          ttgaaTTGAGAGATCTACAACTTTACGGCCAAAATCTGATCCAAACTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCggg          <  2:797584/147‑1 (MQ=255)
                                                                                                                                            cAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGggag        <  1:1260462/150‑1 (MQ=255)
                                                                                                                                                  aGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGgaagaa  >  2:1493216/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATTGCGACCGGAGAAGAAATGCACATTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTGACGTAATGATTCCGGTAACGATGCCAAAGCTTCAATTGAGCGATCTACACCTTTACGGCCAAAATCTGATCCAACCTGCAGTAATAAGTTTTGTTGCTCTTTTATGCCATTTTTCTGGCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAA  >  NC_000913/3806282‑3806577

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: