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.4% 201.6 / NA 62A44V (GCT→GTT) malMmaltose regulon periplasmic protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (34/27);  minor base G (0/1);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGCTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAACGTGCTGATTCTTGATCAGAACATGACCCCATCAGCCTTCTTCCCC  >  NC_000913/4249438‑4249928
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
cATCCGGCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTgttgcg                                                                                                                                                                                                                                                  <  1:438142/251‑1 (MQ=255)
     ggCATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGca                                                                                                                                                                                                                                             <  2:447735/251‑1 (MQ=255)
       cATAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGcaac                                                                                                                                                                                                                                           >  2:475279/1‑251 (MQ=255)
        aTAAAAACAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGcaac                                                                                                                                                                                                                                           >  1:503009/1‑250 (MQ=255)
             aaCAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAAGAAAAGACTCATCGTCCTCTGTTTAGCAGCAGGGTTACTGGCAAGCGCGCCTTGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCa                                                                                                                                                                                                                                     <  1:360252/251‑1 (MQ=255)
              aCAGGTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTGGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATGAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCAc                                                                                                                                                                                                                                    <  2:60861/251‑1 (MQ=255)
               cAGGTTGTCATTATCTGAAAGGGGCTAAAGCCCCTCTGATTAGCGGTTTTAGCGCGCTATAGCCTGGCTACCGCTGAGCTCCAGATTTTTAGGTGAAAACAATGAAAATGAATAAAATGCTCATCGTCCTCTGTTTATCAGCATGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCTCAAAACTCCAGCGACGCGCCAGCCATTCCATCTGTGGCTCTGcaacatttctcc                                                                                                                                                                                                                                   <  2:356398/251‑8 (MQ=255)
                  gTTGTCATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTgg                                                                                                                                                                                                                                >  1:159544/1‑251 (MQ=255)
                       cATTATCTGAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGTGTTACTGTCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCTCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGacac                                                                                                                                                                                                                            <  1:92832/250‑1 (MQ=255)
                               gAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGTTGCAACAACTCACCTGGACACCGGTCGATc                                                                                                                                                                                                                   >  2:98379/1‑251 (MQ=255)
                               gAAAGGGGCGAAAGCCCCTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCg                                                                                                                                                                                                                      <  1:466567/248‑1 (MQ=255)
                                              cccTCTGATTATCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAAccc                                                                                                                                                                                                       >  1:184290/1‑248 (MQ=255)
                                                        aTCGGGTTTAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCGTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCGAAAACCCAGACCACCCAACt                                                                                                                                                                                          <  2:345908/251‑1 (MQ=255)
                                                                tAGCGCGCTATTGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACcg                                                                                                                                                                                  >  1:440787/1‑251 (MQ=255)
                                                                          ttGCCTGGCTAACGCTGAGCTCAAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAACCCCAGACCACCCAACTGGCGACCGGCGGCcaac                                                                                                                                                                         >  2:54967/1‑250 (MQ=255)
                                                                           tGCCTGGCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCcaaca                                                                                                                                                                        >  1:310163/1‑250 (MQ=255)
                                                                               gggCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGa                                                                                                                                                                   <  2:321778/250‑1 (MQ=255)
                                                                                ggCTACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTg                                                                                                                                                                    <  1:475279/249‑1 (MQ=255)
                                                                                   tACCGCTGAGCTCCAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGt                                                                                                                                                               >  1:281657/1‑251 (MQ=255)
                                                                                                cAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCCGACCACCCAACTg                                                                                                                                                                                         >  2:488546/1‑212 (MQ=255)
                                                                                                cAGATTTTGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGGCGATCAATCTAAAACCCAGACCACCCAACTg                                                                                                                                                                                         <  1:488546/212‑1 (MQ=255)
                                                                                                      ttGAGGTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTcc                                                                                                                                             <  1:424533/250‑1 (MQ=255)
                                                                                                           gTGAAAACAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGGCCGGTtgc                                                                                                                                       >  1:108891/1‑251 (MQ=255)
                                                                                                              aaaaCAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTtgctgc                                                                                                                                    >  1:329530/1‑251 (MQ=255)
                                                                                                                 aCAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCCGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCATTATAAAACCCAGACCACCCAACTGGCACCCGGGGCCCAACAACTGAACGTTCCCGGCATAAGTGGTCCGGTTGCGGCGtc                                                                                                                                 >  1:207021/1‑250 (MQ=255)
                                                                                                                 aCAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCgg                                                                                                                                                                                                                         >  1:228692/1‑163 (MQ=255)
                                                                                                                 aCAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCgg                                                                                                                                                                                                                         <  2:228692/163‑1 (MQ=255)
                                                                                                                 aCAATGAAAATGAATAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGt                                                                                                                                  >  2:216314/1‑250 (MQ=255)
                                                                                                                         aatgaaTAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTcc                                                                                                                         <  1:216314/251‑1 (MQ=255)
                                                                                                                         aatgaaTAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTcc                                                                                                                         >  1:156479/1‑251 (MQ=255)
                                                                                                                           tgaaTAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCCTCAGTGGTCCGGTTGCTGCGTACAGCGTCCCg                                                                                                                       >  1:407278/1‑251 (MQ=255)
                                                                                                                              aTAAAAGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGc                                                                                                                     >  2:4779/1‑250 (MQ=255)
                                                                                                                                aaaaGTCTCATCGTCCTCTGTTTATCAGCAGGGTTTCTTGCCAGCGCGCCTGGAATTAGCCTTGCCAATGTTAACTCCGTACCGCAAAAAACCAGCGTCGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCTCCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTTGTCCGGTGGGTGCGTACAGCTTCCCGGCaaa                                                                                                                  <  2:394157/251‑1 (MQ=255)
                                                                                                                                aaaaGTCTCAGCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCaa                                                                                                                   <  2:73025/250‑1 (MQ=255)
                                                                                                                                 aaaGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTAGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAAc                                                                                                                 >  1:142388/1‑251 (MQ=255)
                                                                                                                                  aaGTCTCATCGTCCTCTGGTTAGCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGGGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACa                                                                                                                <  2:281657/251‑1 (MQ=255)
                                                                                                                                   aGTCTCATCGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGACTGGAATTAGCCTTGCCGATGTTAAATACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCACTGCAACAAATAACCTGGACACCGGTCGATCAATCTAACACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCCGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACAt                                                                                                               >  1:585892/1‑251 (MQ=255)
                                                                                                                                       tcATCGTCCTCTGATTCTCAGCAGGGTTACTGGCAAGCGCGCATGGAATTAGCATTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCACTCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGACAACAACTGAACGTTCCCGGCATCAGTGGTCCGATTGCTGCGTACAGCGTCCCGGAAAACATTgcc                                                                                                           >  1:185415/1‑249 (MQ=255)
                                                                                                                                        caacGTCCTCTGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCTCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCGAAAGCCCAGACCACCCAACTGGCTACCGGCGGCCAACAACTGAGCGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTg                                                                                                             <  2:279619/245‑1 (MQ=255)
                                                                                                                                               ctctGTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGAc                                                                                                   >  2:126007/1‑251 (MQ=255)
                                                                                                                                                   gTTTATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGAcc                                                                                                  <  2:421513/248‑1 (MQ=255)
                                                                                                                                                    tttATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTg                                                                                               >  1:519749/1‑250 (MQ=255)
                                                                                                                                                    tttATCAGCAGGGTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGAGTCCCGGCAAACATTGGCGAACTGACCCTgc                                                                                              >  1:318541/1‑250 (MQ=255)
                                                                                                                                                               ggTTACTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGc                                                                                   <  1:4779/251‑1 (MQ=255)
                                                                                                                                                                   aCTGGTAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGGACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGGTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGCCCACCCAACTGGCGGCCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGaa                                                                                <  2:322265/250‑1 (MQ=255)
                                                                                                                                                                   aCTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGaa                                                                                <  1:98379/250‑1 (MQ=255)
                                                                                                                                                                    cTGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGAACAGCGAAgg                                                                              >  2:250787/1‑250 (MQ=255)
                                                                                                                                                                     tGGCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTACCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTg                                                                             >  1:405983/1‑251 (MQ=255)
                                                                                                                                                                      ggCAAGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGa                                                                            <  1:20829/251‑1 (MQ=255)
                                                                                                                                                                         aaGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGaac                                                                          >  1:466994/1‑250 (MQ=255)
                                                                                                                                                                         aaGCGCGCCTGGAATTAGCCTTGCCGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGaac                                                                          >  1:462664/1‑250 (MQ=255)
                                                                                                                                                                                              tGCCGATGTTAACTACGTACCGCAAAACACCAGTGACGCGCCAGCCATTCCAACTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTgc                                                      <  2:421487/249‑1 (MQ=255)
                                                                                                                                                                                                ccGATGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAATCCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTTACCCTGACGCTGACCAGCGAAGGGAACAAACAAACCAGCGTTTTTGCGCCg                                                  <  2:14684/251‑1 (MQ=255)
                                                                                                                                                                                                   aTGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTTGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAAc                                               <  1:238390/251‑1 (MQ=255)
                                                                                                                                                                                                   aTGTTAACTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAAc                                               >  2:167139/1‑251 (MQ=255)
                                                                                                                                                                                                          cTACGTACCGCAAAACACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCAACTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGAGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCCGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAACGTGCTgt                                        >  2:238131/1‑250 (MQ=255)
                                                                                                                                                                                                                       aaCACCAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAACGTGCTGATTCTTGATCAGa                            >  2:162900/1‑250 (MQ=255)
                                                                                                                                                                                                                           ccAGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCCCCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAACGTGCTGATTCTTGATCAGAACATg                       <  2:306175/251‑1 (MQ=255)
                                                                                                                                                                                                                             aGCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAACGTGCTGATTCTTGATCAGAACATGAc                     >  1:558924/1‑251 (MQ=255)
                                                                                                                                                                                                                              gCGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGTTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAACGTGCTGATTCTTGATCAGAACATgtcc                    <  1:517783/251‑4 (MQ=255)
                                                                                                                                                                                                                               cGACGCGCCAGCCATTCCATCTGTTGCGCTGCAACAACTCACCTGGACCCCGGTCGATCAATCTAACACCCAGCCCACACAACTTGCGACCGGCAGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGCTGCTGCGTACAGCGTCCCGGCAACCATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGACCAAACAAACCAGCGTTTTTGCGCCGAACGTGCTGATTCTTGATCAGAACATGAccc                   >  1:263850/1‑251 (MQ=255)
                                                                                                                                                                                                                                               ccATCTGTTGCGCTGCAACAACTCACCTGGACACCGGTCGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACAAACCAGCGTTTTTGCGCCGAACGTGCTGATTCTTGATCAGAACATGACCCCATCAGCCTTCTTcc    <  1:437840/250‑1 (MQ=255)
                                                                                                                                                                                                                                                  tCTGTTGCGCTGCAACAACTCACCTGGACACCGGACGATCAATCTAAAACCCAGACCACCCAACTGGCGACCGGCGGCCAACAACTGAACGTTCCCGGCATCAGTGGTCCGGTTGCTGCGTACAGCGTCCCGGCAAACATTGGCGAACTGACCCTGACGCTGACCAGCGAAGTGAACAAACACACCAGCGATTTTGCGCCGAACGTGCTGATTCTTGATCAGAACATGACCCCATCAGCCTTCTTcccc  >  2:185564/1‑249 (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: