Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,170,282 C→G 7.3% W31C (TGG→TGC)  ldtC ← L,D‑transpeptidase linking Lpp to murein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,170,2820CG 7.3% 270.6 / 11.0 84W31C (TGG→TGC) ldtCL,D‑transpeptidase linking Lpp to murein
Reads supporting (aligned to +/- strand):  ref base C (36/41);  new base G (3/3);  total (39/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.42e-01

GACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCAACGGGATCGTTAACACGCTGCCCGCGCGCGGTACGTAAGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCTGAATTATGGATGGTGACAGTGTCGGATAGTGCAGGGAAGTGCAAAGAATTTGTAAATGTTGCAGATGGGGGCGCAGAAACGCCCCCGATTTACCATTAAGCGATCGCGTT  >  NC_000913/1170038‑1170532
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                            
gaCGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCAACGGGATCGTTAACACGCTGCCCGCGCGCGGTACGTAAGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCAGGTagc‑gccagc                                                                                                                                                                                                                               >  2:31291/1‑275 (MQ=255)
            cAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCAACGGGATCGTTAACACGCTGCCCGCGCGCGGTACGTAAGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTTGCGGCAGCGGCCAGGGGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAA‑GGCAgc                                                                                                                                                                                                                   >  2:15755/1‑274 (MQ=255)
                    aTCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCAACGGGATCGTTAACACGCTGCCCGCGCGCGGTACGTAAGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCTGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACg                                                                                                                                                                                                           >  2:481780/1‑274 (MQ=255)
                       ccagtGCCTTCGGGGGGGGCATCGGGAAGTAGGGTTTGCAACGGGAGCGTTAACACGCGGCCCGCGGGCGGTACGTANGGGTCAAAGCCGGGGTTAGCCTGTAACAGAGCGGGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTaaa                                                                                                                                                                                                       <  3:163180/271‑1 (MQ=255)
                           tGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCAACGGGATCGTTAACACGCTGCCCGCGCGCGGTACGTAAGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAT‑GGCAGCGGCGAACGTAAAAAAc                                                                                                                                                                                                   >  2:489372/1‑275 (MQ=255)
                                tCGCGCGGCGCATCTGGAAGTAGGGTTTGCAACGGGATCGTTAACACGCTGCCCGCGCGCGGTACGTAAGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACAACTGCCCGATGGCGGCAGCGGCCAGGGGTTGGCTTTTTCCGG‑AGCAGCCAGCGCCCCGGGCAGC‑GCGAACGAaaaacaaggttg                                                                                                                                                                                              >  4:494020/1‑267 (MQ=255)
                                 cgcgcgGCGCATCTGGAAGTAGGGTTTGCAACGGGATCGTTAACACGCTGCCCGCGCGCGGTACGTAAGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGACCGTAAAAAACGTTAg                                                                                                                                                                                              >  2:376543/1‑274 (MQ=255)
                                       gcgcATCTGGAAGTAGGGTTTGCAACGGGATCGTTAACACGCTGCCCGCGCGCGGTACGTAAGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTCTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCga                                                                                                                                                                                       >  2:643946/1‑275 (MQ=255)
                                                                       gTTAACACGCTGCCCGCGCGCGGTACGTAAGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGat                                                                                                                                                                      >  2:13028/1‑260 (MQ=255)
                                                                            cccGCTGCCCGCGCGCGGTACGTAAGGATCAACGACGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCataat                                                                                                                                                  <  2:494022/273‑1 (MQ=255)
                                                                            cacGCTGCCCGGGCGCGGTAAGTAATGATCAACGCCCGGGTTAGCCTGTTCCAGAGCGAGCAAGCCGCCGTTGTATTTTTTTGCGATGGCTTCCAGAGAACCACCTTCATTTTCCACCGCGTGAAATTTG‑TTTTCGCCAACCAGACGACGGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCataat                                                                                                                                                  <  4:489370/275‑1 (MQ=255)
                                                                                 tGCCCGCGCGCGGTACGTAAGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCACCCAGACGACTGCCCGCTGGCGGCCGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCCGCGGAGAACGTAAAAAACGCTAGACAGCGACAAACAGGCGTCTTGAGCTGACCCAACAATCAATaataaaca                                                                                                                                             >  4:83303/1‑273 (MQ=255)
                                                                                                 gTAAGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGcc                                                                                                                                                                                                                                                             >  8:463028/1‑149 (MQ=255)
                                                                                                 gTAAGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGc                                                                                                                                                                                                                                                              <  5:334807/148‑1 (MQ=255)
                                                                                                  tAAGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCCGACGACTGCCCGCTGGCGGCAGCGGTCa                                                                                                                                                                                                                                                            <  7:112557/149‑1 (MQ=255)
                                                                                                   aaGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAg                                                                                                                                                                                                                                                           >  8:327415/1‑149 (MQ=255)
                                                                                                   aaGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAg                                                                                                                                                                                                                                                           >  7:636782/1‑149 (MQ=255)
                                                                                                     ggATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAgg                                                                                                                                                                                                                                                          >  6:57456/1‑148 (MQ=255)
                                                                                                      gATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGCACGTCAACAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAg                                                                                                                        >  2:293907/1‑275 (MQ=255)
                                                                                                       tgtaatgCCGGGGTTCGCCTTGAATAGGGCGAGGAAGCCGACGTTTTACTTTTTGGCGATGGCTTCCAGAGAACAACCGTCATTTTCCACCACATGAGATTTG‑TTTTCGCCAACCAGACGACCGCCCGCTGGCGGCAGCTGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGa                                                                                                                       <  4:376541/269‑1 (MQ=255)
                                                                                                        tCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCa                                                                                                                                                      >  2:189723/1‑243 (MQ=255)
                                                                                                        tCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAAAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACCGAt                                                                                                                      >  2:594106/1‑275 (MQ=255)
                                                                                                          aaCGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTgg                                                                                                                                                                                                                                                    >  8:396301/1‑149 (MQ=255)
                                                                                                          aaCGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTgg                                                                                                                                                                                                                                                    >  7:67981/1‑149 (MQ=255)
                                                                                                          aaCGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTgg                                                                                                                                                                                                                                                    <  8:373801/149‑1 (MQ=255)
                                                                                                          aaCGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCa                                                                                                                                                      >  2:271849/1‑241 (MQ=255)
                                                                                                           aCGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTGATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTgg                                                                                                                                                                                                                                                    >  6:431226/1‑148 (MQ=255)
                                                                                                              ccGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCataat                                                                                                                                                  <  4:271848/241‑1 (MQ=255)
                                                                                                                     tAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCTGAATTa                                                                                                          <  2:496025/274‑1 (MQ=255)
                                                                                                                      aGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGt                                                                                                                                                                                                                                         <  5:429320/148‑1 (MQ=255)
                                                                                                                      aGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTa                                                                                                                                                                                                                                        <  7:210166/149‑1 (MQ=255)
                                                                                                                      aGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCAAGAAAAACGCGTTTTGATCATCCCACAAAATcct                                                                                                                                                      >  4:170138/1‑228 (MQ=255)
                                                                                                                          tGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACAACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCataat                                                                                                                                                  <  2:170139/229‑1 (MQ=255)
                                                                                                                               caaagggggaAAGCCGATGTTGTATTATTTGGGAAGGGCTACCAGGGAACCACCGTGATTTTCCACCACATGAAATTTG‑TTTGCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAAGTGTTGGCTTTTGACGGTAGC‑GCCAGCGCAAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCTGAATTATGGATGGTGa                                                                                                <  4:594104/266‑1 (MQ=255)
                                                                                                                               aagagCGATAAATCCGACGTTGTATTGTCTGGAGATAGCTGCCAGAGAACAACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCTGAATTATGGATGGTGa                                                                                                <  2:83304/273‑1 (MQ=255)
                                                                                                                                 gagCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTagc‑gccagcgcc                                                                                                                                                                                                                            >  7:885336/1‑149 (MQ=255)
                                                                                                                                  agCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCa                                                                                                                                                      >  2:814373/1‑217 (MQ=255)
                                                                                                                                     gagaAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTGATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑g                                                                                                                                                                                                                        <  5:431226/148‑1 (MQ=255)
                                                                                                                                     gagaAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑gg                                                                                                                                                                                                                       <  7:396301/149‑1 (MQ=255)
                                                                                                                                      ngaAANCCNACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCNCCCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCataat                                                                                                                                                  <  4:814371/216‑1 (MQ=255)
                                                                                                                                      agaAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGc                                                                                                                                                                                                                                                              <  4:767985/111‑1 (MQ=255)
                                                                                                                                      agaAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGaa                                                                                                                                                                                    >  4:425093/1‑183 (MQ=255)
                                                                                                                                      agaAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCa                                                                                                                                                      <  2:546248/213‑1 (MQ=255)
                                                                                                                                        aaggccgcagtgttATTTTTTTGCGTTTGGTTCCAGAGACCGACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCTGAATTATGGATGGTGACAGTGTCGGa                                                                                      <  4:481778/262‑1 (MQ=255)
                                                                                                                                        aaaaCCTACGTTGTATTTTTTGGCGATGGCTGCCAGAGAGCCACCGTCATTTTCCACCACATGAAATTTG‑TTTACGCCGACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCTGAATTATGGATGGTGACAGTGTCGGa                                                                                      <  4:293906/275‑1 (MQ=255)
                                                                                                                                          aGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAAc                                                                                                                                                                                <  2:425095/183‑1 (MQ=255)
                                                                                                                                                                gATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGc                                                                                                                                                                                             <  5:518308/148‑1 (MQ=255)
                                                                                                                                                                 aTGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGcc                                                                                                                                                                                            <  6:63669/148‑1 (MQ=255)
                                                                                                                                                                   ggCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGc                                                                                                                                                                                         >  7:661845/1‑149 (MQ=255)
                                                                                                                                                                    gCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAg                                                                                                                                                                                          >  6:7076/1‑148 (MQ=255)
                                                                                                                                                                      ttCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCgag                                                                                                                                                                                      >  7:639548/1‑149 (MQ=255)
                                                                                                                                                                            agaACCACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGaaaaa                                                                                                                                                                                 <  5:57456/148‑1 (MQ=255)
                                                                                                                                                                                  accGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGttt                                                                                                                                                                          <  8:254755/149‑1 (MQ=255)
                                                                                                                                                                                     gTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGa                                                                                                                                                                       <  8:639548/149‑1 (MQ=255)
                                                                                                                                                                                      tCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGat                                                                                                                                                                      <  8:636782/149‑1 (MQ=255)
                                                                                                                                                                                         ttttCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGatca                                                                                                                                                                    >  5:149248/1‑148 (MQ=255)
                                                                                                                                                                                           ttCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGAtcatca                                                                                                                                                                 <  8:885336/149‑1 (MQ=255)
                                                                                                                                                                                                ccacATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAg                                                                                                                                         >  4:457449/1‑168 (MQ=255)
                                                                                                                                                                                                    aTGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAAt                                                                                                                                                         <  5:912588/148‑1 (MQ=255)
                                                                                                                                                                                                    aTGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTa                                                                                                                                     <  2:457451/168‑1 (MQ=255)
                                                                                                                                                                                                     tGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATc                                                                                                                                                        >  6:755055/1‑148 (MQ=255)
                                                                                                                                                                                                      gAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCacaacccc                                                                                                                                                       <  6:149248/148‑5 (MQ=255)
                                                                                                                                                                                                      gAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATcc                                                                                                                                                       <  5:755055/148‑1 (MQ=255)
                                                                                                                                                                                                      gAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATcc                                                                                                                                                       <  5:650502/148‑1 (MQ=255)
                                                                                                                                                                                                           ttG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCataat                                                                                                                                                  >  6:365889/1‑148 (MQ=255)
                                                                                                                                                                                                           ttG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCataat                                                                                                                                                  >  5:780976/1‑148 (MQ=255)
                                                                                                                                                                                                                   cGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATaa                                                                                                                                          <  7:463028/149‑1 (MQ=255)
                                                                                                                                                                                                                    gCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAg                                                                                                                                         <  8:67981/149‑1 (MQ=255)
                                                                                                                                                                                                                           caacaccTGGCCGCTGCCGCCAGCGGGCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTg                                                                                                                                  >  8:826173/7‑149 (MQ=255)
                                                                                                                                                                                                                              cgacTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAACGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCTGAATTATGGATGGTGACAGTGTCGGATAGTGCAGGGAAGTGCAAAGAATTTGTAAATGTTGCAGATGGGGGCGCAGAAACGCCCCCGATTTACCATTAAGCGATCGCGtt  <  1:133454/274‑1 (MQ=255)
                                                                                                                                                                                                                                acTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGtaat                                                                                                                              >  5:263399/1‑148 (MQ=255)
                                                                                                                                                                                                                                  tggccGCTGCCGCCAGCGGGCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGtaataa                                                                                                                            >  5:100954/4‑148 (MQ=255)
                                                                                                                                                                                                                                  tGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGTAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATcc                                                                                                                                                       >  1:49577/1‑121 (MQ=255)
                                                                                                                                                                                                                                  tGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGtaataa                                                                                                                            <  5:365889/148‑1 (MQ=255)
                                                                                                                                                                                                                                   ggccGCTGCCGCCAGCGGGCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTg                                                                                                                                                                        >  6:841049/3‑103 (MQ=255)
                                                                                                                                                                                                                                   ggccGCTGCCGCCAGCGGGCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTg                                                                                                                                                                        <  5:841049/101‑1 (MQ=255)
                                                                                                                                                                                                                                    cccGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACa                                                                                                                         <  7:94713/149‑1 (MQ=255)
                                                                                                                                                                                                                                    cccGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATNGTAATAAACAGATAAGCCTGAAt                                                                                                            >  3:99398/1‑162 (MQ=255)
                                                                                                                                                                                                                                      cGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCataa                                                                                                                                                   <  3:49576/121‑1 (MQ=255)
                                                                                                                                                                                                                                       gCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGAt                                                                                                                      <  7:327415/149‑1 (MQ=255)
                                                                                                                                                                                                                                        ctgccgccAGCGGGCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATa                                                                                                                     <  7:826173/142‑1 (MQ=255)
                                                                                                                                                                                                                                        cTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCTGAATTatg                                                                                                        <  1:99399/162‑1 (MQ=255)
                                                                                                                                                                                                                                         tgccgccAGCGGGCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATa                                                                                                                     <  6:100954/142‑1 (MQ=255)
                                                                                                                                                                                                                                         tGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAAc                                                                                                                          >  4:808389/1‑143 (MQ=255)
                                                                                                                                                                                                                                           gcggcAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAg                                                                                                                   <  6:630346/148‑1 (MQ=255)
                                                                                                                                                                                                                                             ggcAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGAt                                                                                                                      <  2:808391/143‑1 (MQ=255)
                                                                                                                                                                                                                                                    gCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCTGAATTa                                                                                                          >  6:163394/1‑148 (MQ=255)
                                                                                                                                                                                                                                                     ccAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCa                                                                                                                                                      >  3:64515/1‑103 (MQ=255)
                                                                                                                                                                                                                                                     ccAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCTGAATTat                                                                                                         >  6:433547/1‑148 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                            
GACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCAACGGGATCGTTAACACGCTGCCCGCGCGCGGTACGTAAGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTG‑TTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGC‑GCCAGCGCCAC‑GGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCTGAATTATGGATGGTGACAGTGTCGGATAGTGCAGGGAAGTGCAAAGAATTTGTAAATGTTGCAGATGGGGGCGCAGAAACGCCCCCGATTTACCATTAAGCGATCGCGTT  >  NC_000913/1170038‑1170532

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑