Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,631,146 G→A A138A (GCG→GCA)  slt → soluble lytic murein transglycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,631,1460GA92.4% 172.9 / ‑5.3 66A138A (GCG→GCA) sltsoluble lytic murein transglycosylase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (32/29);  minor base C (4/1);  total (36/30)
Fisher's exact test for biased strand distribution p-value = 3.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCGAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGCGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CAGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCA  >  NC_000913/4630904‑4631387
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
ggACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCCAAATggg                                                                                                                                                                                                                                               >  1:603647/1‑249 (MQ=255)
   cTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTCACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGCACCGGCGCGGCGTGAAGACTGGCGTTGCTTGTTTACCTTTAGCCCGGGCAAGCCCGGGACCACCGCAGAGCA‑AT‑GTACTTACTAATATTCAAAACGGAacac                                                                                                                                                                                                                                           >  2:651609/1‑251 (MQ=255)
     gaaagATTATCCACTTTATCCCTACCTGGAATCCCGTCAGATCACCGATGACCTTATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTGCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCTGGAAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAACAcc                                                                                                                                                                                                                                          <  1:471997/246‑1 (MQ=255)
      aaGGATTATCCCCCTTATCCCTGCCTGGAATTCCGTCAGATCCCCGATGATCTGATGAAGCAACCGGCGGTGACTGCCACTAACTTTGTTCGCTCTAACCCTACGCTTCCTCCCGCTCGCACGCTGCAATCTCGCTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTTGTCTTTAGCCCGGAAAAGCCAGGAACTGCCGAAGCGCA‑AT‑GTAATTACTACTATGCCAAATggtccaccgt                                                                                                                                                                                                                                        <  2:1980/251‑9 (MQ=255)
        agagtctCCGCTGTATCTCTACCTGGAATACCGACATATCACCGATGACCTGATGAATCATTCGGCGGTGACTGGCACTAACCTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCATTCTCGTTTCGTCAAGGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGTCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAACACCGGGc                                                                                                                                                                                                                                      <  2:181659/245‑1 (MQ=255)
         gATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTTACCTTTAGACCGGAAAAGACCGGGACTAACGAAGAGCACTT‑GTA‑TTAATACTATGCAAAAAGGAaccacgggc                                                                                                                                                                                                                                      >  1:410801/1‑243 (MQ=255)
              tccacTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGGGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGGAGTACAGAAGCGCA‑AT‑TTAATTACTACTATGCCAAATTGAACAACCGGCAGGGTg                                                                                                                                                                                                                                >  1:488034/5‑251 (MQ=255)
                             ccTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCA‑AT‑ATAATTACTACTATGCAAAATGGAACACCGGGCAGAGGGAAAAAACCTGGCa                                                                                                                                                                                                                   >  1:145881/1‑249 (MQ=255)
                                ggAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑gggg                                                                                                                                                                                                             <  1:673734/251‑1 (MQ=255)
                                   aTACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCg                                                                                                                                                                                                           >  1:165437/1‑250 (MQ=255)
                                             atcaacgaagaccTGATGAATCAACCGGCGGTGACTGTCACTAACTGTGTTCGCGCTAACCCCACGCTGCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCtgtg                                                                                                                                                                                                <  2:488307/239‑1 (MQ=255)
                                              tCACCGATGACCTGATGAATCAACCGGCGGTTACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCACTCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGAGTGGCTTGTTTGTATTTAGCCCGGGAAAACCCGGGAATACCGAAACGCA‑ATTGCAAT‑AATAATATGCAAAAAGGAAAACCGGGCAGAGGGAAGAAACCCGGCAAAGGGGGGAAaggctgagg                                                                                                                                                                                                >  2:645864/1‑244 (MQ=255)
                                                    atgacCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAAcc                                                                                                                                                                                         <  1:82426/251‑1 (MQ=255)
                                                      gaccTGATGAATCAACCGGCTGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCGTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGGGCTGTGGCTAACCg                                                                                                                                                                                        <  2:627631/247‑1 (MQ=255)
                                                       accTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCA‑AT‑GGAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAAACggg                                                                                                                                                                                      >  1:566074/3‑250 (MQ=255)
                                                               aaTCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAACCTTTAGCCCGGAAAAACCCGGGACTACCGAAGCGCA‑AT‑GTAATTACTACTATTCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCGGTGGCTAAACGGC‑AAGagccc                                                                                                                                                                             >  1:368735/1‑250 (MQ=255)
                                                                  cAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGagccagcc                                                                                                                                                                          <  1:532176/251‑1 (MQ=255)
                                                                      cggcggTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGGAAAACCCGGAACTACCGCAGCGCA‑AT‑GTAATTACTTATCTGCAAAATGGAACACCGGGCAGAGTGGAGAAGCCTGGGAA‑GGGGGGAAAGAGCTGTGGCTAACCGGC‑AAAAGGCAGCCTaa                                                                                                                                                                       >  1:391693/1‑250 (MQ=255)
                                                                      cggcggTGACGGTCCCTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTTAAGACTGGCGTGGCTTGTTAGCCGTTAGCCCGGAAAAGCCCGGACCTACCGAAGCGCA‑AT‑GTAATTACTACTATGCATAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAAc                                                                                                                                                                      <  1:267581/251‑1 (MQ=255)
                                                                          ggTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCGGCTCGCACGCTGCAATCTCGATTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTTGCCCGGAGAAACCCGGAAAAAACGAAGCGCA‑AT‑GTGATTACTACCATGCAACAAGGAAAAACAGGCAGAGGGAAGACAACTTGCAA‑GGGGCGAAAGAGCTAGGGCCCACCGGC‑CAGagccaggccagagcac                                                                                                                                                                  >  2:122770/1‑241 (MQ=255)
                                                                          ggTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCCAAATGGAACACCGGGCAGAGTGGAGAAGCCTGGCAA‑GGGGCGGAAGAGCTGTGGCTTACCGGC‑AAGAGCCAGCCTAACGccc                                                                                                                                                                  >  1:237614/1‑250 (MQ=255)
                                                                             gactgTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAAt                                                                                                                                                                                                                                                  <  2:532789/166‑1 (MQ=255)
                                                                             gactgTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAAt                                                                                                                                                                                                                                                  >  1:532789/5‑170 (MQ=255)
                                                                                 gtcgcTAACTTTGTTCGCGTTAACCCCCCGCTTCCTCCCTCTCTCTCGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGTCGTGAAGACTGGCGTGGCTTGATAGTCTTTAGCCCTGAAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAACCGTC‑AAGAGACCGCCTGACGTCTTTGACa                                                                                                                                                            <  2:122978/246‑1 (MQ=255)
                                                                                  tCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACa                                                                                                                                                            >  1:220780/1‑249 (MQ=255)
                                                                                    aCTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACGGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGtt                                                                                                                                                        <  2:603647/251‑1 (MQ=255)
                                                                                       aaCTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTAtt                                                                                                                                                     <  1:489045/251‑1 (MQ=255)
                                                                                       aaCTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTAt                                                                                                                                                      <  1:443889/250‑1 (MQ=255)
                                                                                            tGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAACCCCGGGCCGAGTGAAGAAGCCTGGGAA‑GGGGCGAAAGAGCTGGGGCTAACCGGC‑CAGAGCCAGCCTAACGCCTGTGACAAGTTATTTattg                                                                                                                                                >  2:499833/1‑248 (MQ=255)
                                                                                                cgcgCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTg                                                                                                                                            <  1:29728/251‑1 (MQ=255)
                                                                                                  cgcTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTggg                                                                                                                                          >  1:29197/1‑250 (MQ=255)
                                                                                                               cTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTa                                                                                                                              >  1:428555/1‑250 (MQ=255)
                                                                                                                           cGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGGGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCt                                                                                                                 >  1:129965/1‑251 (MQ=255)
                                                                                                                             cACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGCACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAATGGCACACCGGGCCGAGTTAAGAAGCCTGGCAA‑GGGGCGAAAGAACTGTGGCTAACCGGC‑CAGAGCCAGCCTAACGCCTGTGACCAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTgg                                                                                                               >  1:415633/1‑251 (MQ=255)
                                                                                                                                  tGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTTGGCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAAATACTACTATGCAAAATGGAACACCGGGCAAAGAGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGAGGCTAACCGGCCAAG‑GCCAGCCTAACGCCTTGGACAAGTTTTTTAGTGTCTGGCGTGGGGCCGGTAAAAAAAATACGCTGGCGTAt                                                                                                          >  2:184589/1‑251 (MQ=255)
                                                                                                                                    aaagccagttgcGTCAAGGAACCGGCGCGGCGTGAAGACTGGCTGGGCTGGGTCGCCTTTAGCCCGTAAAAGCCCGGAACTACCGAAGCGCA‑ATTGTA‑TTCCTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAACCGGC‑ACGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGAGTGCGTCAGGTAAACAAGATCCGCttgcttagtt                                                                                                        <  2:223674/240‑10 (MQ=255)
                                                                                                                                      aTCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACCAGATCCGCTGGCGTAttt                                                                                                        >  1:224920/1‑249 (MQ=255)
                                                                                                                                           gTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGAA‑AT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGGGTGAAAAAGGCTGGCAA‑GGGGAGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGGGACAAGTTATTTAGTGTCTGGAGTGCGTCAGGTCAACAAGATCCGCTGGCGTGTTTAAAGCGt                                                                                                 >  1:403164/1‑251 (MQ=255)
                                                                                                                                                 tCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGATAGCCTTTTGCCCGGAAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCCAAATGGCACACCGGGCAGAGTGAAGAAGCCAGGCAA‑GGGGCGAAAGAGATGTGGCTAACCGGC‑AAGAGCCCGCCTAACGCCGGTGACCAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGt                                                                                           >  1:414283/1‑251 (MQ=255)
                                                                                                                                                         cTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGa                                                                                    >  1:484003/1‑250 (MQ=255)
                                                                                                                                                         cTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAACCCGGGACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGG‑CGAAAGAGCTGTGGCTAACCGGC‑AAGAGGCCGCCTAACGCCTGTGACCAGTTATTTTGTGTCTGGCGTGCGTCAGGGAAAAAAGATACGCTTGCGTATTTAGAGCGTGTCCGTCTTGCGATg                                                                                  >  1:397503/1‑251 (MQ=255)
                                                                                                                                                          gcgcgcTGCGTGAAGACTGGCGTGGCTGGTTGGCCGTTAGCCCGGAAAAGCCCGGAAGTACCGAAGCGCA‑AT‑GGAATGTCTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGATCCAGCCTAGCGCCTGTGACAAGTTATTTAGTGTCTGGCGGGCGGCAGGTAAACAAGATCCGCTGTAGTATTTAGAGCGTATCCGTCTGGTGATg                                                                                  <  2:318062/249‑1 (MQ=255)
                                                                                                                                                            gcgcGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGaa                                                                                <  1:184589/251‑1 (MQ=255)
                                                                                                                                                                ggcgTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCg                                                                            <  1:499833/251‑1 (MQ=255)
                                                                                                                                                                     gAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCA‑AT‑GGAATTACTACTATGCAAAATGGAACACCGGGCAGAGAGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCCACCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGGAAACGGGTa                                                                        >  1:74645/1‑250 (MQ=255)
                                                                                                                                                                       aGACTGGCGTGGCTTTTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGt                                                                                                                               <  2:279982/193‑1 (MQ=255)
                                                                                                                                                                       aGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGt                                                                                                                               >  1:279982/1‑193 (MQ=255)
                                                                                                                                                                        gACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAaca‑c                                                                   <  1:569744/251‑1 (MQ=255)
                                                                                                                                                                              ggttgCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCt                                                             <  2:224920/247‑1 (MQ=255)
                                                                                                                                                                                  gCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCC‑AT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTgga                                                          >  2:528498/1‑249 (MQ=255)
                                                                                                                                                                                     tGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAAc                                                       <  1:533931/250‑1 (MQ=255)
                                                                                                                                                                                          gCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTATGGCGTGCGTCAGGTAAACAAGATCCGATGGCGTATTTAGAACGAATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCATGGTAACAGGGCt                                                 >  2:316262/1‑251 (MQ=255)
                                                                                                                                                                                              ttAGCCCGGAAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTgg                                               <  1:627777/249‑1 (MQ=255)
                                                                                                                                                                                                  cccGGAAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTagag                                                                                                    >  1:100392/1‑193 (MQ=255)
                                                                                                                                                                                                  cccGGAAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAAAACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCGGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTagag                                                                                                    <  2:100392/193‑1 (MQ=255)
                                                                                                                                                                                                   ccGGAAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATGACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTggcagggca                                        <  1:430180/251‑1 (MQ=255)
                                                                                                                                                                                                         acgCCCGGAACTACCGAAGCGCA‑AT‑GTATTTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGACGCCTGGCAA‑GGGGCGAAAGAGCTGTAGCTAACCGGC‑AAGAGCACGCCTAACGCCTGTGCCAAGTTATTTAGTGTCTGGCGTCCGTCAGGTAAACAAGATCCGCTGGCGTATTTTGAGCGTATCCGTCTGGCGATGAAAGCGGGAAACA‑CTGCCCTGGATACAGTGCTggcagggcagtcgcc                                  >  2:267762/3‑246 (MQ=255)
                                                                                                                                                                                                         aaGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGAtgcc                                  >  1:270930/1‑251 (MQ=255)
                                                                                                                                                                                                          aacccGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTCGTGTCAGGCATGCGTCAGGTAAACCAGATCCGCTGGCGTATTTCGAGCGTATCCGGCTGGCGATGAAAGCGGGTA‑CAGCTGGCCGTGTAACGTTGCTGGCAGGGCCGAtgcct                                 >  2:243136/3‑251 (MQ=255)
                                                                                                                                                                                                            cccGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGAtgcctgc                               <  1:59492/251‑1 (MQ=255)
                                                                                                                                                                                                                  ccTACCGAACCGTA‑AT‑GTAATTACTACTATGCAAAATGGAACACCGGGCATAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCGAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTa                         <  2:559521/250‑1 (MQ=255)
                                                                                                                                                                                                                      ccGAAGCGCA‑AT‑GTAATTACTACTAGGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAg                     <  1:651609/251‑1 (MQ=255)
                                                                                                                                                                                                                          aGCGCA‑AT‑GTAATTACTACTCTGCAAAATGGAACACCGGGCAGAGTGAAGACGCCTAGCAA‑GAGGCGAAAGAGCTGAGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGCATGGCGTGCGTCAGGAAAACAAGATACGCTGGCGTATTTAGAGCGTAACCGTCTGGCGATAAAATCGGGTAACA‑CTGGCCTGGCAACAGTGCTGGCAGGGGAAGTGACTGCCGATTACCacacg                  >  2:206521/1‑246 (MQ=255)
                                                                                                                                                                                                                          aGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACt                  >  1:512929/1‑250 (MQ=255)
                                                                                                                                                                                                                          aGCGCA‑AT‑GTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTATGGCTAACCGGC‑AAGAGCACGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGAGGAGGCCTGCCGATTACCAGACtt                 >  2:86515/1‑250 (MQ=255)
                                                                                                                                                                                                                                             actaTGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCa  <  1:287449/249‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
GGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCA‑AT‑GTAATTACTACTATGCGAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAA‑GGGGCGAAAGAGCTGTGGCTAACCGGC‑AAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGCGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACA‑CAGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCA  >  NC_000913/4630904‑4631387

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