Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,630,912 G→A K60K (AAG→AAA)  slt → soluble lytic murein transglycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,630,9120GA95.0% 164.5 / ‑0.9 60K60K (AAG→AAA) sltsoluble lytic murein transglycosylase
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (28/29);  total (30/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 = 8.39e-01

AGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTT‑AGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCGAAATGGAAC  >  NC_000913/4630671‑4631155
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      
aGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATcca                                                                                                                                                                                                                                               <  1:305138/251‑2 (MQ=255)
   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                                                                                                                                                                                                               <  2:377925/221‑1 (MQ=255)
                                                              gtgGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCAc                                                                                                                                                                                                               >  1:377925/1‑221 (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)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      
AGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTT‑AGCCTTTAGCCCGGAAAA‑GCCCGGAACTACCGAAGCGC‑AATGTAATTACTACTATGCGAAATGGAAC  >  NC_000913/4630671‑4631155

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