Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 547,879 2 bp→GA coding (595‑596/630 nt) h16_A0522 ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM260479547,8790AG100.0% 132.9 / NA 35L199S (CTC→TCC) h16_A0522conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/16);  total (19/16)
*AM260479547,8800GA100.0% 131.0 / NA 35L199S (CTC→TCC) h16_A0522conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/16);  total (19/16)

AACAGTTCGATCAGCAGCGAGTCGGTCATGGGTTGCGACATGAATCGTATCAACAGGTTGTAGCGCGCGCCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGAGCAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCATTGCGCGGCCGCCGCGGCGCTGGCCGCGCCGTTGCTGCCGCCGCGGCTGCCTGCAGCCGGGAAGCGTGCGTC  >  AM260479/547768‑548003
                                                                                                                ||                                                                                                                           
aaCAGTTCGATCAGCAGCGAGTCGGTCATGGGTTGCGACATGAATCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTg                                                                                                        <  2:751709/134‑1 (MQ=255)
         aTCAGCAGCGAGTCGGTCATGGGTTGCGACATGAATCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGc                                                                                                               <  2:197077/118‑1 (MQ=255)
         aTCAGCAGCGAGTCGGTCATGGGTTGCGACATGAATCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGc                                                                                                               <  2:372238/118‑1 (MQ=255)
         aTCAGCAGCGAGTCGGTCATGGGTTGCGACATGAATCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGc                                                                                                               >  1:197077/1‑118 (MQ=255)
         aTCAGCAGCGAGTCGGTCATGGGTTGCGACATGAATCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGc                                                                                                               >  1:372238/1‑118 (MQ=255)
         aTCAGCAGCGAGTCGGTCATGGGTTGCGACATGAATCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGaa                                                                                               >  1:709022/1‑134 (MQ=255)
         aTCAGCAGCGAGTCGGTCATGGGTTGCGACATGAATCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGaa                                                                                               >  1:599297/1‑134 (MQ=255)
            agcagcGAGTCGGTCATGGGTTGCGACATGAATCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCg                                                                                                             <  2:1017672/117‑1 (MQ=255)
            agcagcGAGTCGGTCATGGGTTGCGACATGAATCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCg                                                                                                             >  1:1017672/1‑117 (MQ=255)
             gcagcGAGTCGGTCATGGGTTGCGACATGAATCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGc                                                                                           <  2:339313/134‑1 (MQ=255)
                  gAGTCGGTCATGGGTTGCGACATGAATCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCgt                                                                                      <  2:599297/134‑1 (MQ=255)
                       ggTCATGGGTTGCGACATGAATCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTg                                                                                            <  1:431692/123‑1 (MQ=255)
                       ggTCATGGGTTGCGACATGAATCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTg                                                                                            >  2:431692/1‑123 (MQ=255)
                                      cATGAATCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCca                                                                                                                    <  1:687458/84‑1 (MQ=255)
                                      cATGAATCGTATCCAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCca                                                                                                                    >  2:687458/1‑84 (MQ=255)
                                              gtatccAACAGGTTGTAGCGCGCGCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCATTGCGCGGCCGCcgcg                                                          >  2:573234/6‑134 (MQ=255)
                                                                    gcgccGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCATTGCGCGGCCGCCGCGGCGCTGGCCGCGCCGTtgctg                                     >  2:469071/4‑134 (MQ=255)
                                                                    gcgccGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCATTGCGCGGCCGCCGCGGCGCTGGCCGCGCCGTtgctg                                     >  2:987988/4‑134 (MQ=255)
                                                                       ccGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTg                                                                                                        <  1:1055717/64‑1 (MQ=255)
                                                                       ccGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTg                                                                                                        >  2:1055717/1‑64 (MQ=255)
                                                                        cGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCATTGCGCGGCCGCCGCGGCGCTGGCCGCGCCGTTGCTgccgc                                 >  1:957167/1‑134 (MQ=255)
                                                                         gCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCATTGCGCGGCCGCCGCGGCGCTGGCCGCGCCGTTGCTgccgcc                                >  1:713672/1‑134 (MQ=255)
                                                                         gCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCATTGCGCGGCCGCCGCGGCGCTGGCCGCGCCGTTGCTgccgcc                                <  2:713672/134‑1 (MQ=255)
                                                                                  cagcGCGGCGAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCATTGCGCGGCCGCCGCGGCGCTGGCCGCGCCGTTGCTGCCGCCGCGGctgcc                       <  2:709022/134‑1 (MQ=255)
                                                                                      gcggcgAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCa                                                                                          >  1:844733/1‑63 (MQ=255)
                                                                                      gcggcgAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCa                                                                                          <  2:844733/63‑1 (MQ=255)
                                                                                        ggcgAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCATTGCGCGGCCGCCGCGGCGCTGGCCGCGCCGTTGCTGCCGCCGCGGCTGCCTGCAGc                 <  1:554480/134‑1 (MQ=255)
                                                                                         gcgAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGc                                                                                           >  1:412945/1‑59 (MQ=255)
                                                                                         gcgAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGc                                                                                           <  2:412945/59‑1 (MQ=255)
                                                                                         gcgAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCATTGCGCGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑CAGCGCCGTTGCTGCCGCCGCGGCTGCCTGCAGCCGGGAAGCGTGCGTc  >  1:1119112/1‑134 (MQ=255)
                                                                                           gAAAGATTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCATTGCGCGGCCGCCGCGGCGCTGGCCGCGCCGTTGCTGCCGCCGCGGCTGCCTGCAGCCgg              >  1:629571/1‑134 (MQ=255)
                                                                                                aTTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCATTGCGCGGCCGCCGCGGCGCTGGCCGCGCCGTTGCTGCCGCCGCGGCTGCCTGCAGc                 <  2:699756/126‑1 (MQ=255)
                                                                                                aTTAGGTTCCACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCATTGCGCGGCCGCCGCGGCGCTGGCCGCGCCGTTGCTGCCGCCGCGGCTGCCTGCAGc                 >  1:699756/1‑126 (MQ=255)
                                                                                                        ccACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCgt                                                                                      <  1:1064258/49‑1 (MQ=38)
                                                                                                        ccACTCCGGACAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCgt                                                                                      >  2:1064258/1‑49 (MQ=38)
                                                                                                                ||                                                                                                                           
AACAGTTCGATCAGCAGCGAGTCGGTCATGGGTTGCGACATGAATCGTATCAACAGGTTGTAGCGCGCGCCCGCTTGTCAGCAGCGCGGCGAAAGATTAGGTTCCACTCCGAGCAGCAGCCACAGGCCGTTCTGTTGCCGGAATTGCAGCGTGGTGTAGGCGCATTGCGCGGCCGCCGCGGCGCTGGCCGCGCCGTTGCTGCCGCCGCGGCTGCCTGCAGCCGGGAAGCGTGCGTC  >  AM260479/547768‑548003

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