Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,630,960 T→C D76D (GAT→GAC)  slt → soluble lytic murein transglycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,630,9600TC94.9% 176.2 / ‑1.6 60D76D (GAT→GAC) sltsoluble lytic murein transglycosylase
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (24/32);  total (27/33)
Fisher's exact test for biased strand distribution p-value = 5.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.79e-01

TTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTT‑AGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCGAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTG  >  NC_000913/4630718‑4631197
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                
ttAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCCCTGGATGAGCAGCGTAGCCGTTACGCGCAAATCCAGCGGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCGGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCtgatga                                                                                                                                                                                                                                           >  1:452649/1‑250 (MQ=255)
    agGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCa                                                                                                                                                                                                                                       >  2:118450/1‑250 (MQ=255)
      gATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAAcc                                                                                                                                                                                                                                    >  1:471795/1‑251 (MQ=255)
          cgcTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACcggcg                                                                                                                                                                                                                                <  1:282237/251‑1 (MQ=255)
              tgtgGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGt                                                                                                                                                                                                                              <  1:98517/249‑1 (MQ=255)
                tgGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGAct                                                                                                                                                                                                                          >  2:153390/1‑250 (MQ=255)
                    aaaaGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCa                                                                                                                                                                                                                      <  1:330521/251‑1 (MQ=255)
                    aaaaGCCAAACAAGATAACTGGCGGCTGTTGGCTGACGGTGTATGTCTGCTGACGGTAAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGAAAATCAAGCAGGACTGGGATAATAGACAACTGGAGGTGGGCGAACAAATGATGACTGGACTGAAAGATTATCCACTTTATCCCTAACTGGAATACCGCCAGATAACCGCTACACTGATGAATAAACCGGCGGGGActttc                                                                                                                                                                                                                       >  2:448577/1‑246 (MQ=255)
                     aaaGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCAc                                                                                                                                                                                                                     >  1:73368/1‑251 (MQ=255)
                      aaGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGCCCTGTTGAATCAACCGGCGGGGACTGTCAcc                                                                                                                                                                                                                    >  2:322175/1‑250 (MQ=255)
                                   tACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCATCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAACCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACATATGGATGAGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGTATACCGTCAACTCACCGAGAACCTGATTAATCAAACGACGGTGAccttaatgaaattttttctct                                                                                                                                                                                                       >  2:364007/1‑231 (MQ=255)
                                      cTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCt                                                                                                                                                                                                     >  1:462772/1‑250 (MQ=255)
                                        ggcggcTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAAc                                                                                                                                                                                                  >  1:269812/1‑251 (MQ=255)
                                         ggggcTGTTGGCTGCCGGGGTCTGTCTGTTGACGGTCAGAAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTCCCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAAcc                                                                                                                                                                                                 <  2:330387/249‑1 (MQ=255)
                                           ggcTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAAcc                                                                                                                                                                                                                                    <  1:230161/214‑1 (MQ=255)
                                           ggcTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAAcc                                                                                                                                                                                                                                    >  2:230161/1‑214 (MQ=255)
                                           ggcTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGAGGTGACAGTAACTAACTTTTTTCGCCCTAAcccc                                                                                                                                                                                               >  1:107477/1‑251 (MQ=255)
                                               gTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACg                                                                                                                                                                                            <  2:133864/250‑1 (MQ=255)
                                                 tGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCCCTAACTTTGTTCTCGCTAACCCCACgatt                                                                                                                                                                                         >  2:137657/1‑248 (MQ=255)
                                                 gggCTGCCGGTGTCTGGCGGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCtt                                                                                                                                                                                         <  2:361529/250‑1 (MQ=255)
                                                                 tctgCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGc                                                                                                                                                                         >  1:235764/1‑251 (MQ=255)
                                                                        aCGGTCAACAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGACGTGGTCGAACAAATGATGCCTGGACTGAAGGGTTATCCGCTTTATCCCTACCTGGAAAACCGCCAGATCACCGATGATCAGATGAATCAACCGGCGGTGACGGTCAATAACTTTGATCACGCTAACCCCACGCTTCCTCCCGCTCGCCCGCTGCAAtc                                                                                                                                                                  >  1:111233/1‑251 (MQ=255)
                                                                          ggTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAAtctc                                                                                                                                                                >  1:155982/1‑251 (MQ=255)
                                                                           gTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCg                                                                                                                                                               >  1:211951/1‑251 (MQ=255)
                                                                                  gagTGGCGCGAGCCGAGCCACTGGAGGTGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTTAAAGATTATCCTCTTTACCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACATTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTc                                                                                                                                                        <  2:506566/249‑1 (MQ=255)
                                                                                           gAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGGGACTGGCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTg                                                                                                                                               >  2:21120/1‑251 (MQ=255)
                                                                                                gACTCACTGGATGAGCAGCGTAGTCGTTATGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGTCGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTccc                                                                                                                                                                                   >  2:23962/1‑210 (MQ=255)
                                                                                                gACTCACTGGATGAGCAGCGTAGTCGTTATGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTccc                                                                                                                                                                                   <  1:23962/210‑1 (MQ=255)
                                                                                                gACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGAGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGAGACCGTCACTAACTTTGTTCGCGCTAACCCCCCGCTTCCTCCCGCCCGCCCGCTGCAATCCCGTTTCGTCAATggaacggcgcg                                                                                                                                          >  2:462864/1‑241 (MQ=255)
                                                                                                             agcagcGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTgg                                                                                                                             <  1:199320/251‑1 (MQ=255)
                                                                                                                      gccGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGGCTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTg                                                                                                                    <  2:374047/249‑1 (MQ=255)
                                                                                                                               cgcAAATCAAGCAGGCCTGGGATAAGCGACAAATGGATGTGTTCGAACAAATGATGCCTGGACTGAAAGATTATCCCCTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAAGCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTT‑AGCCttt                                                                                                          <  2:281841/251‑1 (MQ=255)
                                                                                                                                gcAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTT‑AGCCTTTa                                                                                                         <  1:166867/251‑1 (MQ=255)
                                                                                                                                          cAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTT‑AGCCTTTAGCCCGGaaaa                                                                                               <  1:153390/251‑1 (MQ=255)
                                                                                                                                             gCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTT‑AGCCTTTAGCCCGGAAAA‑Gcc                                                                                           <  1:420969/251‑1 (MQ=255)
                                                                                                                                              ccTGGGATAATCGACAAATGGATGTGGTCGAACAAACGATGCCTGGACTGAAAGATTATCCACTTTAGCCCTCCCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCTCGCTGCCTCCCGCTCGCACGCTGCAATCGCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTT‑AGCCTTTAGCCCGGAAAA‑Gccc                                                                                          <  2:204057/251‑1 (MQ=255)
                                                                                                                                                gtggagaaTCGACAAATGGATGTGGTCGAACAAAGGATGCCTGTACTGAAAGATTAGCCACTTTATCCCTACCTGTAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTGCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTT‑AGCCTTTAGCCCGGAAAA‑GCCCgg                                                                                        <  2:258358/245‑1 (MQ=255)
                                                                                                                                                    aTAATCGACAAATGGATGTGGTCGAACAACTGATGCCTGGACTGAAAGATTATCCACTTTCTCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTTACTTTGTTCGAGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCACTCTCGTTTCGTCCATGAACTGGCGCGGCGGGAAGGCTGGCGTGGCTTGTTTGGCCTTT‑TCCCGGGAAA‑GCCCGGagcc                                                                                    >  2:493050/1‑248 (MQ=255)
                                                                                                                                                    aTAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTT‑AGCCTTTAGCCCGGAAAA‑GCCCGGAACt                                                                                    <  1:493050/251‑1 (MQ=255)
                                                                                                                                                         cGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTT‑AGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCg                                                                                <  1:137657/250‑1 (MQ=255)
                                                                                                                                                             aaaTGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGGATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTT‑AGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAgcg                                                                           <  2:172334/251‑1 (MQ=255)
                                                                                                                                                                          cGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCGGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTT‑AGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATtact                                                             <  2:53257/251‑1 (MQ=255)
                                                                                                                                                                          cGAACAAAAGATGCCTGGACTGAAAGATTATCCCCTTTATCCCTTCCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTT‑AGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATtact                                                             <  2:267482/251‑1 (MQ=255)
                                                                                                                                                                                     tgCCAGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTATCCCCACGCTTCCCCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTT‑AGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGcaaa                                                   <  2:442337/250‑4 (MQ=255)
                                                                                                                                                                                     ggCCTGGACTGAAAGATTAGCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCTCTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACGTGCGCGGCGTGAATACTGGCGTGGCTTGTT‑AGTCTTGAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCaaaa                                                  <  1:62715/250‑1 (MQ=255)
                                                                                                                                                                                       ccTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTT‑AGCCTTTAGCCCGGAAAA‑ACCCGGAACTACCGGAGCGCCAAT‑TAATTACTACTATGCAAAATg                                                >  1:63547/1‑251 (MQ=255)
                                                                                                                                                                                        cTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGAGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTT‑TGGCTTTAGCCCGGAAAA‑GGCCAGCACTACCGGAGCGC‑ACTGTAATTACTACTATGCCAAATgg                                               >  2:40670/1‑251 (MQ=255)
                                                                                                                                                                                           gACTGAAAGATTATCTACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTACCCCCACGCTTCCGCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGTCTGGCGTGGCTTGTT‑AGCCTTTAGCTCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCAAAATGGAac                                            <  2:455632/251‑1 (MQ=255)
                                                                                                                                                                                           gACTGAAAGATTATCCCCTTTATCCCTACCTGGAATTCCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGGCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTGCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTT‑AGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCAAAAGGGAac                                            <  2:455641/251‑1 (MQ=255)
                                                                                                                                                                                           gACTGAAAGATTATCCACTTTATCCCTACCTTGAATACCGTCAGATCCCCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACCTTGTTCGCGCTAACCCCACGCTTCCGCTCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTT‑AGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCAAAATGGAac                                            <  1:364007/251‑1 (MQ=255)
                                                                                                                                                                                           gACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTT‑AGCCTTTTGCCCGGAAAAACCCCGGAA‑TACCGAAGCGC‑AATGTAATTACTACTATGCAAAATGGAac                                            >  2:177446/1‑251 (MQ=255)
                                                                                                                                                                                              ggaaagATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGGCACTAACTTTGTTCGCGCTAACCCCACGCTGCCTCCCGCTCGCACGCTGCAATCTCGTTGCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGTCTTGTT‑AGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCAAAATGGAACAcc                                         <  2:471795/246‑1 (MQ=255)
                                                                                                                                                                                                    aTTATCCACTTTATCCCTACCTGGCATACAGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGACCTGGCGCGGCGAGAAGACTGGCGTGGCTTGTT‑AGACTTTAGACCGGAAAA‑GCCCGGGACTACCGAAGCGG‑AAAGTAATTAATACTATGCAAAATGGAAAACCgggg                                     >  1:357877/1‑248 (MQ=255)
                                                                                                                                                                                                              ttATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTT‑AGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTgaagaa                          <  1:262683/250‑1 (MQ=255)
                                                                                                                                                                                                               tATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTT‑AGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGaagaag                         <  1:21120/250‑1 (MQ=255)
                                                                                                                                                                                                                 cccctccctGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTT‑AGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACGACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCt                      <  2:452649/250‑1 (MQ=255)
                                                                                                                                                                                                                        ccggaaacacttcAGATCACCTATGACCTGAGGAACCAACCGGCGGTGCCTGTCACTAACTTTGTTCGCGCTAACCCACCGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATTAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTT‑AGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAgg               <  2:119721/239‑1 (MQ=255)
                                                                                                                                                                                                                                    cAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTT‑AGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATg                                                                     <  2:477644/185‑1 (MQ=255)
                                                                                                                                                                                                                                    cAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTT‑AGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATg                                                                     >  1:477644/1‑185 (MQ=255)
                                                                                                                                                                                                                                     aGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTT‑TGCCTTTAGCCCGGGAAA‑ACCCGGAACTACCGAAGCGC‑AATGGAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAAAAGCCTGGCCAGGGGCGAAagaggtg  >  1:270689/1‑248 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                
TTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTT‑AGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCGAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTG  >  NC_000913/4630718‑4631197

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