Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,225,014 T→A 6.1% intergenic (‑144/+30) 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,0140TA 6.1% 259.8 / 5.4 82intergenic (‑144/+30)pbpG/yohCpeptidoglycan DD‑endopeptidase PbpG/putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (40/37);  new base A (2/3);  total (42/40)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.26e-01

GGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGTTG  >  NC_000913/2224866‑2225161
                                                                                                                                                    |                                                                                                                                                   
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                                                                                                                     >  1:1020175/1‑150 (MQ=255)
                               cgcACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTa                                                                                                                     >  2:174152/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:160005/150‑1 (MQ=255)
                                                          ggATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAg                                                                                          <  1:792657/150‑1 (MQ=255)
                                                          ggATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAg                                                                                          <  2:503417/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                                                                    <  1:1419293/150‑1 (MQ=255)
                                                                                ccATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGagcagc                                                                    >  2:358278/1‑150 (MQ=255)
                                                                                ccATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGagcagc                                                                    <  1:580585/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:358278/150‑1 (MQ=255)
                                                                                                aaaGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCa                                                    >  1:190423/1‑150 (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                                  >  2:882781/1‑150 (MQ=255)
                                                                                                                  gagaATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCg                                  <  1:987606/150‑1 (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)
                                                                                                                                                 aTTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCAAGCCAATCGCTAGCGTCGAACTATATATGTTAAGGCCTTCGTCTttgtt   <  2:846034/150‑1 (MQ=255)
                                                                                                                                                  tttACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTttgttg  <  2:1619070/150‑1 (MQ=255)
                                                                                                                                                  tttACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTttgttg  <  1:201327/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGCATGATGAGCATTCAGATAGTGGTTCACGCGCACGGGTTGCGCACCGCCGGAGTAAGGATTTACTGAGGCTAGCGACGCCATCATAACGAGCAAAAAGTGCGATCGTCAAAGGAGAATCGTGAGGAAATGCTGCATTGCTGACATTTACGCCAGCAATGCAACGTCAAAGAACTTTCTAGAACAAACGATAACCGTAGCCCCATAGTATAACGGTTAGGGCGAGCAGCACTTCCAGTACCAGCACGCCAATCGCCAGCGTCGAACTAGAGAAGCTAAGGCCTTCCTCTTTGTTG  >  NC_000913/2224866‑2225161

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: