Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,225,009 A→G 6.0% intergenic (‑139/+35) pbpG ← / ← yohC peptidoglycan DD‑endopeptidase PbpG/putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,225,0090AG 6.0% 250.5 / 5.0 83intergenic (‑139/+35)pbpG/yohCpeptidoglycan DD‑endopeptidase PbpG/putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (41/37);  new base G (2/3);  total (43/40)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.70e-01

TTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTC  >  NC_000913/2224862‑2225154
                                                                                                                                                   |                                                                                                                                                 
tttCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACa                                                                                                                                                 >  2:689570/1‑150 (MQ=255)
tttCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACa                                                                                                                                                 <  1:1997349/150‑1 (MQ=255)
    ggCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTa                                                                                                                                             >  1:765123/1‑150 (MQ=255)
           tgaGCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGc                                                                                                                                      >  2:92430/1‑150 (MQ=255)
              gCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAAt                                                                                                                                   >  1:503417/1‑150 (MQ=255)
               cATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATg                                                                                                                                  >  1:1585199/1‑150 (MQ=255)
                aTTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGc                                                                                                                                 <  2:1392853/150‑1 (MQ=255)
                  tCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGGCGTAAATGTCAGCAATGCag                                                                                                                               <  1:1061949/150‑2 (MQ=255)
                  tCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGGCGTAAATGTCAGCAATGCag                                                                                                                               <  2:432281/150‑2 (MQ=255)
                   cAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGGCGTAAATGTCAGCAATGCagc                                                                                                                              <  2:1439734/150‑3 (MQ=255)
                    aGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACg                                                                                                                             >  1:961217/1‑150 (MQ=255)
                       tAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCa                                                                                                                          <  1:354555/150‑1 (MQ=255)
                       tAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCa                                                                                                                          <  2:807841/150‑1 (MQ=255)
                        aGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCaa                                                                                                                         <  2:1585199/150‑1 (MQ=255)
                                aCGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACttt                                                                                                                 <  2:236103/150‑1 (MQ=255)
                                 cgcgCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTc                                                                                                                <  1:207243/150‑1 (MQ=255)
                                  gcgcACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCt                                                                                                               <  1:1603093/150‑1 (MQ=255)
                                   cgcACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTa                                                                                                              >  2:174152/1‑150 (MQ=255)
                                   cgcACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTa                                                                                                              >  1:1020175/1‑150 (MQ=255)
                                     cACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGa                                                                                                            <  1:1362318/150‑1 (MQ=255)
                                     cACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAACAATGCAACGTCAAAGAACTTTCTAGa                                                                                                            <  1:1844450/150‑1 (MQ=255)
                                         ggTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACaa                                                                                                        <  2:1052392/150‑1 (MQ=255)
                                           ttGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAAc                                                                                                      <  2:961217/150‑1 (MQ=255)
                                             gcgcACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGa                                                                                                    >  2:987606/1‑150 (MQ=255)
                                              cgcACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGAt                                                                                                   <  1:846030/150‑1 (MQ=255)
                                               gcACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATAGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATa                                                                                                  >  1:1304906/1‑150 (MQ=255)
                                                cACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTa                                                                                                              <  2:1623912/137‑1 (MQ=255)
                                                cACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTa                                                                                                              >  1:1623912/1‑137 (MQ=255)
                                                    gccgGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGt                                                                                             <  1:929292/150‑1 (MQ=255)
                                                    gccgGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAAAGATAACCGt                                                                                             >  1:865717/1‑150 (MQ=255)
                                                      cgGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAg                                                                                           >  2:1168732/1‑150 (MQ=255)
                                                       ggAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGc                                                                                          <  1:1374326/150‑1 (MQ=255)
                                                        gAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGcc                                                                                         >  1:1319154/1‑150 (MQ=255)
                                                         aGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGccc                                                                                        <  1:69204/150‑1 (MQ=255)
                                                          gTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGcccc                                                                                       >  2:1419293/1‑150 (MQ=255)
                                                            aaGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCAt                                                                                     >  2:1182768/1‑150 (MQ=255)
                                                              ggATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAg                                                                                   <  1:792657/150‑1 (MQ=255)
                                                              ggATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAg                                                                                   <  2:503417/150‑1 (MQ=255)
                                                              ggATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAg                                                                                   <  1:160005/150‑1 (MQ=255)
                                                               gATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGt                                                                                  <  2:1900366/150‑1 (MQ=255)
                                                                  ttACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAATGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGtata                                                                               <  1:982254/150‑1 (MQ=255)
                                                                   tACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATaa                                                                              >  2:1005121/1‑150 (MQ=255)
                                                                       gAGGCTAGCGCCGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGt                                                                          >  2:217497/1‑150 (MQ=255)
                                                                       gAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGt                                                                          >  1:319149/1‑150 (MQ=255)
                                                                        aGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGtt                                                                         <  2:80302/150‑1 (MQ=255)
                                                                            tAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAggg                                                                     <  2:1738689/150‑1 (MQ=255)
                                                                                gACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGag                                                                 >  2:489578/1‑150 (MQ=255)
                                                                                    ccATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGagcagc                                                             >  2:358278/1‑150 (MQ=255)
                                                                                    ccATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGagcagc                                                             <  1:580585/150‑1 (MQ=255)
                                                                                    ccATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGagcagc                                                             <  1:1419293/150‑1 (MQ=255)
                                                                                      atcatAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCAc                                                           >  1:194691/1‑150 (MQ=255)
                                                                                       tcatAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACt                                                          >  2:59042/1‑150 (MQ=255)
                                                                                         atAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTAAGGGCGAGCAGCACTTc                                                        >  2:1122807/1‑150 (MQ=255)
                                                                                          tAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTcc                                                       <  2:1629583/150‑1 (MQ=255)
                                                                                           aaCGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCa                                                      >  2:584193/1‑150 (MQ=255)
                                                                                            aCGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAg                                                     >  1:846034/1‑150 (MQ=255)
                                                                                               aGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTAc                                                  >  1:301944/1‑150 (MQ=255)
                                                                                                gCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTAcc                                                 <  2:1020175/150‑1 (MQ=255)
                                                                                                   aaaaGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGc                                              >  1:1478337/1‑150 (MQ=255)
                                                                                                    aaaGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCa                                             >  1:190423/1‑150 (MQ=255)
                                                                                                    aaaGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCa                                             <  1:358278/150‑1 (MQ=255)
                                                                                                      aGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACg                                           <  2:323141/150‑1 (MQ=255)
                                                                                                      aGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACg                                           >  2:993698/1‑150 (MQ=255)
                                                                                                          cGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCaa                                       >  1:1619070/1‑150 (MQ=255)
                                                                                                           gATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAAt                                      >  1:1530693/1‑150 (MQ=255)
                                                                                                                  aaaGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGc                               >  2:24112/1‑150 (MQ=255)
                                                                                                                   aaGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACAAGCACGCCAATCGCCAGCg                              >  1:630542/1‑150 (MQ=255)
                                                                                                                     ggAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTc                            <  2:1213415/150‑1 (MQ=255)
                                                                                                                      gagaATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCg                           <  1:987606/150‑1 (MQ=255)
                                                                                                                      gagaATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCg                           >  2:882781/1‑150 (MQ=255)
                                                                                                                       agaATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGa                          >  2:1293515/1‑150 (MQ=255)
                                                                                                                          aTCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACt                       <  2:1464987/150‑1 (MQ=255)
                                                                                                                             gTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTaga                    >  2:1386782/1‑150 (MQ=255)
                                                                                                                             gTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTaga                    >  1:1513099/1‑150 (MQ=255)
                                                                                                                              tGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTagag                   <  1:59042/150‑1 (MQ=255)
                                                                                                                                     aTGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTaa            <  2:190423/150‑1 (MQ=255)
                                                                                                                                       gctgcATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAgg          <  1:689570/150‑1 (MQ=255)
                                                                                                                                        ttgcATTGCTGGCGTAAATGTCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGc         >  1:1639435/2‑150 (MQ=255)
                                                                                                                                          gcATTGCTGGCGTAAATGTCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCt       >  1:928062/1‑150 (MQ=255)
                                                                                                                                             ttGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCttgc    <  1:174152/150‑3 (MQ=255)
                                                                                                                                               gCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCctc  >  1:204376/1‑150 (MQ=255)
                                                                                                                                               gCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCctc  >  1:788664/1‑150 (MQ=255)
                                                                                                                                               gCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCctc  >  2:1100682/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TTTCGGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTC  >  NC_000913/2224862‑2225154

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