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→TGCldtCL,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

GACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCAACGGGATCGTTAACACGCTGCCCGCGCGCGGTACGTAAGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCTGAATTATGGATGGTGACAGTGTCGGATAGTGCAGGGAAGTGCAAAGAATTTGTAAATGTTGCAGATGGGGGCGCAGAAACGCCCCCGATTTACCATTAAGCGATCGCGTT  >  NC_000913/1170038‑1170532
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                            
gaCGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCAACGGGATCGTTAACACGCTGCCCGCGCGCGGTACGTAAGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCAGGTagcgccagc                                                                                                                                                                                                                               >  2:31291/1‑275 (MQ=255)
            cAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCAACGGGATCGTTAACACGCTGCCCGCGCGCGGTACGTAAGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTTGCGGCAGCGGCCAGGGGTTGGCTTTTGCCGGTAGCGCCAGCGCCAAGGCAgc                                                                                                                                                                                                                   >  2:15755/1‑274 (MQ=255)
                    aTCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCAACGGGATCGTTAACACGCTGCCCGCGCGCGGTACGTAAGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCTGGTAGCGCCAGCGCCACGGCAGCGGCGAACg                                                                                                                                                                                                           >  2:481780/1‑274 (MQ=255)
                       ccagtGCCTTCGGGGGGGGCATCGGGAAGTAGGGTTTGCAACGGGAGCGTTAACACGCGGCCCGCGGGCGGTACGTANGGGTCAAAGCCGGGGTTAGCCTGTAACAGAGCGGGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTaaa                                                                                                                                                                                                       <  3:163180/271‑1 (MQ=255)
                           tGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCAACGGGATCGTTAACACGCTGCCCGCGCGCGGTACGTAAGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCATGGCAGCGGCGAACGTAAAAAAc                                                                                                                                                                                                   >  2:489372/1‑275 (MQ=255)
                                tCGCGCGGCGCATCTGGAAGTAGGGTTTGCAACGGGATCGTTAACACGCTGCCCGCGCGCGGTACGTAAGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACAACTGCCCGATGGCGGCAGCGGCCAGGGGTTGGCTTTTTCCGGAGCAGCCAGCGCCCCGGGCAGCGCGAACGAaaaacaaggttg                                                                                                                                                                                              >  4:494020/1‑267 (MQ=255)
                                 cgcgcgGCGCATCTGGAAGTAGGGTTTGCAACGGGATCGTTAACACGCTGCCCGCGCGCGGTACGTAAGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGACCGTAAAAAACGTTAg                                                                                                                                                                                              >  2:376543/1‑274 (MQ=255)
                                       gcgcATCTGGAAGTAGGGTTTGCAACGGGATCGTTAACACGCTGCCCGCGCGCGGTACGTAAGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTCTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCga                                                                                                                                                                                       >  2:643946/1‑275 (MQ=255)
                                                                       gTTAACACGCTGCCCGCGCGCGGTACGTAAGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGat                                                                                                                                                                      >  2:13028/1‑260 (MQ=255)
                                                                            cccGCTGCCCGCGCGCGGTACGTAAGGATCAACGACGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCataat                                                                                                                                                  <  2:494022/273‑1 (MQ=255)
                                                                            cacGCTGCCCGGGCGCGGTAAGTAATGATCAACGCCCGGGTTAGCCTGTTCCAGAGCGAGCAAGCCGCCGTTGTATTTTTTTGCGATGGCTTCCAGAGAACCACCTTCATTTTCCACCGCGTGAAATTTGTTTTCGCCAACCAGACGACGGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCataat                                                                                                                                                  <  4:489370/275‑1 (MQ=255)
                                                                                 tGCCCGCGCGCGGTACGTAAGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCACCCAGACGACTGCCCGCTGGCGGCCGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCCGCGGAGAACGTAAAAAACGCTAGACAGCGACAAACAGGCGTCTTGAGCTGACCCAACAATCAATaataaaca                                                                                                                                             >  4:83303/1‑273 (MQ=255)
                                                                                                 gTAAGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGcc                                                                                                                                                                                                                                                             >  8:463028/1‑149 (MQ=255)
                                                                                                 gTAAGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGc                                                                                                                                                                                                                                                              <  5:334807/148‑1 (MQ=255)
                                                                                                  tAAGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCCGACGACTGCCCGCTGGCGGCAGCGGTCa                                                                                                                                                                                                                                                            <  7:112557/149‑1 (MQ=255)
                                                                                                   aaGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAg                                                                                                                                                                                                                                                           >  8:327415/1‑149 (MQ=255)
                                                                                                   aaGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAg                                                                                                                                                                                                                                                           >  7:636782/1‑149 (MQ=255)
                                                                                                     ggATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAgg                                                                                                                                                                                                                                                          >  6:57456/1‑148 (MQ=255)
                                                                                                      gATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGCACGTCAACAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAg                                                                                                                        >  2:293907/1‑275 (MQ=255)
                                                                                                       tgtaatgCCGGGGTTCGCCTTGAATAGGGCGAGGAAGCCGACGTTTTACTTTTTGGCGATGGCTTCCAGAGAACAACCGTCATTTTCCACCACATGAGATTTGTTTTCGCCAACCAGACGACCGCCCGCTGGCGGCAGCTGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGa                                                                                                                       <  4:376541/269‑1 (MQ=255)
                                                                                                        tCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCa                                                                                                                                                      >  2:189723/1‑243 (MQ=255)
                                                                                                        tCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAAAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACCGAt                                                                                                                      >  2:594106/1‑275 (MQ=255)
                                                                                                          aaCGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTgg                                                                                                                                                                                                                                                    >  8:396301/1‑149 (MQ=255)
                                                                                                          aaCGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTgg                                                                                                                                                                                                                                                    >  7:67981/1‑149 (MQ=255)
                                                                                                          aaCGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTgg                                                                                                                                                                                                                                                    <  8:373801/149‑1 (MQ=255)
                                                                                                          aaCGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCa                                                                                                                                                      >  2:271849/1‑241 (MQ=255)
                                                                                                           aCGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTGATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTgg                                                                                                                                                                                                                                                    >  6:431226/1‑148 (MQ=255)
                                                                                                              ccGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCataat                                                                                                                                                  <  4:271848/241‑1 (MQ=255)
                                                                                                                     tAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCTGAATTa                                                                                                          <  2:496025/274‑1 (MQ=255)
                                                                                                                      aGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGt                                                                                                                                                                                                                                         <  5:429320/148‑1 (MQ=255)
                                                                                                                      aGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTa                                                                                                                                                                                                                                        <  7:210166/149‑1 (MQ=255)
                                                                                                                      aGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCAAGAAAAACGCGTTTTGATCATCCCACAAAATcct                                                                                                                                                      >  4:170138/1‑228 (MQ=255)
                                                                                                                          tGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACAACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCataat                                                                                                                                                  <  2:170139/229‑1 (MQ=255)
                                                                                                                               caaagggggaAAGCCGATGTTGTATTATTTGGGAAGGGCTACCAGGGAACCACCGTGATTTTCCACCACATGAAATTTGTTTGCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAAGTGTTGGCTTTTGACGGTAGCGCCAGCGCAACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCTGAATTATGGATGGTGa                                                                                                <  4:594104/266‑1 (MQ=255)
                                                                                                                               aagagCGATAAATCCGACGTTGTATTGTCTGGAGATAGCTGCCAGAGAACAACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCTGAATTATGGATGGTGa                                                                                                <  2:83304/273‑1 (MQ=255)
                                                                                                                                 gagCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTagcgccagcgcc                                                                                                                                                                                                                            >  7:885336/1‑149 (MQ=255)
                                                                                                                                  agCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCa                                                                                                                                                      >  2:814373/1‑217 (MQ=255)
                                                                                                                                     gagaAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTGATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACg                                                                                                                                                                                                                        <  5:431226/148‑1 (MQ=255)
                                                                                                                                     gagaAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACgg                                                                                                                                                                                                                       <  7:396301/149‑1 (MQ=255)
                                                                                                                                      ngaAANCCNACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCNCCCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCataat                                                                                                                                                  <  4:814371/216‑1 (MQ=255)
                                                                                                                                      agaAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGc                                                                                                                                                                                                                                                              <  4:767985/111‑1 (MQ=255)
                                                                                                                                      agaAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGaa                                                                                                                                                                                    >  4:425093/1‑183 (MQ=255)
                                                                                                                                      agaAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCa                                                                                                                                                      <  2:546248/213‑1 (MQ=255)
                                                                                                                                        aaggccgcagtgttATTTTTTTGCGTTTGGTTCCAGAGACCGACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCTGAATTATGGATGGTGACAGTGTCGGa                                                                                      <  4:481778/262‑1 (MQ=255)
                                                                                                                                        aaaaCCTACGTTGTATTTTTTGGCGATGGCTGCCAGAGAGCCACCGTCATTTTCCACCACATGAAATTTGTTTACGCCGACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCTGAATTATGGATGGTGACAGTGTCGGa                                                                                      <  4:293906/275‑1 (MQ=255)
                                                                                                                                          aGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAAc                                                                                                                                                                                <  2:425095/183‑1 (MQ=255)
                                                                                                                                                                gATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGc                                                                                                                                                                                             <  5:518308/148‑1 (MQ=255)
                                                                                                                                                                 aTGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGcc                                                                                                                                                                                            <  6:63669/148‑1 (MQ=255)
                                                                                                                                                                   ggCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGc                                                                                                                                                                                         >  7:661845/1‑149 (MQ=255)
                                                                                                                                                                    gCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAg                                                                                                                                                                                          >  6:7076/1‑148 (MQ=255)
                                                                                                                                                                      ttCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCgag                                                                                                                                                                                      >  7:639548/1‑149 (MQ=255)
                                                                                                                                                                            agaACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGaaaaa                                                                                                                                                                                 <  5:57456/148‑1 (MQ=255)
                                                                                                                                                                                  accGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGttt                                                                                                                                                                          <  8:254755/149‑1 (MQ=255)
                                                                                                                                                                                     gTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGa                                                                                                                                                                       <  8:639548/149‑1 (MQ=255)
                                                                                                                                                                                      tCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGat                                                                                                                                                                      <  8:636782/149‑1 (MQ=255)
                                                                                                                                                                                         ttttCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGatca                                                                                                                                                                    >  5:149248/1‑148 (MQ=255)
                                                                                                                                                                                           ttCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGAtcatca                                                                                                                                                                 <  8:885336/149‑1 (MQ=255)
                                                                                                                                                                                                ccacATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAg                                                                                                                                         >  4:457449/1‑168 (MQ=255)
                                                                                                                                                                                                    aTGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAAt                                                                                                                                                         <  5:912588/148‑1 (MQ=255)
                                                                                                                                                                                                    aTGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTa                                                                                                                                     <  2:457451/168‑1 (MQ=255)
                                                                                                                                                                                                     tGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATc                                                                                                                                                        >  6:755055/1‑148 (MQ=255)
                                                                                                                                                                                                      gAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCacaacccc                                                                                                                                                       <  6:149248/148‑5 (MQ=255)
                                                                                                                                                                                                      gAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATcc                                                                                                                                                       <  5:755055/148‑1 (MQ=255)
                                                                                                                                                                                                      gAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATcc                                                                                                                                                       <  5:650502/148‑1 (MQ=255)
                                                                                                                                                                                                           ttGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCataat                                                                                                                                                  >  6:365889/1‑148 (MQ=255)
                                                                                                                                                                                                           ttGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCataat                                                                                                                                                  >  5:780976/1‑148 (MQ=255)
                                                                                                                                                                                                                   cGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATaa                                                                                                                                          <  7:463028/149‑1 (MQ=255)
                                                                                                                                                                                                                    gCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAg                                                                                                                                         <  8:67981/149‑1 (MQ=255)
                                                                                                                                                                                                                           caacaccTGGCCGCTGCCGCCAGCGGGCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTg                                                                                                                                  >  8:826173/7‑149 (MQ=255)
                                                                                                                                                                                                                              cgacTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAACGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCTGAATTATGGATGGTGACAGTGTCGGATAGTGCAGGGAAGTGCAAAGAATTTGTAAATGTTGCAGATGGGGGCGCAGAAACGCCCCCGATTTACCATTAAGCGATCGCGtt  <  1:133454/274‑1 (MQ=255)
                                                                                                                                                                                                                                acTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGtaat                                                                                                                              >  5:263399/1‑148 (MQ=255)
                                                                                                                                                                                                                                  tggccGCTGCCGCCAGCGGGCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGtaataa                                                                                                                            >  5:100954/4‑148 (MQ=255)
                                                                                                                                                                                                                                  tGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGTAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATcc                                                                                                                                                       >  1:49577/1‑121 (MQ=255)
                                                                                                                                                                                                                                  tGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGtaataa                                                                                                                            <  5:365889/148‑1 (MQ=255)
                                                                                                                                                                                                                                   ggccGCTGCCGCCAGCGGGCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTg                                                                                                                                                                        >  6:841049/3‑103 (MQ=255)
                                                                                                                                                                                                                                   ggccGCTGCCGCCAGCGGGCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTg                                                                                                                                                                        <  5:841049/101‑1 (MQ=255)
                                                                                                                                                                                                                                    cccGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACa                                                                                                                         <  7:94713/149‑1 (MQ=255)
                                                                                                                                                                                                                                    cccGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATNGTAATAAACAGATAAGCCTGAAt                                                                                                            >  3:99398/1‑162 (MQ=255)
                                                                                                                                                                                                                                      cGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCataa                                                                                                                                                   <  3:49576/121‑1 (MQ=255)
                                                                                                                                                                                                                                       gCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGAt                                                                                                                      <  7:327415/149‑1 (MQ=255)
                                                                                                                                                                                                                                        ctgccgccAGCGGGCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATa                                                                                                                     <  7:826173/142‑1 (MQ=255)
                                                                                                                                                                                                                                        cTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCTGAATTatg                                                                                                        <  1:99399/162‑1 (MQ=255)
                                                                                                                                                                                                                                         tgccgccAGCGGGCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATa                                                                                                                     <  6:100954/142‑1 (MQ=255)
                                                                                                                                                                                                                                         tGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAAc                                                                                                                          >  4:808389/1‑143 (MQ=255)
                                                                                                                                                                                                                                           gcggcAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAg                                                                                                                   <  6:630346/148‑1 (MQ=255)
                                                                                                                                                                                                                                             ggcAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGAt                                                                                                                      <  2:808391/143‑1 (MQ=255)
                                                                                                                                                                                                                                                    gCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCTGAATTa                                                                                                          >  6:163394/1‑148 (MQ=255)
                                                                                                                                                                                                                                                     ccAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCa                                                                                                                                                      >  3:64515/1‑103 (MQ=255)
                                                                                                                                                                                                                                                     ccAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCTGAATTat                                                                                                         >  6:433547/1‑148 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                            
GACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCAACGGGATCGTTAACACGCTGCCCGCGCGCGGTACGTAAGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCTGAATTATGGATGGTGACAGTGTCGGATAGTGCAGGGAAGTGCAAAGAATTTGTAAATGTTGCAGATGGGGGCGCAGAAACGCCCCCGATTTACCATTAAGCGATCGCGTT  >  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: