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→AGCsltsoluble 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

AATCCCAGCCGGTCGATCACCCTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGG  >  NC_000913/4630594‑4631077
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
aaccccagttggTCGTTCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTa                                                                                                                                                                                                                                               <  2:276682/239‑1 (MQ=255)
aaTCCCAGTTGGTCGATCCCCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTAcg                                                                                                                                                                                                                                             <  1:220361/251‑1 (MQ=255)
 aTCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTAcgc                                                                                                                                                                                                                                            <  1:137423/251‑1 (MQ=255)
              gATCACCTTCGCAGTAAAAGTAACCAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACATGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGGGGCGAGAGCCGACTCCCTGGTTGAGACACGTAGCCGTTACGCGCaattccagccg                                                                                                                                                                                                                               >  2:261614/1‑242 (MQ=255)
                   ccttCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTg                                                                                                                                                                                                                         >  1:374047/4‑251 (MQ=255)
                        gCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGAt                                                                                                                                                                                                                     >  2:98517/1‑250 (MQ=255)
                           gtaaaagtaaaaAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATc                                                                                                                                                                                                                 >  2:237803/1‑251 (MQ=255)
                                          tgtCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATACTCGCCAAATGGAttggg                                                                                                                                                                                                 >  2:104878/1‑247 (MQ=255)
                                             ccGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTc                                                                                                                                                                                               <  1:412175/250‑1 (MQ=255)
                                                 aaCGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAAc                                                                                                                                                                                           >  1:361529/1‑250 (MQ=255)
                                                  aCGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAACCAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAGTCGACAAATGGATGTGGTCGAACa                                                                                                                                                                                          >  2:199320/1‑250 (MQ=255)
                                                      ggcttAAACGGCAATGACTGTTTAGCTTAACTCTATTACGCTGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAATGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGCCTGCTGACGGTCAGCATCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCATGCCTGGGATAATCGACAAATGGATTTGGTCGAACAAatga                                                                                                                                                                                     <  1:439940/247‑1 (MQ=255)
                                                      ggCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAatga                                                                                                                                                                                     <  1:237803/251‑1 (MQ=255)
                                                        cGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATTAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAatgatg                                                                                                                                                                                   <  1:460676/251‑1 (MQ=255)
                                                             aCGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTgg                                                                                                                                                                              >  1:53257/1‑251 (MQ=255)
                                                              cggaactGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGACTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACCAATGGATGTGGTCGACAAAATGATGCATgg                                                                                                                                                                              >  2:155165/7‑250 (MQ=255)
                                                                   ctGACTGGTTAGCATAAATCTATTACGAGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGGGGAAAAAGCCAAAAAAGTTACCTGGCGGCTGTTGGCTGCCGGGGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGGGCCGACTCACTGGAAGAGCAGCGTAGCCGTTACGCGCAAATACAGCAGACCTAGGGAAATCGACAAATGGAGGTGGGCGAACACATGATGcagggagtga                                                                                                                                                                         >  2:109584/2‑241 (MQ=255)
                                                                       cTGGTTAGCATAAATCTATTACGCGGCCTGACACTTCATTGCTGTATTTACACTTAGAGGATGCGCTTGTGGAAAACGCCAAACAAGTTACATGGCGGCTGGTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGAGCGAGACGACCCACTGGATCAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATACTCGCCAAACGGATGTGGTCGAACacattgtacccagactgaa                                                                                                                                                                        >  2:27222/1‑229 (MQ=255)
                                                                        tGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGATGTTGTCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGAGCGGGCCGACTCACTGGATGAGCCGCGTAGCCGATACGCGCAAATCAAGCCGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCAGGACTGaaaaatt                                                                                                                                                                   >  2:360772/1‑246 (MQ=255)
                                                                            tAGCATAAATCTATTAAGCGGCATGACGCTGCATTGATGTATTTACACATAGAGGATGCGCTTGTGAAAAAAGCCAAACAAGTTACCTGGAGGATGTTGGCTGCCGGTGTCTGTCTGATGACGGTAAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGAAAATAAAGCAGGCCTGGGATAATCGAAAAAATGATGTGGTAGAACCAATGTTGCCTGGACGGaaagcatatcc                                                                                                                                                               >  2:348043/1‑242 (MQ=255)
                                                                             aGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATcca                                                                                                                                                              <  1:305138/251‑2 (MQ=255)
                                                                                aTAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACtt                                                                                                                                                           >  1:322880/1‑251 (MQ=255)
                                                                                   aaTCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTAt                                                                                                                                                        <  1:187957/251‑1 (MQ=255)
                                                                                      cTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGGGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTCACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGGCGAACAAATGATGCCTGGACTGAAAGCTTATACACTTTATccc                                                                                                                                                     >  2:179769/1‑251 (MQ=255)
                                                                                      cTATTACGCGGCATGACGCTGCATTGATGAATTTACCCTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGAGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATACTCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGACAGATTATCCACTTTATccc                                                                                                                                                     >  1:204057/1‑251 (MQ=255)
                                                                                              cgGCATGACGCTGCATTGGTGAATTTACACTTGGAGGATACGCTTGTGTAAAAAGCCAAACACGTTACCTAGAGACTGTTGGCTGACGGTGTCTGACTGCTGCAGGTAAGCAGCGTGGCGCGAACCGACTCAGCGGCTGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGAGCTAATCGACTAAATGATGTGGTAGAGCAACTGGTGCAACGAAGGAAAGATTATCCACTTTATCCCTATCTga                                                                                                                                              >  2:55373/1‑249 (MQ=255)
                                                                                                        cTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCgtc                                                                                                                                     >  1:396692/1‑247 (MQ=255)
                                                                                                                            ttAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCCCTGGATGAGCAGCGTAGCCGTTACGCGCAAATCCAGCGGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCGGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCtgatga                                                                                                                >  1:452649/1‑250 (MQ=255)
                                                                                                                                agGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCa                                                                                                            >  2:118450/1‑250 (MQ=255)
                                                                                                                                  gATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAAcc                                                                                                         >  1:471795/1‑251 (MQ=255)
                                                                                                                                      cgcTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACcggcg                                                                                                     <  1:282237/251‑1 (MQ=255)
                                                                                                                                          tgtgGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGt                                                                                                   <  1:98517/249‑1 (MQ=255)
                                                                                                                                           gtgGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCAc                                                                                                                              <  2:377925/221‑1 (MQ=255)
                                                                                                                                           gtgGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCAc                                                                                                                              >  1:377925/1‑221 (MQ=255)
                                                                                                                                            tgGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGAct                                                                                               >  2:153390/1‑250 (MQ=255)
                                                                                                                                                aaaaGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCa                                                                                           <  1:330521/251‑1 (MQ=255)
                                                                                                                                                aaaaGCCAAACAAGATAACTGGCGGCTGTTGGCTGACGGTGTATGTCTGCTGACGGTAAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGAAAATCAAGCAGGACTGGGATAATAGACAACTGGAGGTGGGCGAACAAATGATGACTGGACTGAAAGATTATCCACTTTATCCCTAACTGGAATACCGCCAGATAACCGCTACACTGATGAATAAACCGGCGGGGActttc                                                                                            >  2:448577/1‑246 (MQ=255)
                                                                                                                                                 aaaGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCAc                                                                                          >  1:73368/1‑251 (MQ=255)
                                                                                                                                                  aaGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGCCCTGTTGAATCAACCGGCGGGGACTGTCAcc                                                                                         >  2:322175/1‑250 (MQ=255)
                                                                                                                                                               tACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCATCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAACCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACATATGGATGAGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGTATACCGTCAACTCACCGAGAACCTGATTAATCAAACGACGGTGAccttaatgaaattttttctct                                                                            >  2:364007/1‑231 (MQ=255)
                                                                                                                                                                  cTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCt                                                                          >  1:462772/1‑250 (MQ=255)
                                                                                                                                                                    ggcggcTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAAc                                                                       >  1:269812/1‑251 (MQ=255)
                                                                                                                                                                     ggggcTGTTGGCTGCCGGGGTCTGTCTGTTGACGGTCAGAAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTCCCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAAcc                                                                      <  2:330387/249‑1 (MQ=255)
                                                                                                                                                                       ggcTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAAcc                                                                                                         <  1:230161/214‑1 (MQ=255)
                                                                                                                                                                       ggcTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAAcc                                                                                                         >  2:230161/1‑214 (MQ=255)
                                                                                                                                                                       ggcTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGAGGTGACAGTAACTAACTTTTTTCGCCCTAAcccc                                                                    >  1:107477/1‑251 (MQ=255)
                                                                                                                                                                           gTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACg                                                                 <  2:133864/250‑1 (MQ=255)
                                                                                                                                                                             tGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCCCTAACTTTGTTCTCGCTAACCCCACgatt                                                              >  2:137657/1‑248 (MQ=255)
                                                                                                                                                                             gggCTGCCGGTGTCTGGCGGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCtt                                                              <  2:361529/250‑1 (MQ=255)
                                                                                                                                                                                             tctgCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGc                                              >  1:235764/1‑251 (MQ=255)
                                                                                                                                                                                                    aCGGTCAACAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGACGTGGTCGAACAAATGATGCCTGGACTGAAGGGTTATCCGCTTTATCCCTACCTGGAAAACCGCCAGATCACCGATGATCAGATGAATCAACCGGCGGTGACGGTCAATAACTTTGATCACGCTAACCCCACGCTTCCTCCCGCTCGCCCGCTGCAAtc                                       >  1:111233/1‑251 (MQ=255)
                                                                                                                                                                                                      ggTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAAtctc                                     >  1:155982/1‑251 (MQ=255)
                                                                                                                                                                                                       gTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCg                                    >  1:211951/1‑251 (MQ=255)
                                                                                                                                                                                                              gagTGGCGCGAGCCGAGCCACTGGAGGTGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTTAAAGATTATCCTCTTTACCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACATTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTc                             <  2:506566/249‑1 (MQ=255)
                                                                                                                                                                                                                       gAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGGGACTGGCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTg                    >  2:21120/1‑251 (MQ=255)
                                                                                                                                                                                                                            gACTCACTGGATGAGCAGCGTAGTCGTTATGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGTCGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTccc                                                        >  2:23962/1‑210 (MQ=255)
                                                                                                                                                                                                                            gACTCACTGGATGAGCAGCGTAGTCGTTATGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTccc                                                        <  1:23962/210‑1 (MQ=255)
                                                                                                                                                                                                                            gACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGAGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGAGACCGTCACTAACTTTGTTCGCGCTAACCCCCCGCTTCCTCCCGCCCGCCCGCTGCAATCCCGTTTCGTCAATggaacggcgcg               >  2:462864/1‑241 (MQ=255)
                                                                                                                                                                                                                                         agcagcGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTgg  <  1:199320/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
AATCCCAGCCGGTCGATCACCCTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGG  >  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: