Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 396,286 A→G 6.4% M1S (TTG→TCG) † ampH ← D‑alanyl‑D‑alanine‑ carboxypeptidase/endopeptidase; penicillin‑binding protein; weak beta‑lactamase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913396,2860AG 6.4% 195.1 / 3.8 63M1S (TTG→TCG) †ampHD‑alanyl‑D‑alanine‑ carboxypeptidase/endopeptidase; penicillin‑binding protein; weak beta‑lactamase
Reads supporting (aligned to +/- strand):  ref base A (20/39);  new base G (2/2);  total (22/41)
Fisher's exact test for biased strand distribution p-value = 6.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.23e-01

TCATTAGCTTGGTGAGGGAAGCAATACGCACGACCGAATCCAGCTGCGGGCGAACGTTATTACCAGGTCGCGTTTCGCCATAACTGCGAAAGACGCGCTGGTTACCGTCGATAACAACCAGAGCCATCCCCGTGGCACCGCTGCCGTAAAAAATATGATCGGCATAACGATCGACAATATCAGAGGCAAACTCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCAACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAGCGTATTTATACTACTCATTGATAGATTGCAAAGGGCGAATTAGTGCAATAATAGTGTTAAAAACGTAACGGCGCGTAAACACGATCTTTTTCCTGCTCACCAATGCCTGCATTTGTGATCCTGCGCGTGTTAATCCGCTTATGTAAACAATCTTTGGTTTGTGAAATCCGTCGAAACAATTCTTATGGTCAGGCTGGCGAACTAAGCGCCTTGCTATGGGTCACAATGGGCGGTTTCA  >  NC_000913/396011‑396560
                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                   
aaatgagcggtgagAGGGAAAAAAAACGCACNCCCAAATCCAGCTGCGGGAGAACGTNATTACCAGGTCGCGTGTCGACATAACCGCGAAAGACGCGCTGGTTAACATCGATAAAAACCAGAGCCATCCCCGTGGCACNGCTGGCGTAAAAAATATGATCGGCATAACGATCGACAATATCAGAGGCAAACTCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCaa                                                                                                                                                                                                                                                                                    <  3:102032/264‑1 (MQ=255)
    tAGCTTGGCGAGGGAAGCAATACGCACGACCGAATCCAGCTGCGGGCGAACGTTATTACCAGGTCGCGTTTCGCCATAACTGCGAAAGACGCGCTGGTTACCGTCGATAACAACCAGAGCCATCCCCGTGGCACCGCTGCCGTAAAAAATATGATCGGCATAACGATCGACAATATCAGAGGCAAACTCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGCAAAAAGCAGACTACGTTTCAACg                                                                                                                                                                                                                                                                                  >  2:120124/1‑275 (MQ=255)
      gCTTGGTGAGGGAAGCAATACGCACGACCGAATCCAGCTGCGGGCGAACGTTATTACCAGGTCGCGTTTCGCCATAACTGCGAAAGACGCGCTGGTTACCGTCGATAACAACCAGAGCCATCCCCGTGGCACCGCTGCCGTAAAAAATATGATCGGCATAACGATCGACAATATCAGAGGCAAACTCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCAACggt                                                                                                                                                                                                                                                                                <  4:409535/275‑1 (MQ=255)
       cTTGGTGAGGGAAGCAATACGCACGACCGAATCCAGCTGCGGGCGAACGTTATTACCAGGTCGCGTTTCGCCATAACTGCGAAAGACGCGCTGGTTACCGTCGATAACAACCAGAGCCATCACAGTGGCACCGCTGCCNTAAAAAATATGATCGGCATAACGATCGACAATATCAGAGGCAAACTACGGTTCGGTGATCGGCTGCGACGCATTNAAGAGGGGCAATGACGCCGCACAACGCACCNNAGACAAAAGCAGACCACGCTTAAACcgcggtng                                                                                                                                                                                                                                                                           >  3:73768/1‑277 (MQ=255)
         tGGTGAGGGAAGCAATACGCACGACCGAATCCAGCTGCGGGCGAACGTTATTACCAGGTCGCGTTTCGCCATAACTGCGAAAGACGCGCTGGTTACCGTCGACAACAACCACAGCCATCCCCGTGACGCCGCGGCCGTAAAAAATCTGAACGGCCTACCTAACGACAATATCAGATGCAAACTCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAAGAACGCAGCACCCACAACGGCAGACAAAAACAGACTACGTTTCAACggtggt                                                                                                                                                                                                                                                                             >  4:382563/1‑275 (MQ=255)
               gggAAGCAATACGCACGACCGAATCCAGCTGCGGGCGAACGTTATTACCAGGTCGCGTTTCGCCATAACTGCGAAAGACGCGCTGGTTACCGTCGATAACAACCAGAGCCATCCCCGTGGCACCGCTGCCGTAAAAAATATGATCGGCATAACGATCGACAATATCAGAGGCAAACTCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCAACGGTGGTGTCCa                                                                                                                                                                                                                                                                        <  2:510895/274‑1 (MQ=255)
                     cAATCCGAACGACCGAAGCCGGATGGGGGCGANCGTTATTCCCAGTTCGCTTGTCACCNAAACTGAGAAAGACGCGCTGGTTACCGTCGAGAACAACCAGAGCCATCCCCGTGGCACCGCTGCCGTAAAAAATATGATGNGCATAACGATCGACAATATCAGAGGCAAACGCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCAACGGTGGTGTCCATAGATtgaat                                                                                                                                                                                                                                                              <  3:102517/278‑1 (MQ=255)
                     cAATACGCACGACCGAATCCAGCTGCGGGCGAACGTTATTACCAGGTCGCGTTTCGCCATAACTGCGAAAGACGCGCTGGTTACCGTCGATAACAACCAGAGCCATCCCCGTGGCACCGCTGCCGTAAAAAATATGATCGGCATAACGATCGACAATATCAGAGGCAAACTCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCAACGGTGGTGTCCATAGATtg                                                                                                                                                                                                                                                                 >  4:514009/1‑275 (MQ=255)
                                 ccGAATCCAGCTGCGGGCTAACGTTATTCCCAGGTCGCGTTTCGCCATAACTGCGAAAGACGCGCTGGTTACCGTCGATAACAACCAGAGCCATCCCCGTGGCACCGCTGCCGTAAAAAATATGATCGGCATAACGATCGACAATATCAGAGGCAAACTCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCAACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAg                                                                                                                                                                                                                                                 <  1:41208/279‑1 (MQ=255)
                                    aaTCCAGCTGCGGGCGAACCTTATTACCAGGTCGCGTTTCGCCATAACTGCGAAAGACGCTCTGGTTACCGTCGATAACAACCAGAGCCATCCCCGTGGCACCGCTGCCGTAAAAAATATGATCGGCATAACGATCGACAATATCAGAGGCAAACTCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCAACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAGCGt                                                                                                                                                                                                                                              <  1:73768/279‑1 (MQ=255)
                                     gtccgtatGCGGGCGGACCGTATGACCAGGGGTGGTTTCGACATAACTGCGGAAGACGCGCTGGTTACCGTCGATAACAACCAGAGCCATCTCCGTGGCACCGCTGCCGTAAAAAATATGATCGGTATAACGATCGACAATATCAGAGGCAAACTCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCAACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAg                                                                                                                                                                                                                                                 <  4:182344/268‑1 (MQ=255)
                                     aaccggcgtcGGGCGAACGTGTTTCCCAGGTCTCGTTTGGCCATAACTGCGAAAGACGGGCTGGTTACCGTGGATAACAACGAGAGCCATCCCCGTGGCACCGCTGCCGTAAAAAATATGATCGGCATAACGATCGACAATATCAGAGGCAAACTCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAAGGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCAACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAg                                                                                                                                                                                                                                                 <  4:414899/266‑1 (MQ=255)
                                      tCCAGCTGCGGGCGACCGTTATTACCAGGGCGCGTTTCGCCAGAACTGCGAAAGACGCGCTGGTTACCGTCGATAACAACCAGAGCCATCCCCGTGGCACCGCTGCCGTAAAAAATATGATCGGCATAACGATCGACAATATCAGAGGCAAACTCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCAACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAGCGTAt                                                                                                                                                                                                                                            <  1:165250/279‑1 (MQ=255)
                                        cAGCTGCGGGCGAACGTTATTACCAGGTCGCGTTTCGCCATAACTGCGAAAGACGCGCTGGTTACCGTCGATAACAACCAGAGCCATCCCCGTGGCACCGCTGCCGTAAAAAATATGATCGGCATAACGATCGACAATATCAGAGGCAAACTCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCAACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAGCGt                                                                                                                                                                                                                                              <  2:147638/275‑1 (MQ=255)
                                          gCTCCGGGCGAACGTTATTACCAGGTCGCGTTTCGCCATAACTGCGAAAGACGCGCTGGTTACCGTCGATAACAACCAGAGCCATCCCCGTGGCACCGCTGCCGTAAAAAATATGATCGGCATAACGATCGACAATATCAGAGGCAAACTCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCAACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAGCGTAt                                                                                                                                                                                                                                            <  2:27053/275‑1 (MQ=255)
                                           cTGCGGGCGAACGTTATTACCAGGTCGCGTTTCGCCATAACTGCGAAAGACGCGCTGGTTACCGTCGATAACAACCAGAGCCATCCCCGTGGCACCGCTGCCGTAAAAAATATGATCGGCATAACGATCGACAATATCAGAGGCAAACTCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCAACGGTGGTGTCCATAGATTGAATGAAGGt                                                                                                                                                                                                                                                        >  2:330188/1‑262 (MQ=255)
                                           cTGCGGGCGAACGTTATTACCAGGTCGCGTTTCGCCATAACTGCGAAAGACGCGCTGGTTACCGTCGATAACAACCAGAGCCATCCCCGTGGCACCGCTGCCGTAAAAAATATGATCGGCATAACGATCGACAATATCAGAGGCAAACTCCGGTTCAGTGATCGGCTGCGCCGAATGGACGGAGGTCAACGACCCCGCACACAGCCCGGCAGAAAAAAGCAGACGACGTTACAAAGGAGGCGCCCAACGATGGCATGAAGGTAACGGAACATAtt                                                                                                                                                                                                                                           >  4:8880/1‑275 (MQ=255)
                                               ggggGAACGTTATTAACAGGTCGCGTTTCGCCATAACTGCGAAAGACGGGCTGGTTACCGTCGATAACAACCAGAGCCATCCCCGTGGCACCGCTGCCGTAAAAAATATGATCGGCATAACGATCGACAATATCAGAGGCAAACTCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCAACGGTGGTGTCCATAGATTGAATGAAGGTAACt                                                                                                                                                                                                                                                    <  4:330188/262‑1 (MQ=255)
                                               gggCGAACGTTATTTCCAGGTCGCGTTTCGCCATAACTGCGAAAGACGCGCTGGTTACCGTCGATAACAACCAGAGCCATCCCCGTGGCACCGCTGCCGTAAAAAATATGATCGGCATAACGATCGACAATATCAGAGGCAAACGCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCAACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAGCGTATTTAtacta                                                                                                                                                                                                                                    <  1:3836/278‑1 (MQ=255)
                                               gggCGAACGTTATTACCAGGTCGCGTTTCGCCATAACTGCGAAAGACGCGCTGGTTACCGTCGATAACAACCAGAGCCATCCCCGTGGCACCGCTGCCGTAAAAAATATGATCGGCATAACGATCGACAATATCAGAGGCAAACTCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCCGCAAAAAGCCGACTCAGTTTCAAAGGTGGTTTCCATAAATGTAATGAAGGTAACTGAACGTATTtata                                                                                                                                                                                                                                       >  4:97188/1‑275 (MQ=255)
                                                ggCGAACGTTATTACCAGGTCGCGTTTCGCCATAACTGCGAAAGACGCGCTGGTTACCGTCGATAACAACCAGAGCCATCCCCGTGGCACCGCTGCCGTAAAAAATATGATCGGCATAACGATCGACAATATCAGAGGCAAACTCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCAACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAGCGTATTTAtac                                                                                                                                                                                                                                      <  2:119961/275‑1 (MQ=255)
                                                  cGAACGTTATTACCAGGTCGCGTTTCGCCATAACTGCGAAAGACGCGCTGGTTACCGTCGATAACAACCAGAGCCATCCCCGTGGCACCGCTGCCGTAAAAAATATGATCGGCATAACGATCGACAATATCAGAGGCAAACTCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCAACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAGCGTATTTAtacta                                                                                                                                                                                                                                    <  2:10387/275‑1 (MQ=255)
                                                               ttgtctacgTTTGGCCACAACTGCCAAAGGATCCCGGGTTACCGTCCATAACAACCAGAGCCATCCCCGTGGAACCCCTGCNGTAAAAAATATGTTTGGCAGAACGAGCGACAATATCAGAGGCAAACTCGGGTCCAGGTAGCGGCTGCGCCGAATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCAACGGTGGTTTCCATAGATTTAATGAAGGTAACTGAGCGTATTTATACTACTCATTGATAGATTGCa                                                                                                                                                                                                                   <  3:131894/272‑1 (MQ=255)
                                                                 ggTCGCGTTTCGCCATAACTGCGAAAGACGCGCTGGTTACCGTCGATAACAACCAGAGCCATCCCCGTGGCACCGCTGCCGTAAAAAATATGATCGGCATAACGATCGACAATATCAGAGGCAAACTCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCAAc                                                                                                                                                                                                                                                                                   >  2:402462/1‑213 (MQ=255)
                                                                   tggcagttCGCGATAACTGCGAAAAGCGGGAAGGGTACGCTGGATAAAAACCAGGGCCATCCCCGTGGAACCGCCGCCGTAAAAAAAATGATCGGCGTAACGATCGACAATATCAGAGGCAAACGCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCAACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAGCGTATTTATACTACTCATTGATAGATTGCa                                                                                                                                                                                                                   <  4:262035/269‑1 (MQ=255)
                                                                   tCGCGTTTCGCCATAACTGCGAAAGACGCGCTGGTTACCGTCGATAACAACCAGAGCCATCCCCGTGGCACCGCTGCCGTAAAAAATATGATCGGCATAACGATCGACAATATCAGAGGCAAACTCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCAACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAGCGTATTTATACTACTCATTGATAGATTGCa                                                                                                                                                                                                                   <  2:514009/275‑1 (MQ=255)
                                                                      cgTTTCGCCATAACTGCGAAAGACGCGCTGGTTACCGTCGATAACAACCAGAGCCATCCCCGTGGCACCGCTGCCGTAAAAAATATGGTCGGCATAACGATCGACAATATCAGAGGCAAACTCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCAACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAGCGt                                                                                                                                                                                                                                              <  2:32304/245‑1 (MQ=255)
                                                                             ccATAACTGCGAAAGACGCGCTGGTTACCGTCGATAACAACCAGAGCCATCCCCGTGGCACCGCTGCCGTAAAAAATATGATCGGCATAACGATCGACAATATCAGAGGCAAACTCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCAACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAGCGTATTTATACTACTCATTGATAGATTGCAAAGGGCGAAt                                                                                                                                                                                                         <  4:359623/275‑1 (MQ=255)
                                                                             ccATAACTGCGAAAGACGCGCTGGTTACCGTCGATAACAACCAGAGCCATCCCCGTGGCACCGCTGCCGTAAAAAATATGAACGGCATAACGATCGACAATATCAGAGGCAAACTCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCAACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAGCGTATTTATACTACTCATTGATAGATTGCAAAGGGCGAAt                                                                                                                                                                                                         <  2:578993/275‑1 (MQ=255)
                                                                                 aaCTGCGAAAGACGCGCTGGTTACCGTCGATAACAACCAGAGCCATCCCCGTGGCACCGCGGCCGTAAACAATATGATCGGCATAACGATCGACAATATCAGAGGCAAACTCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCAACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAGCGTATTTATACTACTCATTGATAGATTGCAAAGGGCGAATTAg                                                                                                                                                                                                      <  4:116516/274‑1 (MQ=255)
                                                                                                                       agagCCATCCCCGTGGCACCGCTGCCGTAAAAAATATGATCGGCATAACGATCGACAATATCAGAGGCAAACTCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCAACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAGCGTATTTATACTACTCATTGATAGATTGCAAAGGGCGAATTAGTGCAATAATAGTGTTAAAAACGTAACGGCGCGTAAacac                                                                                                                                                               <  2:466601/275‑1 (MQ=255)
                                                                                                                       aaggCCATCCCCGTGGCACCGCTGCCGTAAAAAATATGATCGGCATAACGATCGACAATATCAGAGGCAAACTCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCAACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAGCGTATTTATACTACTCATTGATAGATTGCAAAGGGCGAATTAGTGCAATAATAGTGTTAAAAACGTAACGGCGCGTAAacac                                                                                                                                                               <  2:494187/272‑1 (MQ=255)
                                                                                                                                     cgCACCTCGGCCAGAAAAAATATGATCGGCATAACGATCGACAATATCGGATGCAAACTCCGGTTCAGTGATCGGCTTCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCAACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAGCGTATTTATACTACTCATTGATAGATTGCAAAGGGCGAATTAGTGCAATAATAGTGTTAAAAACGTAACGGCGCGTAAACACGATCTTTTTCCTGCTCAc                                                                                                                                             <  1:62775/278‑1 (MQ=255)
                                                                                                                                          cGCTGCCGTAAAAAATATGATCGGCATAACGATCGACAATATCAGAGGCAAACTCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCAACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAGCGTATTTATACTACTCATTGATAGATTGCAAAGGGCg                                                                                                                                                                                                            >  4:69507/1‑211 (MQ=255)
                                                                                                                                            cTGCCGTAAAAAATATGATCGGCATAACGATCGACAATATCAGAGGCAAACTCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGATTACGTTTCAACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAGCGTATTTATACTACTCATTGATAGATTGCAAAGGGc                                                                                                                                                                                                             >  2:203863/1‑208 (MQ=255)
                                                                                                                                              gCCGTAAAAAATATGATCGGCATAACGATCGACAATATCAGAGGCAAACTCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCAACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAGCGTATTTATACTACTCATTGATAGATTGCAAAGGGCGAAtt                                                                                                                                                                                                        <  2:69507/211‑1 (MQ=255)
                                                                                                                                               ccGTAAAAAATATGATCGGCATAACGATCGACAATATCAGAGGCAAACTCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGNTTCAACGGTGGTTTCCATAGATGGAATGAAGGTAACTGAGCGTATTTATACTAATCATTGATAGATTGANAAGGGCGAATTAGTGCATAAATAGTGTTAAAAACGAAAAGGCGCGNANACAAGAAATTTTTCCTGCCACACAATTCCTGc                                                                                                                                  >  3:61701/1‑279 (MQ=255)
                                                                                                                                                cGTAAAAAATATGATCGGCATAACGATCGACAATATCAGAGGCAAACTCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCAACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAGCGTATTTATACTACTCATTGATAGATTGCAAAGGGCGAATTAGTGCAATAATAGTGTTAAAAACGTAACGGCGCGTAAACACGATCTTTTTCCTGCTCACCAATGcc                                                                                                                                     <  2:8880/275‑1 (MQ=255)
                                                                                                                                                      aaaTATGATAGGCATAACGAGCGACAATATCAGAGGCAAACGCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCAACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAGCGTATTTATACTACTCATTGATAGATTGCAAAGGGCGAATTAGTGCAATAATAGTGTTAAAAACGTAACGGCGCGTAAACACGATCTTTTTCCTGCTCACCAATGCCTGCAt                                                                                                                                <  2:382563/274‑1 (MQ=255)
                                                                                                                                                         tatGATCGGCATAACGATCGACAATATCAGAGGCAAACTCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCAACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAGCGTATTTATACTACTCATTGATAGATTGCAAAGGGCGAATTAGTGCAATAATAGTGTTAAAAACGTAACGGCGCGTAAACACGATCTTTTTCCTGCTCACCAATGCCTGCATTtgt                                                                                                                            >  2:327095/1‑275 (MQ=255)
                                                                                                                                                          atGAGCGGCAGAACGGACGGCAATATCAGAGGCAAACTCCGGGTAAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCAGACAGCACGGCAGAAAAAAGCAGACGACGTTTAAACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAGCGTATTTATACTACTCATTGATAGATTGCAAAGGGCGAATTAGTGCAATAATAGTGTTAAAAACGTAACGGAGCGTAAACACGATCTTTTTCATGCTCACCAATGCCTGCATTtgt                                                                                                                            <  2:67185/274‑1 (MQ=255)
                                                                                                                                                                 gCATAACGATCGACAATATCAGAGGCAAACTCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCGACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAGCGTAttt                                                                                                                                                                                                                                          >  4:581880/1‑158 (MQ=255)
                                                                                                                                                                     aaCGATCGACAATATCAGAGGCAAACTCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCGACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAGCGTATTTAtac                                                                                                                                                                                                                                      <  2:581880/158‑1 (MQ=255)
                                                                                                                                                                            gACAATATCAGAGGCAAACTCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCAACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAGCGTATTTATACTACTCATTGATAGATTGCAAAGGGCGAATTAGTGCAATAATAGTGTTAAAAACGTAACGGCGCGTAAACACGATCtttt                                                                                                                                                       <  1:2145/230‑1 (MQ=255)
                                                                                                                                                                            gACAATATCAGAGGCAAACTCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCAACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAGCGTATTTATACTACTCATTGATAGATTGCAAAGGGCGAATTAGTGCAATAATAGTGTTAAAAACGTAACGGCGCGTAAACACGATCTTTTTCCTGCTCACCAATGCCTGCATTTGTGATCCTGCGCGTGTTAATc                                                                                                         >  2:556719/1‑275 (MQ=255)
                                                                                                                                                                                        ggCAAACTCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCAACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAGCGTATTTATACTACTCATTGATAGATTGCAAAGGGCGAATTAGTGCAATAATAGTGTTAAAAACGTAACGGCGCGTAAACACGATCTTTTTCCTGCTCACCAATGCCTGCATTTGTGATCCTGCGCGTGTTAATCCGCTTATGTaaa                                                                                             >  2:225938/1‑275 (MQ=255)
                                                                                                                                                                                         gCAAACTCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCAACGGTGGTGTCCATAGATtgaat                                                                                                                                                                                                                                                              >  1:92927/1‑114 (MQ=255)
                                                                                                                                                                                         gCAAACTCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCAACGGTGGTGTCCATAGATtgaat                                                                                                                                                                                                                                                              <  3:92927/114‑1 (MQ=255)
                                                                                                                                                                                         gCAAACTCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCAACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAGCGTATTTATACTACTCATTGATAGATTGCa                                                                                                                                                                                                                   <  1:141159/157‑1 (MQ=255)
                                                                                                                                                                                         gCAAACTCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCAACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAGCGTATTTATACTACTCATTGATAGATTGCa                                                                                                                                                                                                                   >  3:141159/1‑157 (MQ=255)
                                                                                                                                                                                                         tGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCGACGGTGGTGTCCATAGATTGAATGAAGGTAAc                                                                                                                                                                                                                                                     >  4:359904/1‑107 (MQ=255)
                                                                                                                                                                                                             cGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCGACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAg                                                                                                                                                                                                                                                 <  2:359904/107‑1 (MQ=255)
                                                                                                                                                                                                               gCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTCCGTTTCAACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAGCGTATTTATACTACTCATTGATAGATTGCAAAGGGCGAATTAGTGCAATAATAGTGTTAAAAACGTAACGGCGCGTAAACACGATCTTTTTCCTGCTCACCAATGCCTGCATTTGTGATCCTGCGCGTGTTAATCCGCTTATGTAAACAATCTTTGGTTTGTGAAATCCGTCGa                                                                  <  1:175586/279‑1 (MQ=255)
                                                                                                                                                                                                                   cgccgcATGGACGGCGGTCAATGATGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCAACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAGCGTATTTATACTACTCATTGATAGATTGCAAAGGGCGAATTAGTGCAATAATAGTGTTAAAAACGTAACGGCGCGTAAACACGATCTTTTTCCTGCTCACCAATGCCTGCATTTGTGATCCTGCGCGTGTTAATCCGCTTATGTAAACAATCTTTGGTTTGTGAAATCCGTCGAcaca                                                              >  1:39236/1‑279 (MQ=255)
                                                                                                                                                                                                                             gcggcggGCAAAGAAGCCGTGCACAACAGAGCAGGAGAAAGCAGACCAGGTTTGAAAGGTTGTGTCCATGAAATGAATGAAGGTAACTGGGCGTATATATACTACTGATTGATAGATTGCACAAGGCGACGTAGTGCAATACGAGTGTTAAAAACGTAACGGCGCGTAAACACGATTTTTTTCCTGCGCACCAAATCCTGCAATTGTGATACTGCGCGTGTTAATCCGCTTATGTAGACAATCTATCGTTTGAGAAATCCGTCGAAACAGTTCt                                                         <  4:120124/273‑1 (MQ=255)
                                                                                                                                                                                                                                                      gCACGGCAGAAAAAAGCAGACTACGTTTCAACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAGCGTATTTATACTACTCATTGATAGATTGCAAAGGGCGAATTAGTGCAATAATAGTGTTAAAAACGTAACGGCGCGTAAACACGATCTTTTTCCTGCTCACCAATGCCTGCATTTGTGATCCTgcg                                                                                                                   >  4:408082/1‑191 (MQ=255)
                                                                                                                                                                                                                                                      gCACGGCAGAAAAAAGCAGACTACGTTTCAACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAGCGTATTTATACTACTCATTGATAGATTGCAAAGGGCGAATTAGTGCAATAATAGTGTTAAAAACGTAACGGCGCGTAAACACGATCTTTTTCCTGCTCACCAATGCCTGCATTTGTGATCCTGCGCGTGTTAATCCGCTTATGTAAACAATCTTTGGTTTGTGAAATCCGTCGAAACAATTCTTATGGTCAGGCTGGCGAACTAAgcgc                               <  4:556719/275‑1 (MQ=255)
                                                                                                                                                                                                                                                        acgtaaaaaaaaaTCATTCTACGGTTCAAAGGNGNGGTCAATAGATTGCATGTTGGTGACTGAGTGTTTGTATACGCCTTNTTGATAGAATTCAAAGGGCTAATTAGTGTAATAATAGTGTTAAAACAGTAACGGCGCGTAAACACGATCNTTGTCCTGCTCACCAATGCCTGCATTTGTGATCCTGCGCGTGTTAATCCGCTTATGTAAACAATCTTTGGTTTGTGAAATCCGTCGAAACAATTCTTATGGTCATGCTGGCGAACTAAGCGCCTTGCt                         <  3:39236/272‑1 (MQ=255)
                                                                                                                                                                                                                                                          ggCAGAAAAAAGCAGACTACGTTTCAACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAGCGTATTTATACTACTCATTGATAGATTGCAAAGGGCGAATTAGTGCAATAATAGTGTTAAAAACGTAACGGCGCGTAAACACGATCTTTTTCCTGCTCACCAATGCCTGCATTTGTGATCCTGCGCgtg                                                                                                               <  2:408082/191‑1 (MQ=255)
                                                                                                                                                                                                                                                           gaaGAAAAAAGCAGACTGCTTTTCAACGGTGGTTTCCATAGATTGAATGAAGGTAACTGAGCGTATTTATACTACTCATTGATAGATTGCAAAGGGCGAATTAGTGCAATAATAGTGTTAAAAACGTAACGGCGCGTAAACACGATCTTTTTCCTGCTCACCAATGCCTGCATTTGTGATCCTGCGCGTGTTAATCCGCTTATGTAAACAATCTTTGGTTTGTGAAATCCGTCGAAACAATTCTTATGGTCAGGCTGGCGAACTAAGCGCCTTGCTAt                       <  1:61701/276‑1 (MQ=255)
                                                                                                                                                                                                                                                                      cAGACTACGTTTCAACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAGCGTATTTATACTACTCATTGATAGATTGCAAAGGGCGAATTAg                                                                                                                                                                                                      <  3:74650/93‑1 (MQ=255)
                                                                                                                                                                                                                                                                      cAGACTACGTTTCAACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAGCGTATTTATACTACTCATTGATAGATTGCAAAGGGCGAATTAg                                                                                                                                                                                                      >  1:74650/1‑93 (MQ=255)
                                                                                                                                                                                                                                                                         aCTACGTTTCAACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAGCGTATTTATACTACTCATTGATAGATTGCAAAGGGCGAAtt                                                                                                                                                                                                        <  1:159116/88‑1 (MQ=255)
                                                                                                                                                                                                                                                                         aCTACGTTTCAACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAGCGTATTTATACTACTCATTGATAGATTGCAAAGGGCGAAtt                                                                                                                                                                                                        >  3:159116/1‑88 (MQ=255)
                                                                                                                                                                                                                                                                                   aaCGGTGGTGTCCATAGATTGAATGAAGGTAACTGAGCGTATTTATACTACTCATTGATAGATTGCAAAGGGCGAATTAGTGCAATAATAGTGTTAAAAACGTAACGGCGCGTAAACACGATCTTTTTCCTGCTCACCAATGCCTGCATTTGTGATCCTGCGCGTGTTAATCCGCTTATGTAAACAATCTTTGGTTTGTGAAATCCGTCGAAACAATTCTTATGGTCAGGCTGGCGAACTAAGCGCCTTGCTATGGGTCACAATGGGCGGTTTCa  >  2:466973/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                                   aaCGGTGGTGTCCATAGATTGAATGAAGGTAACTGAGCGTATTTATACTACTAATTGAGAGATTGCAAAGGGCGCATGAGTGCAATAATAGTGTTATAAACGTAACGGCGCGTAAACACGATATTTTTCCTGCTAATCAATGCAGACATTTGTGTTCCTGCGCGTCTTAATCCGCTTATGTAACCAATGTTCGGTTTGACAAATCCGTCGAAACAATTCTTATGGGCAGGCTAGCCGACTGAGGTGTTTGTTCTGGAGCACAATGGGAGGtgtaa  >  4:589217/1‑271 (MQ=255)
                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                   
TCATTAGCTTGGTGAGGGAAGCAATACGCACGACCGAATCCAGCTGCGGGCGAACGTTATTACCAGGTCGCGTTTCGCCATAACTGCGAAAGACGCGCTGGTTACCGTCGATAACAACCAGAGCCATCCCCGTGGCACCGCTGCCGTAAAAAATATGATCGGCATAACGATCGACAATATCAGAGGCAAACTCCGGTTCAGTGATCGGCTGCGCCGCATGGACGGAGGTCAATGACGCCGCACACAGCACGGCAGAAAAAAGCAGACTACGTTTCAACGGTGGTGTCCATAGATTGAATGAAGGTAACTGAGCGTATTTATACTACTCATTGATAGATTGCAAAGGGCGAATTAGTGCAATAATAGTGTTAAAAACGTAACGGCGCGTAAACACGATCTTTTTCCTGCTCACCAATGCCTGCATTTGTGATCCTGCGCGTGTTAATCCGCTTATGTAAACAATCTTTGGTTTGTGAAATCCGTCGAAACAATTCTTATGGTCAGGCTGGCGAACTAAGCGCCTTGCTATGGGTCACAATGGGCGGTTTCA  >  NC_000913/396011‑396560

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