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.6% 244.4 / NA 71A44V (GCT→GTT) malMmaltose regulon periplasmic protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (47/23);  minor base A (1/0);  total (48/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGCTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAACGTGCTGATTCTTGATCAGAACATGACCCCATCAGCCTTCTTCCCC  >  NC_000913/4249438‑4249928
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
cATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTgttgcg                                                                                                                                                                                                                                                  >  2:9352/1‑251 (MQ=255)
   ccGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTgttgcgctg                                                                                                                                                                                                                                               <  1:271339/251‑1 (MQ=255)
   ccGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTgttgcgctg                                                                                                                                                                                                                                               <  1:181628/251‑1 (MQ=255)
            aaaCAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTc                                                                                                                                                                                                                                      >  2:270284/1‑251 (MQ=255)
             aaCAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTc                                                                                                                                                                                                                                      >  2:41016/1‑250 (MQ=255)
                    tGTCATTATCTGCAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCACCAACTCACCTGGac                                                                                                                                                                                                                              >  1:484600/1‑251 (MQ=255)
                       cATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGCATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACAcc                                                                                                                                                                                                                           >  1:371523/1‑251 (MQ=255)
                       cATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGacac                                                                                                                                                                                                                            <  2:479129/250‑1 (MQ=255)
                           aTCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTc                                                                                                                                                                                                                       >  1:351974/1‑251 (MQ=255)
                            tCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCACCAACTCACCTGGACACCGGTCg                                                                                                                                                                                                                      >  2:462298/1‑251 (MQ=255)
                               gAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCAGTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTc                                                                                                                                                                                                                       <  2:60448/247‑1 (MQ=255)
                                  aGGGGCGACAGCCCCTCTGATTATCGGGTTTAGCGCGCTCTTGCCTGGCTACCGCTGATCTCAAGATTTTGAGGTGACAACAATGAAAATGACTAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAGGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCAGATGCGCTGCAACAACTCACCTGGACACCGATCGAACAAt                                                                                                                                                                                                                >  2:544394/1‑251 (MQ=255)
                                       cGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTaa                                                                                                                                                                                                            >  2:359798/1‑250 (MQ=255)
                                                    gATTATAGGGTTGAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTGGATGTGAAAACAATGAAAATGAAAAAAAGTCTCATCGTCCTCTTTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATGAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCACCTACTCGCCTGGACACCTGTCGATCAATCTAAACCCCAGaccacc                                                                                                                                                                                               <  2:300936/250‑1 (MQ=255)
                                                                      gcTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGcc                                                                                                                                                                            >  2:117550/1‑251 (MQ=255)
                                                                            gCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGGTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCcaacaa                                                                                                                                                                       <  2:263786/250‑1 (MQ=255)
                                                                                ggCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAAGAAAAGGCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCCATCGCGCCTGGAATTATCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCTGACAACAACTGaa                                                                                                                                                                  <  2:25927/251‑1 (MQ=255)
                                                                                     ccGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTc                                                                                                                                                             >  2:496357/1‑251 (MQ=255)
                                                                                      cGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAATAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACGCCCCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACGGAACGTTcc                                                                                                                                                            <  1:419066/251‑1 (MQ=255)
                                                                                       gCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTccc                                                                                                                                                           <  2:495136/251‑1 (MQ=255)
                                                                                       gCTGAGATCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTccc                                                                                                                                                           >  1:469674/1‑251 (MQ=255)
                                                                                                cAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGt                                                                                                                                                  >  1:482768/1‑251 (MQ=255)
                                                                                                           gTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTtgc                                                                                                                                       >  1:177013/1‑251 (MQ=255)
                                                                                                            tGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTtgct                                                                                                                                      >  1:245396/1‑251 (MQ=255)
                                                                                                             gAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTtgctg                                                                                                                                     >  1:446736/1‑251 (MQ=255)
                                                                                                                aaCAATGAAAATGAATAAAAGTCTCATCGTCGTCTGTTTATCAGCAGCGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGCTGTTAACTACGTAACGCAAAACACCAGCGAAGCGACAGCCATTCCATCTGTTGCGATGCAACAACTCACCTGGACACCGGTCCATCAATCTAAAACCCAGACCACCCAACAGGAGTCCGGCGACCAACAAATGAACGTTCCCGGCATCAGAGGTCCGGTTGCTGCGt                                                                                                                                  >  1:244412/1‑251 (MQ=255)
                                                                                                                       aaaatgaaTAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGt                                                                                                                           >  1:190561/1‑251 (MQ=255)
                                                                                                                       aaaatgaaTAAAAGTCTCATAGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGt                                                                                                                           >  1:224264/1‑251 (MQ=255)
                                                                                                                           tgaaTAAAAGTCTCATCGACCTCTGTTTATCAGCAGGGTTACTGACAAGCCCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAAAACCCGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGCCACCGGTCGATCAATCTAAAACCCAGACCACACAACTGGCGACCGGCGGCCAACAACTGAACGTCCCCGCCATCAGTGGTCAGGTTGCTGCGTACAGCGTCCCg                                                                                                                       >  1:72313/1‑251 (MQ=255)
                                                                                                                             aaTAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGc                                                                                                                     >  1:114563/1‑251 (MQ=255)
                                                                                                                              aTAAAAGTCACATCGTCCTCTGTTTAGCAGCAGGGTTACTGGCAAGCTCGCCTTGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCAGTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTTGAGCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGGACGTTCCCGGCATGAGTGGTCCGGGTGCTGCGTACAGCGTCCCgtca                                                                                                                    <  2:484600/251‑4 (MQ=255)
                                                                                                                                  aaGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAAc                                                                                                                 >  1:266341/1‑250 (MQ=255)
                                                                                                                                   aGTCTCATCCTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTCCGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACACAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCCAACAt                                                                                                               >  2:377923/1‑251 (MQ=255)
                                                                                                                                    gTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCa                                                                                                                                                                                                                  <  1:434755/151‑1 (MQ=255)
                                                                                                                                    gTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCa                                                                                                                                                                                                                  >  2:434755/1‑151 (MQ=255)
                                                                                                                                         aTCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGcc                                                                                                                                                                            <  2:314907/184‑1 (MQ=255)
                                                                                                                                         aTCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGcc                                                                                                                                                                            >  1:314907/1‑184 (MQ=255)
                                                                                                                                          tcgtcaactgtttcaccgcagGGTTACTGGCACGCGCACCTGGAATTACACTTGCCGCTGCTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGCTCAATCTAAAACACCGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGAAATCAGTGGTCCGGTGGCTGCGTACAGCGTCCCGACAAACATTGGCGa                                                                                                         >  2:551576/18‑250 (MQ=255)
                                                                                                                                                      tATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACg                                                                                            >  1:145743/1‑251 (MQ=255)
                                                                                                                                                         cagcagGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTg                                                                                         <  1:558889/251‑1 (MQ=255)
                                                                                                                                                             agGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGAcc                                                                                      >  1:206292/1‑250 (MQ=255)
                                                                                                                                                                   aCTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAg                                                                               <  1:543401/251‑1 (MQ=255)
                                                                                                                                                                   aCTGGCAAGCGCGCCTGGAATTAGCCTTGACGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGGCCCGGCAAACATGGGCGAACTGACCCTGAAGCTGCCCAGCGaa                                                                                >  2:204665/1‑250 (MQ=255)
                                                                                                                                                                      ggCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTg                                                                             <  1:270284/250‑1 (MQ=255)
                                                                                                                                                                        cAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACAAGCGAAGTGaac                                                                          >  1:458936/1‑251 (MQ=255)
                                                                                                                                                                               gcCTGGAATTAGCCTTGACGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCAAGCCCTTCAATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATAAGTGGTCCGGTTGCTGAGTCAAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAacaaacaa                                                                    >  2:572532/1‑250 (MQ=255)
                                                                                                                                                                                ccTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACaaacaaac                                                                  <  2:190561/251‑1 (MQ=255)
                                                                                                                                                                                  tGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAAcc                                                                 >  1:437286/1‑250 (MQ=255)
                                                                                                                                                                                     aaTTAGCCTTGCCGATGTTAACTACGTACCGCAAAACCCCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTGCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCg                                                             <  2:351974/251‑1 (MQ=255)
                                                                                                                                                                                         aGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGtttt                                                         <  2:216633/251‑1 (MQ=255)
                                                                                                                                                                                            cTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTgc                                                      >  1:231621/1‑251 (MQ=255)
                                                                                                                                                                                            cTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGGTGCTGCGTACAGCGTCCCGGCACACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCACCGTTTTTgc                                                      >  2:21346/1‑251 (MQ=255)
                                                                                                                                                                                                  gATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGa                                                 >  1:292859/1‑250 (MQ=255)
                                                                                                                                                                                                     gTTAACTACGTACCGCAAAACGCCAGCGACGCTCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTTCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAACGt                                             <  1:287151/251‑1 (MQ=255)
                                                                                                                                                                                                     gTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGt                                                                              >  2:566407/1‑218 (MQ=255)
                                                                                                                                                                                                     gTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGt                                                                              <  1:566407/218‑1 (MQ=255)
                                                                                                                                                                                                       tAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCACCAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAACGTGc                                           >  1:37545/1‑251 (MQ=255)
                                                                                                                                                                                                        aaCTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGCCCCTGACGCTGACCAGCGAAGTGAACAAACACACCAGCGTTTTTGCGCCGAACGTGCt                                          >  1:364781/1‑251 (MQ=255)
                                                                                                                                                                                                        aaCTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAACGTGCt                                          >  2:433430/1‑251 (MQ=255)
                                                                                                                                                                                                          cTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAACGTGCTg                                         >  1:104504/1‑250 (MQ=255)
                                                                                                                                                                                                            aCGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAACGTGCTGAt                                       <  2:159861/250‑1 (MQ=255)
                                                                                                                                                                                                               tACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAACGTGCTGATTCtt                                   >  1:258177/1‑251 (MQ=255)
                                                                                                                                                                                                                 ccGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGCAGTGCACAAACAAACCAGCGTTTTTGCGCCGAACGTGCTGATTCTTg                                  >  2:372833/1‑250 (MQ=255)
                                                                                                                                                                                                                  cGCAAAACACAAGCGACGCGCCAGCCATTACATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACAACACAACGGGCGACCGGCGGCCAACAACTGACCGTTCCCGGCATCAGTGGTCCTGTTGCTGCGTACAGAGTCCGGGCAAACATGGGAGAACTGACCCTGACGATGACCAGCGAAGTGAACAAAAAAACCAGCGTTTTTGCACCGGACGAGCAGAGTCTTGAt                                >  1:480597/1‑251 (MQ=255)
                                                                                                                                                                                                                   gCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGGACAAACAAACCAGCGTTTTTGCGCCGAACGTGCTGATTCTTGATc                               >  2:138050/1‑251 (MQ=255)
                                                                                                                                                                                                                    cAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCACCAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAACGTGCTGATTCTTGATc                               >  2:147766/1‑250 (MQ=255)
                                                                                                                                                                                                                     aaaaCACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTTCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTGTTGCGCCGAACGTGCTGATTCTTGATCAg                             <  1:430126/251‑1 (MQ=255)
                                                                                                                                                                                                                      aaaCACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGTTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAATTGAACAAACAAACCAGCGTTTTTGCGCCGAACGTGCTGATTCTTGATCAGa                            <  1:361028/251‑1 (MQ=255)
                                                                                                                                                                                                                                 aCGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAACGTGCTGATTCTTGATCAGAACATGAcccc                  >  1:547503/1‑250 (MQ=255)
                                                                                                                                                                                                                                    cgcCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAACGTGCTGATTCTTGATCAGAACATGACCCCATCa              <  1:496357/251‑1 (MQ=255)
                                                                                                                                                                                                                                     gccagccaTTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTGTTTGCGCCGAACGTGCTGATTCTTGATCAGAACATGACCCCATCAg             <  2:99112/251‑1 (MQ=255)
                                                                                                                                                                                                                                             ttCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACAGGCGGCCAACAACTGAACGTTACCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGCACAAACAAACCAGCGTTTTTGCGCCGAACGTGCTGCTTCTTGATCAGAACATGACCCCATCAGCCttcttc     >  2:512125/1‑251 (MQ=255)
                                                                                                                                                                                                                                               ccATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCAAACGTGCTGATTCTTGATCAGAACATGACCCCATCAGCCTTCTTccc   >  2:445801/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 aTCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTACAACCCAGACCACCCAACTGGCGACCGGCGGCCCACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGCCCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAACGTGCTGATTCTTGATCAGAACATGACCCCATCAGCCTTCTTcccc  >  1:578488/1‑250 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
CATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGCTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAACGTGCTGATTCTTGATCAGAACATGACCCCATCAGCCTTCTTCCCC  >  NC_000913/4249438‑4249928

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