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,6840CT96.6% 191.5 / NA 59A44V (GCT→GTT) malMmaltose regulon periplasmic protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (31/26);  minor base G (0/2);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 2.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGCTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAACGTGCTGATTCTTGATCAGAACATGACCCCATCAGCCT  >  NC_000913/4249437‑4249920
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                            
gcATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCtgttgc                                                                                                                                                                                                                                           <  1:182173/251‑1 (MQ=255)
           aaaatCAGGTTGTTATTATCTGAAAGGGGCGAAAGCCCCTCTGAATATCGGGTTTAGCGCGCTGTTTCCCGGCTCCCGCTGAGCGTCAGATGTTGAGGTTAAAACAATGAAAAGGAATTAAAGTCGTAGCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGTATTAGCCGTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCAGGCTATCGGGTGCGCTGcaacaa                                                                                                                                                                                                                                 <  2:219749/250‑1 (MQ=255)
             aaaCAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTc                                                                                                                                                                                                                              >  1:58443/1‑251 (MQ=255)
                      gTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGac                                                                                                                                                                                                                      <  2:390877/250‑1 (MQ=255)
                      gTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCGGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGGTAACTACGTACCGCAAATCTCCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGac                                                                                                                                                                                                                      <  1:20916/250‑1 (MQ=255)
                        cATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACAcc                                                                                                                                                                                                                   <  1:460276/251‑1 (MQ=255)
                           aaTCTGAAAGGGGCGAAAGCCCCTCTGATTTTCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTCAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGt                                                                                                                                                                                                                <  2:58671/250‑1 (MQ=255)
                            atccggaaGGGGCGAAAGCCCCCCTGAATATCGGGGTTATGGCCGCATTTCCTGGGTACCGCTGAGCTCCAGATTTTGAAGTGAAAACCAATGAAATGAATAAAAGTCGCCTCCTCCTCCGTTTATCAGCCGGGTTACTGGCAAGCGCGCCTTGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTc                                                                                                                                                                                                               <  1:408289/245‑1 (MQ=255)
                                   aGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAAt                                                                                                                                                                                                        <  1:38397/251‑1 (MQ=255)
                                      ggCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGGCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCTCCTGGAATTAGCCTGGCCGATGTTAACTACGTACCGCAAAACACCAGCGGCGCGCCAGCCATTCCATCTGTTGCGCTGCAACATCTCTCCTGGACGCTGGTCGATCGATCTa                                                                                                                                                                                                     <  2:213987/251‑1 (MQ=255)
                                         gAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTaaa                                                                                                                                                                                                   >  1:133130/1‑250 (MQ=255)
                                             gCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCCGATTTTGCGGTGAACACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTCAAAccc                                                                                                                                                                                               >  2:450555/1‑250 (MQ=255)
                                                      aTTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCa                                                                                                                                                                                     <  2:147380/251‑1 (MQ=255)
                                                                         cTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTACCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGccc                                                                                                                                                                   >  2:399494/1‑250 (MQ=255)
                                                                              gCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCaacaac                                                                                                                                                              >  1:364803/1‑251 (MQ=255)
                                                                                cTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGAAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTg                                                                                                                                                            >  1:291485/1‑251 (MQ=255)
                                                                                  ggCTACCGCTGAGCTACAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTGATCAGCAGGGTTACTTGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCGCCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGa                                                                                                                                                           <  2:233864/250‑1 (MQ=255)
                                                                                                 ccAGATTTTGAGGGGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGTCCACCCAACTGGCGACCTGCGGCCAACAACTTAACGTTCCCGGCATCAg                                                                                                                                           <  2:419506/251‑1 (MQ=255)
                                                                                                  ccgatttagcggtcacaaCAATGAAAATGAATAAACGTCTCATTGTCCTATGTATATCAGCAGGGTTCCTGGCAAGAGCACCTGGAATTAGCCTTGCCGCTGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTACATCTGTTGCGCTGCAACAAAGCACCTGGACACCGGTCGATCACTCTAAAACCCAGACCACCCACCTGACGACCGGCGCCCAACAAATGACCGTTCCCGGCATCAGt                                                                                                                                          >  2:29922/17‑251 (MQ=255)
                                                                                                       tttGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGAGACCGGCGGCCAACAACTGCACGTTCCCGGCATCAGCGGTcc                                                                                                                                     >  1:231212/1‑251 (MQ=255)
                                                                                                         tGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCgg                                                                                                                                   >  1:101834/1‑251 (MQ=255)
                                                                                                               gAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTtgctg                                                                                                                             >  1:291027/1‑251 (MQ=255)
                                                                                                                       atgaaactgcaTAAAAGTCTAATCGTCCTCTGTTTATCAACAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGCTCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAAAGTTCCCGGCATCCGTGGACCGGTTGCTGCGGACa                                                                                                                       >  2:270823/1‑250 (MQ=255)
                                                                                                                          aaaatgaaTACAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGt                                                                                                                   >  2:196568/1‑251 (MQ=255)
                                                                                                                           aaatgaaTAAAAGTCTCAGCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATGAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCGGCGTACAGCGTc                                                                                                                  <  2:431283/251‑1 (MQ=255)
                                                                                                                             atgaaTAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTccc                                                                                                                >  2:391454/1‑251 (MQ=255)
                                                                                                                                   aaaaGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGATGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCaa                                                                                                           >  2:271605/1‑250 (MQ=255)
                                                                                                                                    aaaGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAAc                                                                                                         <  1:391454/251‑1 (MQ=255)
                                                                                                                                          tcagcatCCTCTGTTAATCAGCAAGGTTACTGGCAAGCCCACCAGGCATTAGCCTTGCCGATGCTCCATACGTCCCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGCCGATCAATCTAAAACCCAGACCACCCCACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCCTCAGGGGTACGGTTGCTGCGTACAGCggacaagcaaacatggcc                                                                                                   >  2:195543/7‑234 (MQ=255)
                                                                                                                                           cATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCg                                                                                                  >  2:129744/1‑251 (MQ=255)
                                                                                                                                           cATCGACCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCATGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGCTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCg                                                                                                  >  1:357881/1‑251 (MQ=255)
                                                                                                                                                tCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTg                                                                                             >  2:382063/1‑251 (MQ=255)
                                                                                                                                                  cgctGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGGCAACAACTGAACGGTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCATACATTGGCGAACTGAc                                                                                           <  2:211371/249‑1 (MQ=255)
                                                                                                                                                          aTCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTTGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACg                                                                                    <  1:75265/250‑1 (MQ=255)
                                                                                                                                                               cagGGTTACTGGCAAGCGCGCATTGAATTAGCCTGGCCGATGTTAACTACGTACCGCAAAACACCACCGACGCGCCAGCCATTCCAGCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTAGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGAcc                                                                              <  2:112478/251‑1 (MQ=255)
                                                                                                                                                                  ggTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGc                                                                           <  2:400968/251‑1 (MQ=255)
                                                                                                                                                                                 cgcCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCAGTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAacaaacaa                                                            <  1:55477/251‑1 (MQ=255)
                                                                                                                                                                                  gcCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGCACGTTCCCGGCATCAg                                                                                                                                           >  1:85348/1‑171 (MQ=255)
                                                                                                                                                                                  gcCTGGAATTAGCCTTGCCGATGTTAACTACGGACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGTCATCAg                                                                                                                                           <  2:85348/171‑1 (MQ=255)
                                                                                                                                                                                       gAATTAGCCTTGCCGATGTGAACTACGTACCGCAAAACACCAGAGACTGGACAGCCATTCCATCTGTTGCGCTGCAACAACTCTGCTGGACACCGGTCGATCAAGCTAAAACCCAGACCTCCCAACGGGCGACCGGCGGCCAACAACGGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGGCCCGGCAAACATTGGCGAACGGACCCGGACGCTGACCAGCGAAGTGATCAAACATACCAGc                                                      <  2:40447/251‑1 (MQ=255)
                                                                                                                                                                                        acttAGCCTTCCCGTTGTCACCTACGCACCACACACCACCTGCGCCGCGCCAGCCATTCCATCTGTTGCGCTGCCACAAATCACCTGGCCACCGGTCGATCAATCTACAACCCAGACCACCCCACTCGCGACCGGCGCCCAACAACTGACCGATCCCGGCATCAGTGGTCCGGTTGCCCCGTACCGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGATGACCAGCGACACGACCAAACAAACCAGCg                                                     >  2:54803/3‑251 (MQ=255)
                                                                                                                                                                                            aGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGtttt                                                 <  1:196568/251‑1 (MQ=255)
                                                                                                                                                                                                  gCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGCGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGcc                                           >  2:364970/1‑251 (MQ=255)
                                                                                                                                                                                                      aTGTTAACTACGTACCGCAAAACACCAGCGACGAGCCCGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGACTTTGCGCCGAAc                                       >  2:300651/1‑251 (MQ=255)
                                                                                                                                                                                                        gttatctacggtccGCAAAACACCAGCGACGTGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAAACTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTGTTGCGCCGAACGt                                     <  2:189894/239‑1 (MQ=255)
                                                                                                                                                                                                        gTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAACgg                                     >  1:164338/1‑250 (MQ=255)
                                                                                                                                                                                                         atAACTACGTGCCGCAAAACACCAGCGACGCGCCAGCCATTCCAGCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAACCCCAGACCACCCAACGGGCGACCGGCTGCCTACAACTGAACGTTCCCGGCAGCAGTGGTCCGGTGGCTGCGTACGGCTTCCCGGCAAACATGGGCGAACTGACCCTGACGCCGATCATCTAAGTGAACAAACAAACCAGCGTTTTTGAGCCGAACGTg                                    <  2:166806/250‑1 (MQ=255)
                                                                                                                                                                                                          tAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAACGTGc                                   >  1:402249/1‑251 (MQ=255)
                                                                                                                                                                                                            aCTACGTACCGCACAACACCAGCGACACGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCACTCTAAAACCCAGACCACCCAACTGACGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCACGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAACGTGCTg                                 >  2:320809/1‑251 (MQ=255)
                                                                                                                                                                                                            aCTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAACGTGCTg                                 <  2:133130/251‑1 (MQ=255)
                                                                                                                                                                                                                 gTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAACGTGCTGATTCt                            <  1:303708/251‑1 (MQ=255)
                                                                                                                                                                                                                          aaCACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAACGTGCTGATTCTTGATCAGa                    >  2:49606/1‑250 (MQ=255)
                                                                                                                                                                                                                             acCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTccc                                                                                                                >  1:203796/1‑155 (MQ=255)
                                                                                                                                                                                                                             acCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTccc                                                                                                                <  2:203796/155‑1 (MQ=255)
                                                                                                                                                                                                                                aGCGACGCGCCAGCCATTCCATCTGTTGCGATGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAACGTGCTGATTCTTGATCAGAACATGAc             <  1:271605/251‑1 (MQ=255)
                                                                                                                                                                                                                                   gACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACACAGACCACCCAACTCGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCCGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAACGTGCTGATTCTTGATCAGAACATGAcccc          >  1:454375/1‑251 (MQ=255)
                                                                                                                                                                                                                                   gACGCGCCAGCCATTCCATCTGGTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCCCCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAACGTGCTGATTCTTGATCAGAACATGAcccc          <  2:140419/251‑1 (MQ=255)
                                                                                                                                                                                                                                     cgcgCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAACGTTTTTGCGCCGAACGTGCTGATTCTTGATCAGAACATGACCCCAt        >  2:435392/1‑251 (MQ=255)
                                                                                                                                                                                                                                          cagccaTTCCATCTGTTGCGATGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGACGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAACGTGCTGATTATTGATCAGAACATGACCCCATCAGcc   >  1:160283/1‑251 (MQ=255)
                                                                                                                                                                                                                                           agccaTTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCACTCAATCTAAAAACCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGAACAGCGTACCGGCAAAAATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAAACAGCGTTTTTGCGCCGAACGTGCAGATTCTTGATCAGAACATGACCCCATCAGCct  >  1:305979/1‑251 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                            
GCATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGCTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAACGTGCTGATTCTTGATCAGAACATGACCCCATCAGCCT  >  NC_000913/4249437‑4249920

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: