Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,630,837 T→C S35S (AGT→AGC)  slt → soluble lytic murein transglycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,630,8370TC93.2% 162.7 / ‑1.9 59S35S (AGT→AGC) sltsoluble lytic murein transglycosylase
Reads supporting (aligned to +/- strand):  ref base T (3/1);  new base C (36/19);  total (39/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

AATCCCAGCCGGTCGATCACCCTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCT‑GGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGG  >  NC_000913/4630594‑4631077
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
aaccccagttggTCGTTCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTTa                                                                                                                                                                                                                                               <  2:276682/239‑1 (MQ=255)
aaTCCCAGTTGGTCGATCCCCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTTAcg                                                                                                                                                                                                                                             <  1:220361/251‑1 (MQ=255)
 aTCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTTAcgc                                                                                                                                                                                                                                            <  1:137423/251‑1 (MQ=255)
              gATCACCTTCGCAGTAAAAGTAACCAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACATGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGGGGCGAGAGCCGACTCCCTGGTTGAGACA‑CGTAGCCGTTACGCGCaattccagccg                                                                                                                                                                                                                               >  2:261614/1‑242 (MQ=255)
                   ccttCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCT‑g                                                                                                                                                                                                                         >  1:374047/4‑251 (MQ=255)
                        gCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCT‑GGGAt                                                                                                                                                                                                                     >  2:98517/1‑250 (MQ=255)
                           gtaaaagtaaaaAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCT‑GGGATAATc                                                                                                                                                                                                                 >  2:237803/1‑251 (MQ=255)
                                          tgtCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCT‑GGGATACTCGCC‑AAATGGAttggg                                                                                                                                                                                                 >  2:104878/1‑247 (MQ=255)
                                             ccGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCT‑GGGATAATCGAC‑AAATGGATGTGGTc                                                                                                                                                                                               <  1:412175/250‑1 (MQ=255)
                                                 aaCGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCT‑GGGATAATCGAC‑AAATGGATGTGGTCGAAc                                                                                                                                                                                           >  1:361529/1‑250 (MQ=255)
                                                  aCGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAACCAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCT‑GGGATAGTCGAC‑AAATGGATGTGGTCGAACa                                                                                                                                                                                          >  2:199320/1‑250 (MQ=255)
                                                      ggcttAAACGGCAATGACTGTTTAGCTTAACTCTATTACGCTGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAATGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGCCTGCTGACGGTCAGCATCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTTACGCGCAAATCAAGCATGCCT‑GGGATAATCGAC‑AAATGGATTTGGTCGAACAAatga                                                                                                                                                                                     <  1:439940/247‑1 (MQ=255)
                                                      ggCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCT‑GGGATAATCGAC‑AAATGGATGTGGTCGAACAAatga                                                                                                                                                                                     <  1:237803/251‑1 (MQ=255)
                                                        cGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATTAG‑CAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCT‑GGGATAATCGAC‑AAATGGATGTGGTCGAACAAatgatg                                                                                                                                                                                   <  1:460676/251‑1 (MQ=255)
                                                             aCGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCT‑GGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTgg                                                                                                                                                                              >  1:53257/1‑251 (MQ=255)
                                                              cggaactGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGACTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCT‑GGGATAATCGAC‑CAATGGATGTGGTCGACAAAATGATGCATgg                                                                                                                                                                              >  2:155165/7‑250 (MQ=255)
                                                                   ctGACTGGTTAGCATAAATCTATTACGAGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGGGGAAAAAGCCAAAAAAGTTACCTGGCGGCTGTTGGCTGCCGGGGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGGGCCGACTCACTGGAAGAG‑CAGCGTAGCCGTTACGCGCAAATACAGCAGACCTAGGGA‑AATCGAC‑AAATGGAGGTGGGCGAACACATGATGcagggagtga                                                                                                                                                                         >  2:109584/2‑241 (MQ=255)
                                                                       cTGGTTAGCATAAATCTATTACGCGGCCTGACACTTCATTGCTGTATTTACACTTAGAGGATGCGCTTGTGGAAAACGCCAAACAAGTTACATGGCGGCTGGTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGAGCGAGACGACCCACTGGATCAG‑CAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCT‑GGGATACTCGCC‑AAACGGATGTGGTCGAACacattgtacccagactgaa                                                                                                                                                                        >  2:27222/1‑229 (MQ=255)
                                                                        tGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGATGTTGTCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGAGCGGGCCGACTCACTGGATGAG‑CCGCGTAGCCGATACGCGCAAATCAAGCCGGCCT‑GGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCAGGACTGaaaaatt                                                                                                                                                                   >  2:360772/1‑246 (MQ=255)
                                                                            tAGCATAAATCTATTAAGCGGCATGACGCTGCATTGATGTATTTACACATAGAGGATGCGCTTGTGAAAAAAGCCAAACAAGTTACCTGGAGGATGTTGGCTGCCGGTGTCTGTCTGATGACGGTAAGCAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTTACGCGAAAATAAAGCAGGCCT‑GGGATAATCGAA‑AAAATGATGTGGTAGAACCAATGTTGCCTGGACGGaaagcatatcc                                                                                                                                                               >  2:348043/1‑242 (MQ=255)
                                                                             aGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCT‑GGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATcca                                                                                                                                                              <  1:305138/251‑2 (MQ=255)
                                                                                aTAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCT‑GGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACtt                                                                                                                                                           >  1:322880/1‑251 (MQ=255)
                                                                                   aaTCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCT‑GGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTAt                                                                                                                                                        <  1:187957/251‑1 (MQ=255)
                                                                                      cTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGGGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTCACGCGCAAATCAAGCAGGCCT‑GGGATAATCGAC‑AAATGGATGTGGGCGAACAAATGATGCCTGGACTGAAAGCTTATACACTTTATccc                                                                                                                                                     >  2:179769/1‑251 (MQ=255)
                                                                                      cTATTACGCGGCATGACGCTGCATTGATGAATTTACCCTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGAGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCT‑GGGATACTCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGACAGATTATCCACTTTATccc                                                                                                                                                     >  1:204057/1‑251 (MQ=255)
                                                                                              cgGCATGACGCTGCATTGGTGAATTTACACTTGGAGGATACGCTTGTGTAAAAAGCCAAACACGTTACCTAGAGACTGTTGGCTGACGGTGTCTGACTGCTGCAGGTAAGCAGCGTGGCGCGAACCGACTCAGCGGCTGAG‑CAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCT‑GAGCTAATCGACTAAATG‑ATGTGGTAGAGCAACTGGTGCAACGAAGGAAAGATTATCCACTTTATCCCTATCTga                                                                                                                                              >  2:55373/1‑249 (MQ=255)
                                                                                                        cTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCT‑GGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCgtc                                                                                                                                     >  1:396692/1‑247 (MQ=255)
                                                                                                                            ttAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCCCTGGATGAG‑CAGCGTAGCCGTTACGCGCAAATCCAGCGGGCCT‑GGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCGGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCtgatga                                                                                                                >  1:452649/1‑250 (MQ=255)
                                                                                                                                agGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCT‑GGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCa                                                                                                            >  2:118450/1‑250 (MQ=255)
                                                                                                                                  gATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCT‑GGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAAcc                                                                                                         >  1:471795/1‑251 (MQ=255)
                                                                                                                                      cgcTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCT‑GGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACcggcg                                                                                                     <  1:282237/251‑1 (MQ=255)
                                                                                                                                          tgtgGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCT‑GGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGt                                                                                                   <  1:98517/249‑1 (MQ=255)
                                                                                                                                           gtgGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCT‑GGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCAc                                                                                                                              <  2:377925/221‑1 (MQ=255)
                                                                                                                                           gtgGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCT‑GGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCAc                                                                                                                              >  1:377925/1‑221 (MQ=255)
                                                                                                                                            tgGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCT‑GGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGAct                                                                                               >  2:153390/1‑250 (MQ=255)
                                                                                                                                                aaaaGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCT‑GGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCa                                                                                           <  1:330521/251‑1 (MQ=255)
                                                                                                                                                aaaaGCCAAACAAGATAACTGGCGGCTGTTGGCTGACGGTGTATGTCTGCTGACGGTAAGCAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTTACGCGAAAATCAAGCAGGACT‑GGGATAATAGAC‑AACTGGAGGTGGGCGAACAAATGATGACTGGACTGAAAGATTATCCACTTTATCCCTAACTGGAATACCGCCAGATAACCGCTACACTGATGAATAAACCGGCGGGGActttc                                                                                            >  2:448577/1‑246 (MQ=255)
                                                                                                                                                 aaaGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCT‑GGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCAc                                                                                          >  1:73368/1‑251 (MQ=255)
                                                                                                                                                  aaGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCT‑GGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGCCCTGTTGAATCAACCGGCGGGGACTGTCAcc                                                                                         >  2:322175/1‑250 (MQ=255)
                                                                                                                                                               tACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCATCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAACCGTTACGCGCAAATCAAGCAGGCCT‑GGGATAATCGAC‑ATATGGATGAGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGTATACCGTCAACTCACCGAGAACCTGATTAATCAAACGACGGTGAccttaatgaaattttttctct                                                                            >  2:364007/1‑231 (MQ=255)
                                                                                                                                                                  cTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCT‑GGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCt                                                                          >  1:462772/1‑250 (MQ=255)
                                                                                                                                                                    ggcggcTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCT‑GGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAAc                                                                       >  1:269812/1‑251 (MQ=255)
                                                                                                                                                                     ggggcTGTTGGCTGCCGGGGTCTGTCTGTTGACGGTCAGAAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCT‑GGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTCCCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAAcc                                                                      <  2:330387/249‑1 (MQ=255)
                                                                                                                                                                       ggcTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCT‑GGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAAcc                                                                                                         <  1:230161/214‑1 (MQ=255)
                                                                                                                                                                       ggcTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCT‑GGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAAcc                                                                                                         >  2:230161/1‑214 (MQ=255)
                                                                                                                                                                       ggcTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCT‑GGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGAGGTGACAGTAACTAACTTTTTTCGCCCTAAcccc                                                                    >  1:107477/1‑251 (MQ=255)
                                                                                                                                                                           gTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCT‑GGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACg                                                                 <  2:133864/250‑1 (MQ=255)
                                                                                                                                                                             tGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCT‑GGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCCCTAACTTTGTTCTCGCTAACCCCACgatt                                                              >  2:137657/1‑248 (MQ=255)
                                                                                                                                                                             gggCTGCCGGTGTCTGGCGGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCT‑GGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCtt                                                              <  2:361529/250‑1 (MQ=255)
                                                                                                                                                                                             tctgCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCT‑GGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGc                                              >  1:235764/1‑251 (MQ=255)
                                                                                                                                                                                                    aCGGTCAACAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCT‑GGGATAATCGAC‑AAATGGACGTGGTCGAACAAATGATGCCTGGACTGAAGGGTTATCCGCTTTATCCCTACCTGGAAAACCGCCAGATCACCGATGATCAGATGAATCAACCGGCGGTGACGGTCAATAACTTTGATCACGCTAACCCCACGCTTCCTCCCGCTCGCCCGCTGCAAtc                                       >  1:111233/1‑251 (MQ=255)
                                                                                                                                                                                                      ggTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCT‑GGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAAtctc                                     >  1:155982/1‑251 (MQ=255)
                                                                                                                                                                                                       gTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCT‑GGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCg                                    >  1:211951/1‑251 (MQ=255)
                                                                                                                                                                                                              gagTGGCGCGAGCCGAGCCACTGGAGGTG‑CAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCT‑GGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTTAAAGATTATCCTCTTTACCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACATTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTc                             <  2:506566/249‑1 (MQ=255)
                                                                                                                                                                                                                       gAGCCGACTCACTGGATGAG‑CAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCT‑GGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGGGACTGGCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTg                    >  2:21120/1‑251 (MQ=255)
                                                                                                                                                                                                                            gACTCACTGGATGAG‑CAGCGTAGTCGTTATGCGCAAATCAAGCAGGCCT‑GGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGTCGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTccc                                                        >  2:23962/1‑210 (MQ=255)
                                                                                                                                                                                                                            gACTCACTGGATGAG‑CAGCGTAGTCGTTATGCGCAAATCAAGCAGGCCT‑GGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTccc                                                        <  1:23962/210‑1 (MQ=255)
                                                                                                                                                                                                                            gACTCACTGGATGAG‑CAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCT‑GGGATAATCGAC‑AAATGGATGAGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGAGACCGTCACTAACTTTGTTCGCGCTAACCCCCCGCTTCCTCCCGCCCGCCCGCTGCAATCCCGTTTCGTCAATggaacggcgcg               >  2:462864/1‑241 (MQ=255)
                                                                                                                                                                                                                                         ag‑cagcGTAGCCGTTACGCGCAAATCAAGCAGGCCT‑GGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTgg  <  1:199320/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
AATCCCAGCCGGTCGATCACCCTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAG‑CAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCT‑GGGATAATCGAC‑AAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGG  >  NC_000913/4630594‑4631077

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