Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,249,684 C→T A44V (GCT→GTT)  malM → maltose regulon periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,249,6840CT98.2% 171.3 / NA 58A44V (GCT→GTT) malMmaltose regulon periplasmic protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (33/23);  minor base G (1/0);  total (34/24)
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

GCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGCTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAACGTGCTGATTCTTGATCAGAACATGACCCCATCAGCCTTCTTCCCC  >  NC_000913/4249437‑4249928
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
gcATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTGACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCtgttgc                                                                                                                                                                                                                                                   <  1:38746/251‑1 (MQ=255)
 cATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTgttgcg                                                                                                                                                                                                                                                  <  1:349660/251‑1 (MQ=255)
      ggCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGca                                                                                                                                                                                                                                             <  1:442289/251‑1 (MQ=255)
               aCAGGATGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACCATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGCTGTTAACTACGTACCGCAAAACACCAGCGAAGCGCCAGCCAATACATATGTTGCGCTGCAAAAACTCAc                                                                                                                                                                                                                                    >  2:412150/1‑251 (MQ=255)
                cAGGTTGTCAGTATCTGAAAGGGGCGAAAGCCCCGCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACGCAcc                                                                                                                                                                                                                                   <  2:56252/251‑1 (MQ=255)
                   gTTGTCATTATCTGAAATGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTgg                                                                                                                                                                                                                                <  1:80615/251‑1 (MQ=255)
                       tCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTAACACAACTCACCTGGacac                                                                                                                                                                                                                            >  2:274601/1‑251 (MQ=255)
                         aTTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCg                                                                                                                                                                                                                          <  2:119256/251‑1 (MQ=255)
                                          aaaGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCACTCTAAAAc                                                                                                                                                                                                         >  2:306619/1‑251 (MQ=255)
                                                      aTTATCTGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCGTCGTCATCTGTTTATCAGCAGGGTTACTGGCAAGAGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGAACACCGGTCGATCAATCTACAACCCATAcccccc                                                                                                                                                                                              >  2:304958/1‑250 (MQ=255)
                                                            ggggttAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGc                                                                                                                                                                                       <  2:474793/247‑1 (MQ=255)
                                                                  aGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACcgg                                                                                                                                                                                 <  1:568911/251‑1 (MQ=255)
                                                                      cgcaaTTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGTTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCATGTTTACTGGCAAGCGCGCCGGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCGGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGTCCggcggc                                                                                                                                                                             <  2:548467/247‑1 (MQ=255)
                                                                        cTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGCAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCAATTCCATCTGGTGCGCTGCAACAACTCACCTGGACACCGGTCGATCCATCTAACACCCCAACCACCCAACTGGCGACCGGCGGcc                                                                                                                                                                            >  2:149331/1‑250 (MQ=255)
                                                                             gCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAACAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCCAAACACCAGCGACGCCCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAAACCACCCAACTGGCGACCGCAGGCCaacaac                                                                                                                                                                      >  2:540247/1‑251 (MQ=255)
                                                                                     cccgatctgCTCCAGAATTTGAGGTGAAAACAATTAAAATGAATAAAATTCTCAGCGTCCTCTGTGTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATGAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGCCCACCCAACTGGCGACCGGCGGCCAACATCTGAACGt                                                                                                                                                               <  2:333473/242‑1 (MQ=255)
                                                                                       cGCTGAGCTCCAGATTTTGCGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTcc                                                                                                                                                            >  2:285863/1‑251 (MQ=255)
                                                                                             gCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGCACGt                                                                                                                                                               >  1:333473/1‑242 (MQ=255)
                                                                                               tCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAATCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCAGTTGCGCTGCAACAACTCACCTGCACACCGGTCGATAAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACATAACATTCCCGGCAtac                                                                                                                                                    >  2:21367/1‑249 (MQ=255)
                                                                                                     ttttGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTc                                                                                                                                              >  1:519603/1‑251 (MQ=255)
                                                                                                          agctGAAAACAATGAAAATGAATGAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACTCAGACCACCCAACCGGCCACCGCCGGCAAACAACTGAACGTTCCCGGCATCCACGGTCCGGtt                                                                                                                                         >  2:149268/4‑251 (MQ=255)
                                                                                                               aaaaCAATGAAAATGAAAAAAAGTCCCATCGTCCTCTGGTTATCAGCAGGGTTAATGGCAAGCGCGGCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCCCCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTtgctgc                                                                                                                                    <  2:535982/251‑1 (MQ=255)
                                                                                                                   cAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGACGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACAACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTAc                                                                                                                                >  1:265198/1‑251 (MQ=255)
                                                                                                                     atgaaactgaaTAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGCCACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACACCTGAACGTTCCCGCCATAAGTGGTCCGGTTGCTGCGAACa                                                                                                                               >  2:536091/1‑250 (MQ=255)
                                                                                                                          aatgaaTAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCCACCACTGAACGTTCCCGGCACCAGTGGTCCGATTGCTGAGTACACCGTcc                                                                                                                         >  2:197493/1‑251 (MQ=255)
                                                                                                                                    aGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGCTGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTACGGTTGCTGCGTACAGCGTCCCGCCAAACAt                                                                                                               >  2:467578/1‑251 (MQ=255)
                                                                                                                                                 tctGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTgccc                                                                                                  >  1:360113/1‑248 (MQ=255)
                                                                                                                                                  ctGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGAccc                                                                                                 >  1:35764/1‑251 (MQ=255)
                                                                                                                                                          cagcagGGTTACTGGCAAGCGCGCCTGGACTTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCACACATTGGCGAACTGACCCTGACGCtt                                                                                         >  1:399422/1‑250 (MQ=255)
                                                                                                                                                           agcaaGGTTCCTTGCAACCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCTGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGACCGGGTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTg                                                                                         <  2:368860/250‑1 (MQ=255)
                                                                                                                                                             cagGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGAcc                                                                                      <  1:306619/251‑1 (MQ=255)
                                                                                                                                                               gggTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGAcc                                                                                      >  1:35969/1‑249 (MQ=255)
                                                                                                                                                                  ttACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTTGTCCGGTTGCTGCGAACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCg                                                                                  >  2:427580/1‑250 (MQ=255)
                                                                                                                                                                      tGGCAAGCGCGACTCGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACAGAACGTTACCGGCATCAGTGGTCCGGCTGCTGCGTACAGCGACACGGCAACCATTGACGAACTGACCCTGACGCTGACCAGAGAAGTg                                                                             >  2:6518/1‑251 (MQ=255)
                                                                                                                                                                           aGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGaaca                                                                         >  1:110580/1‑250 (MQ=255)
                                                                                                                                                                                   tGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAACCAAAcc                                                                 >  1:220425/1‑250 (MQ=255)
                                                                                                                                                                                    ggAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCa                                                                <  1:304958/250‑1 (MQ=255)
                                                                                                                                                                                      aaTTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGCTCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAAACAGCg                                                             >  2:119632/1‑251 (MQ=255)
                                                                                                                                                                                       atttgCCTTGCCGTTGTTAACTCCGTACCGCAAAACTCCAGCGCCGCGCCTGCCATTCCATCTGTTTCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCCCCCAACTGGCGACCGGCGTCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATGGGCGAACTGCCCCTGACGCTGACCAGCGAAGTGAACAAACAACCCAGCGt                                                            <  2:229578/247‑1 (MQ=255)
                                                                                                                                                                                           cccTTGACGATGTTAACTACGTACCGCAAAACTCCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCGGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGttttt                                                        <  2:519603/250‑1 (MQ=255)
                                                                                                                                                                                              ttGCCGATGTTAACTACGTACCGCAAAACACCAGCGCCGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGACATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAACCATTGGCGAACTGACCCTGACGCTGCCCAGCGAAGGGAACAAACAAACCAGCCTTTTTgcg                                                     >  2:447060/1‑251 (MQ=255)
                                                                                                                                                                                                    acgtctgtgtcGTACCGAAAACCCCCAGCGACCCACCCGTCATTCCATCTGTTGCCCTGCAACAACGCGCCTAGACACCGGTCGACCAATCCAAAACCCAGACCACCCAACGGGCGAACGGCGGCCAACATCTGAACGTTCCCGGCATCATTGGTCCGGTTGCGGCGTACAGCGACCTGGCAAACATTGGCGAACTGACCCTGACGCTGTCCAGCTAAGTGTACAACCAAAGCAGCGTTTTTGCGCCCAAc                                               <  2:239709/241‑1 (MQ=255)
                                                                                                                                                                                                        tAACTACGTACCGCACAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACACAGACCACCCACCTGGCGACCGGCGGCCAACAAATGAACGTTCACGGCATCAGTggac                                                                                                                                              >  2:515356/1‑150 (MQ=255)
                                                                                                                                                                                                        tAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTTCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTc                                                                                                                                              <  1:515356/152‑1 (MQ=255)
                                                                                                                                                                                                             aCGTCCCGCAAAACACCAGCGACGCGCCAGCCATGCCATCTGTGGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCTCCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCAGCATTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACAGTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAACGTGCTGa                                        <  2:220425/249‑1 (MQ=255)
                                                                                                                                                                                                             aCGTACCGCAAAACCCCAGCGGCGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAACGTGCTGAtt                                      <  1:536091/251‑1 (MQ=255)
                                                                                                                                                                                                              cGTACCGCAAAACACCAGCGACACGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAACGTGCTGATTc                                     >  2:113347/1‑251 (MQ=255)
                                                                                                                                                                                                                tACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAGGTGAACAACCAAACCAGCGATTTTACGCGTCACGTGCTGATTCtt                                   >  2:102595/1‑251 (MQ=255)
                                                                                                                                                                                                                tACCGAAAAACGCCAGCGTCGCGACATCCGTTTCATCTGATGTGCTGCAACACCTCACCTGGAAACCGGACGATCTATCTAAAACCCAGCCCACCCAACTGGCGACCGGAGGGCAACAGCTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGTACTGTCCCTGACGCTGACCAGCGAAGTGGACAAACAAACCAGCTTTTTTTCGTCGAACGTGCTGATTCt                                    <  2:168963/250‑1 (MQ=255)
                                                                                                                                                                                                                      aaaaCACCATCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAACGTGCTGATTCTTGATCAg                             <  1:540247/251‑1 (MQ=255)
                                                                                                                                                                                                                      aaaaCACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAACGTGATGCTTCTTGATc                               >  1:89889/1‑249 (MQ=255)
                                                                                                                                                                                                                          cacCAGCGACACGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAACGTGCTGATTCTTGATCAGAACa                         <  1:113347/251‑1 (MQ=255)
                                                                                                                                                                                                                            ccAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAACGGGCTGATTCATGATCAGAACATg                       >  2:102523/1‑251 (MQ=255)
                                                                                                                                                                                                                             cAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAACGTGCTGCTTCTTGATCAGAACATg                       >  2:99695/1‑250 (MQ=255)
                                                                                                                                                                                                                                cGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAACGTGCTGATTCTTGATCAGAACATGAccc                   <  1:56044/251‑1 (MQ=255)
                                                                                                                                                                                                                                    gcgcCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAACGTGCTGATTCTTGATCAGAACATGACCCCATc               <  1:285863/251‑1 (MQ=255)
                                                                                                                                                                                                                                     cgcCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGGCCGGTTGCTGCGTACAGAGTCCCGGCAAACATTGGCGCACTGACCCTGACGCTGACCAGCGTAGGGAACAAACAAACCAGCGTTTTTGCGCCGAACGTGCTGATTCTTGATCAGAACATGACCCCATc               >  2:180791/1‑250 (MQ=255)
                                                                                                                                                                                                                                        cagccaTTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGGTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAAATGACCCTGACGCTGCCCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAACGTGCTGATACTTGATCAGAACATGACCCCATCAGcc           >  2:517825/1‑251 (MQ=255)
                                                                                                                                                                                                                                            aaTTCCATCTGTTGTGCTGCAACAACTCCCCTGGACTCCGGTCGTACAATCTAAAACCCATACCTCCCAACTGGCGACCGGCGGCCAACAACTGAACGTGCCCGGCATCAGTGGTCCGGTTGCTTCGTACAGCGGCCCGGCAAATATAGGCGAACTGACCCTTACGCTGTCCAGCGAAGTGAACAAACAACCCAGCGTCTTTTCGCCGAACGTGCTTATTCTTGTTCAGAACCTTTCCCCATCAGCcttct       <  2:343483/250‑1 (MQ=255)
                                                                                                                                                                                                                                                  aTCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGAACACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACAAGCGTTTTTGCGCCGAACGTGCTGATTCTTGATCAGAACATGACACCATCAGCCTTCTTcccc  >  1:60705/1‑250 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
GCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGCTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAACGTGCTGATTCTTGATCAGAACATGACCCCATCAGCCTTCTTCCCC  >  NC_000913/4249437‑4249928

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