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→GCAsltsoluble 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

GGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCGAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGCGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACAGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCA  >  NC_000913/4630904‑4631387
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
ggACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCCAAATggg                                                                                                                                                                                                                                               >  1:603647/1‑249 (MQ=255)
   cTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTCACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGCACCGGCGCGGCGTGAAGACTGGCGTTGCTTGTTTACCTTTAGCCCGGGCAAGCCCGGGACCACCGCAGAGCAATGTACTTACTAATATTCAAAACGGAacac                                                                                                                                                                                                                                           >  2:651609/1‑251 (MQ=255)
     gaaagATTATCCACTTTATCCCTACCTGGAATCCCGTCAGATCACCGATGACCTTATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTGCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCTGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACAcc                                                                                                                                                                                                                                          <  1:471997/246‑1 (MQ=255)
      aaGGATTATCCCCCTTATCCCTGCCTGGAATTCCGTCAGATCCCCGATGATCTGATGAAGCAACCGGCGGTGACTGCCACTAACTTTGTTCGCTCTAACCCTACGCTTCCTCCCGCTCGCACGCTGCAATCTCGCTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTTGTCTTTAGCCCGGAAAAGCCAGGAACTGCCGAAGCGCAATGTAATTACTACTATGCCAAATggtccaccgt                                                                                                                                                                                                                                        <  2:1980/251‑9 (MQ=255)
        agagtctCCGCTGTATCTCTACCTGGAATACCGACATATCACCGATGACCTGATGAATCATTCGGCGGTGACTGGCACTAACCTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCATTCTCGTTTCGTCAAGGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGTCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGc                                                                                                                                                                                                                                      <  2:181659/245‑1 (MQ=255)
         gATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTTACCTTTAGACCGGAAAAGACCGGGACTAACGAAGAGCACTTGTATTAATACTATGCAAAAAGGAaccacgggc                                                                                                                                                                                                                                      >  1:410801/1‑243 (MQ=255)
              tccacTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGGGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGGAGTACAGAAGCGCAATTTAATTACTACTATGCCAAATTGAACAACCGGCAGGGTg                                                                                                                                                                                                                                >  1:488034/5‑251 (MQ=255)
                             ccTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATATAATTACTACTATGCAAAATGGAACACCGGGCAGAGGGAAAAAACCTGGCa                                                                                                                                                                                                                   >  1:145881/1‑249 (MQ=255)
                                ggAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAgggg                                                                                                                                                                                                             <  1:673734/251‑1 (MQ=255)
                                   aTACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCg                                                                                                                                                                                                           >  1:165437/1‑250 (MQ=255)
                                             atcaacgaagaccTGATGAATCAACCGGCGGTGACTGTCACTAACTGTGTTCGCGCTAACCCCACGCTGCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCtgtg                                                                                                                                                                                                <  2:488307/239‑1 (MQ=255)
                                              tCACCGATGACCTGATGAATCAACCGGCGGTTACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCACTCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGAGTGGCTTGTTTGTATTTAGCCCGGGAAAACCCGGGAATACCGAAACGCAATTGCAATAATAATATGCAAAAAGGAAAACCGGGCAGAGGGAAGAAACCCGGCAAAGGGGGGAAaggctgagg                                                                                                                                                                                                >  2:645864/1‑244 (MQ=255)
                                                    atgacCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAAcc                                                                                                                                                                                         <  1:82426/251‑1 (MQ=255)
                                                      gaccTGATGAATCAACCGGCTGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCGTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGGGCTGTGGCTAACCg                                                                                                                                                                                        <  2:627631/247‑1 (MQ=255)
                                                       accTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGGAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAAACggg                                                                                                                                                                                      >  1:566074/3‑250 (MQ=255)
                                                               aaTCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAACCTTTAGCCCGGAAAAACCCGGGACTACCGAAGCGCAATGTAATTACTACTATTCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCGGTGGCTAAACGGCAAGagccc                                                                                                                                                                             >  1:368735/1‑250 (MQ=255)
                                                                  cAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGagccagcc                                                                                                                                                                          <  1:532176/251‑1 (MQ=255)
                                                                      cggcggTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGGAAAACCCGGAACTACCGCAGCGCAATGTAATTACTTATCTGCAAAATGGAACACCGGGCAGAGTGGAGAAGCCTGGGAAGGGGGGAAAGAGCTGTGGCTAACCGGCAAAAGGCAGCCTaa                                                                                                                                                                       >  1:391693/1‑250 (MQ=255)
                                                                      cggcggTGACGGTCCCTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTTAAGACTGGCGTGGCTTGTTAGCCGTTAGCCCGGAAAAGCCCGGACCTACCGAAGCGCAATGTAATTACTACTATGCATAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAAc                                                                                                                                                                      <  1:267581/251‑1 (MQ=255)
                                                                          ggTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCGGCTCGCACGCTGCAATCTCGATTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTTGCCCGGAGAAACCCGGAAAAAACGAAGCGCAATGTGATTACTACCATGCAACAAGGAAAAACAGGCAGAGGGAAGACAACTTGCAAGGGGCGAAAGAGCTAGGGCCCACCGGCCAGagccaggccagagcac                                                                                                                                                                  >  2:122770/1‑241 (MQ=255)
                                                                          ggTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCCAAATGGAACACCGGGCAGAGTGGAGAAGCCTGGCAAGGGGCGGAAGAGCTGTGGCTTACCGGCAAGAGCCAGCCTAACGccc                                                                                                                                                                  >  1:237614/1‑250 (MQ=255)
                                                                             gactgTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAAt                                                                                                                                                                                                                                                  <  2:532789/166‑1 (MQ=255)
                                                                             gactgTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAAt                                                                                                                                                                                                                                                  >  1:532789/5‑170 (MQ=255)
                                                                                 gtcgcTAACTTTGTTCGCGTTAACCCCCCGCTTCCTCCCTCTCTCTCGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGTCGTGAAGACTGGCGTGGCTTGATAGTCTTTAGCCCTGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGTCAAGAGACCGCCTGACGTCTTTGACa                                                                                                                                                            <  2:122978/246‑1 (MQ=255)
                                                                                  tCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACa                                                                                                                                                            >  1:220780/1‑249 (MQ=255)
                                                                                    aCTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACGGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGtt                                                                                                                                                        <  2:603647/251‑1 (MQ=255)
                                                                                       aaCTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTAtt                                                                                                                                                     <  1:489045/251‑1 (MQ=255)
                                                                                       aaCTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTAt                                                                                                                                                      <  1:443889/250‑1 (MQ=255)
                                                                                            tGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACCCCGGGCCGAGTGAAGAAGCCTGGGAAGGGGCGAAAGAGCTGGGGCTAACCGGCCAGAGCCAGCCTAACGCCTGTGACAAGTTATTTattg                                                                                                                                                >  2:499833/1‑248 (MQ=255)
                                                                                                cgcgCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTg                                                                                                                                            <  1:29728/251‑1 (MQ=255)
                                                                                                  cgcTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTggg                                                                                                                                          >  1:29197/1‑250 (MQ=255)
                                                                                                               cTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTa                                                                                                                              >  1:428555/1‑250 (MQ=255)
                                                                                                                           cGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGGGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCt                                                                                                                 >  1:129965/1‑251 (MQ=255)
                                                                                                                             cACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGCACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGCACACCGGGCCGAGTTAAGAAGCCTGGCAAGGGGCGAAAGAACTGTGGCTAACCGGCCAGAGCCAGCCTAACGCCTGTGACCAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTgg                                                                                                               >  1:415633/1‑251 (MQ=255)
                                                                                                                                  tGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTTGGCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAAATACTACTATGCAAAATGGAACACCGGGCAAAGAGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGAGGCTAACCGGCCAAGGCCAGCCTAACGCCTTGGACAAGTTTTTTAGTGTCTGGCGTGGGGCCGGTAAAAAAAATACGCTGGCGTAt                                                                                                          >  2:184589/1‑251 (MQ=255)
                                                                                                                                    aaagccagttgcGTCAAGGAACCGGCGCGGCGTGAAGACTGGCTGGGCTGGGTCGCCTTTAGCCCGTAAAAGCCCGGAACTACCGAAGCGCAATTGTATTCCTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCACGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGAGTGCGTCAGGTAAACAAGATCCGCttgcttagtt                                                                                                        <  2:223674/240‑10 (MQ=255)
                                                                                                                                      aTCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACCAGATCCGCTGGCGTAttt                                                                                                        >  1:224920/1‑249 (MQ=255)
                                                                                                                                           gTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGAAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGGGTGAAAAAGGCTGGCAAGGGGAGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGGGACAAGTTATTTAGTGTCTGGAGTGCGTCAGGTCAACAAGATCCGCTGGCGTGTTTAAAGCGt                                                                                                 >  1:403164/1‑251 (MQ=255)
                                                                                                                                                 tCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGATAGCCTTTTGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCCAAATGGCACACCGGGCAGAGTGAAGAAGCCAGGCAAGGGGCGAAAGAGATGTGGCTAACCGGCAAGAGCCCGCCTAACGCCGGTGACCAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGt                                                                                           >  1:414283/1‑251 (MQ=255)
                                                                                                                                                         cTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGa                                                                                    >  1:484003/1‑250 (MQ=255)
                                                                                                                                                         cTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAACCCGGGACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGGCCGCCTAACGCCTGTGACCAGTTATTTTGTGTCTGGCGTGCGTCAGGGAAAAAAGATACGCTTGCGTATTTAGAGCGTGTCCGTCTTGCGATg                                                                                  >  1:397503/1‑251 (MQ=255)
                                                                                                                                                          gcgcgcTGCGTGAAGACTGGCGTGGCTGGTTGGCCGTTAGCCCGGAAAAGCCCGGAAGTACCGAAGCGCAATGGAATGTCTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGATCCAGCCTAGCGCCTGTGACAAGTTATTTAGTGTCTGGCGGGCGGCAGGTAAACAAGATCCGCTGTAGTATTTAGAGCGTATCCGTCTGGTGATg                                                                                  <  2:318062/249‑1 (MQ=255)
                                                                                                                                                            gcgcGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGaa                                                                                <  1:184589/251‑1 (MQ=255)
                                                                                                                                                                ggcgTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCg                                                                            <  1:499833/251‑1 (MQ=255)
                                                                                                                                                                     gAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGGAATTACTACTATGCAAAATGGAACACCGGGCAGAGAGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCCACCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGGAAACGGGTa                                                                        >  1:74645/1‑250 (MQ=255)
                                                                                                                                                                       aGACTGGCGTGGCTTTTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGt                                                                                                                               <  2:279982/193‑1 (MQ=255)
                                                                                                                                                                       aGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGt                                                                                                                               >  1:279982/1‑193 (MQ=255)
                                                                                                                                                                        gACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAacac                                                                   <  1:569744/251‑1 (MQ=255)
                                                                                                                                                                              ggttgCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCt                                                             <  2:224920/247‑1 (MQ=255)
                                                                                                                                                                                  gCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCCATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTgga                                                          >  2:528498/1‑249 (MQ=255)
                                                                                                                                                                                     tGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAAc                                                       <  1:533931/250‑1 (MQ=255)
                                                                                                                                                                                          gCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTATGGCGTGCGTCAGGTAAACAAGATCCGATGGCGTATTTAGAACGAATCCGTCTGGCGATGAAAGCGGGTAACACTGGCATGGTAACAGGGCt                                                 >  2:316262/1‑251 (MQ=255)
                                                                                                                                                                                              ttAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTgg                                               <  1:627777/249‑1 (MQ=255)
                                                                                                                                                                                                  cccGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTagag                                                                                                    >  1:100392/1‑193 (MQ=255)
                                                                                                                                                                                                  cccGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAAAACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCGGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTagag                                                                                                    <  2:100392/193‑1 (MQ=255)
                                                                                                                                                                                                   ccGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATGACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTggcagggca                                        <  1:430180/251‑1 (MQ=255)
                                                                                                                                                                                                         acgCCCGGAACTACCGAAGCGCAATGTATTTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGACGCCTGGCAAGGGGCGAAAGAGCTGTAGCTAACCGGCAAGAGCACGCCTAACGCCTGTGCCAAGTTATTTAGTGTCTGGCGTCCGTCAGGTAAACAAGATCCGCTGGCGTATTTTGAGCGTATCCGTCTGGCGATGAAAGCGGGAAACACTGCCCTGGATACAGTGCTggcagggcagtcgcc                                  >  2:267762/3‑246 (MQ=255)
                                                                                                                                                                                                         aaGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGAtgcc                                  >  1:270930/1‑251 (MQ=255)
                                                                                                                                                                                                          aacccGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTCGTGTCAGGCATGCGTCAGGTAAACCAGATCCGCTGGCGTATTTCGAGCGTATCCGGCTGGCGATGAAAGCGGGTACAGCTGGCCGTGTAACGTTGCTGGCAGGGCCGAtgcct                                 >  2:243136/3‑251 (MQ=255)
                                                                                                                                                                                                            cccGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGAtgcctgc                               <  1:59492/251‑1 (MQ=255)
                                                                                                                                                                                                                  ccTACCGAACCGTAATGTAATTACTACTATGCAAAATGGAACACCGGGCATAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCGAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTa                         <  2:559521/250‑1 (MQ=255)
                                                                                                                                                                                                                      ccGAAGCGCAATGTAATTACTACTAGGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAg                     <  1:651609/251‑1 (MQ=255)
                                                                                                                                                                                                                          aGCGCAATGTAATTACTACTCTGCAAAATGGAACACCGGGCAGAGTGAAGACGCCTAGCAAGAGGCGAAAGAGCTGAGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGCATGGCGTGCGTCAGGAAAACAAGATACGCTGGCGTATTTAGAGCGTAACCGTCTGGCGATAAAATCGGGTAACACTGGCCTGGCAACAGTGCTGGCAGGGGAAGTGACTGCCGATTACCacacg                  >  2:206521/1‑246 (MQ=255)
                                                                                                                                                                                                                          aGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACt                  >  1:512929/1‑250 (MQ=255)
                                                                                                                                                                                                                          aGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTATGGCTAACCGGCAAGAGCACGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGAGGAGGCCTGCCGATTACCAGACtt                 >  2:86515/1‑250 (MQ=255)
                                                                                                                                                                                                                                             actaTGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCa  <  1:287449/249‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
GGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCGAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGCGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACAGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTACCAGACTATCGCCTCGGCAATCA  >  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: