Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,630,921 G→A P63P (CCG→CCA)  slt → soluble lytic murein transglycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,630,9210GA94.7% 154.1 / ‑0.6 61P63P (CCG→CCA) sltsoluble lytic murein transglycosylase
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (29/25);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.32e-01

ATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTT‑AGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCGAAATGGAACACCGGGC  >  NC_000913/4630674‑4631162
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
aTAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACtt                                                                                                                                                                                                                                                   >  1:322880/1‑251 (MQ=255)
   aaTCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTAt                                                                                                                                                                                                                                                <  1:187957/251‑1 (MQ=255)
      cTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGGGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTCACGCGCAAATCAAGCAGGCCTGGGATAATCGAC‑AAATGGATGTGGGCGAACAAATGATGCCTGGACTGAAAGCTTATACACTTTATccc                                                                                                                                                                                                                                             >  2:179769/1‑251 (MQ=255)
      cTATTACGCGGCATGACGCTGCATTGATGAATTTACCCTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGAGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATACTCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGACAGATTATCCACTTTATccc                                                                                                                                                                                                                                             >  1:204057/1‑251 (MQ=255)
              cgGCATGACGCTGCATTGGTGAATTTACACTTGGAGGATACGCTTGTGTAAAAAGCCAAACACGTTACCTAGAGACTGTTGGCTGACGGTGTCTGACTGCTGCAGGTAAGCAGCGTGGCGCGAACCGACTCAGCGGCTGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGAGCTAATCGACTAAATG‑ATGTGGTAGAGCAACTGGTGCAACGAAGGAAAGATTATCCACTTTATCCCTATCTga                                                                                                                                                                                                                                      >  2:55373/1‑249 (MQ=255)
                        cTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCgtc                                                                                                                                                                                                                             >  1:396692/1‑247 (MQ=255)
                                            ttAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCCCTGGATGAGCAGCGTAGCCGTTACGCGCAAATCCAGCGGGCCTGGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCGGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCtgatga                                                                                                                                                                                                        >  1:452649/1‑250 (MQ=255)
                                                agGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCa                                                                                                                                                                                                    >  2:118450/1‑250 (MQ=255)
                                                  gATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAAcc                                                                                                                                                                                                 >  1:471795/1‑251 (MQ=255)
                                                      cgcTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACcggcg                                                                                                                                                                                             <  1:282237/251‑1 (MQ=255)
                                                          tgtgGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGt                                                                                                                                                                                           <  1:98517/249‑1 (MQ=255)
                                                           gtgGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCAc                                                                                                                                                                                                                      >  1:377925/1‑221 (MQ=255)
                                                           gtgGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCAc                                                                                                                                                                                                                      <  2:377925/221‑1 (MQ=255)
                                                            tgGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGAct                                                                                                                                                                                       >  2:153390/1‑250 (MQ=255)
                                                                aaaaGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCa                                                                                                                                                                                   <  1:330521/251‑1 (MQ=255)
                                                                aaaaGCCAAACAAGATAACTGGCGGCTGTTGGCTGACGGTGTATGTCTGCTGACGGTAAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGAAAATCAAGCAGGACTGGGATAATAGAC‑AACTGGAGGTGGGCGAACAAATGATGACTGGACTGAAAGATTATCCACTTTATCCCTAACTGGAATACCGCCAGATAACCGCTACACTGATGAATAAACCGGCGGGGActttc                                                                                                                                                                                    >  2:448577/1‑246 (MQ=255)
                                                                 aaaGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCAc                                                                                                                                                                                  >  1:73368/1‑251 (MQ=255)
                                                                  aaGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGCCCTGTTGAATCAACCGGCGGGGACTGTCAcc                                                                                                                                                                                 >  2:322175/1‑250 (MQ=255)
                                                                               tACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCATCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAACCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGAC‑ATATGGATGAGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGTATACCGTCAACTCACCGAGAACCTGATTAATCAAACGACGGTGAccttaatgaaattttttctct                                                                                                                                                                    >  2:364007/1‑231 (MQ=255)
                                                                                  cTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCt                                                                                                                                                                  >  1:462772/1‑250 (MQ=255)
                                                                                    ggcggcTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAAc                                                                                                                                                               >  1:269812/1‑251 (MQ=255)
                                                                                     ggggcTGTTGGCTGCCGGGGTCTGTCTGTTGACGGTCAGAAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTCCCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAAcc                                                                                                                                                              <  2:330387/249‑1 (MQ=255)
                                                                                       ggcTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAAcc                                                                                                                                                                                                 <  1:230161/214‑1 (MQ=255)
                                                                                       ggcTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAAcc                                                                                                                                                                                                 >  2:230161/1‑214 (MQ=255)
                                                                                       ggcTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGAGGTGACAGTAACTAACTTTTTTCGCCCTAAcccc                                                                                                                                                            >  1:107477/1‑251 (MQ=255)
                                                                                           gTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACg                                                                                                                                                         <  2:133864/250‑1 (MQ=255)
                                                                                             tGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCCCTAACTTTGTTCTCGCTAACCCCACgatt                                                                                                                                                      >  2:137657/1‑248 (MQ=255)
                                                                                             gggCTGCCGGTGTCTGGCGGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCtt                                                                                                                                                      <  2:361529/250‑1 (MQ=255)
                                                                                                             tctgCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGc                                                                                                                                      >  1:235764/1‑251 (MQ=255)
                                                                                                                    aCGGTCAACAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGAC‑AAATGGACGTGGTCGAACAAATGATGCCTGGACTGAAGGGTTATCCGCTTTATCCCTACCTGGAAAACCGCCAGATCACCGATGATCAGATGAATCAACCGGCGGTGACGGTCAATAACTTTGATCACGCTAACCCCACGCTTCCTCCCGCTCGCCCGCTGCAAtc                                                                                                                               >  1:111233/1‑251 (MQ=255)
                                                                                                                      ggTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAAtctc                                                                                                                             >  1:155982/1‑251 (MQ=255)
                                                                                                                       gTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCg                                                                                                                            >  1:211951/1‑251 (MQ=255)
                                                                                                                              gagTGGCGCGAGCCGAGCCACTGGAGGTGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTTAAAGATTATCCTCTTTACCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACATTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTc                                                                                                                     <  2:506566/249‑1 (MQ=255)
                                                                                                                                       gAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGGGACTGGCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTg                                                                                                            >  2:21120/1‑251 (MQ=255)
                                                                                                                                            gACTCACTGGATGAGCAGCGTAGTCGTTATGCGCAAATCAAGCAGGCCTGGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGTCGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTccc                                                                                                                                                >  2:23962/1‑210 (MQ=255)
                                                                                                                                            gACTCACTGGATGAGCAGCGTAGTCGTTATGCGCAAATCAAGCAGGCCTGGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTccc                                                                                                                                                <  1:23962/210‑1 (MQ=255)
                                                                                                                                            gACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGAC‑AAATGGATGAGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGAGACCGTCACTAACTTTGTTCGCGCTAACCCCCCGCTTCCTCCCGCCCGCCCGCTGCAATCCCGTTTCGTCAATggaacggcgcg                                                                                                       >  2:462864/1‑241 (MQ=255)
                                                                                                                                                         agcagcGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTgg                                                                                          <  1:199320/251‑1 (MQ=255)
                                                                                                                                                                  gccGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGGCTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTg                                                                                 <  2:374047/249‑1 (MQ=255)
                                                                                                                                                                           cgcAAATCAAGCAGGCCTGGGATAAGCGAC‑AAATGGATGTGTTCGAACAAATGATGCCTGGACTGAAAGATTATCCCCTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAAGCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTT‑AGCCttt                                                                       <  2:281841/251‑1 (MQ=255)
                                                                                                                                                                            gcAAATCAAGCAGGCCTGGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTT‑AGCCTTTa                                                                      <  1:166867/251‑1 (MQ=255)
                                                                                                                                                                                      cAGGCCTGGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTT‑AGCCTTTAGCCCGGaaaa                                                            <  1:153390/251‑1 (MQ=255)
                                                                                                                                                                                         gCCTGGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTT‑AGCCTTTAGCCCGGAAAA‑Gcc                                                        <  1:420969/251‑1 (MQ=255)
                                                                                                                                                                                          ccTGGGATAATCGAC‑AAATGGATGTGGTCGAACAAACGATGCCTGGACTGAAAGATTATCCACTTTAGCCCTCCCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCTCGCTGCCTCCCGCTCGCACGCTGCAATCGCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTT‑AGCCTTTAGCCCGGAAAA‑Gccc                                                       <  2:204057/251‑1 (MQ=255)
                                                                                                                                                                                            gtggagaaTCGAC‑AAATGGATGTGGTCGAACAAAGGATGCCTGTACTGAAAGATTAGCCACTTTATCCCTACCTGTAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTGCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTT‑AGCCTTTAGCCCGGAAAA‑GCCCgg                                                     <  2:258358/245‑1 (MQ=255)
                                                                                                                                                                                                aTAATCGAC‑AAATGGATGTGGTCGAACAACTGATGCCTGGACTGAAAGATTATCCACTTTCTCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTTACTTTGTTCGAGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCACTCTCGTTTCGTCCATGAACTGGCGCGGCGGGAAGGCTGGCGTGGCTTGTTTGGCCTTT‑TCCCGGGAAA‑GCCCGGagcc                                                 >  2:493050/1‑248 (MQ=255)
                                                                                                                                                                                                aTAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTT‑AGCCTTTAGCCCGGAAAA‑GCCCGGAACt                                                 <  1:493050/251‑1 (MQ=255)
                                                                                                                                                                                                     cGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTT‑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)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
ATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTT‑AGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCGAAATGGAACACCGGGC  >  NC_000913/4630674‑4631162

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: ‑