Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 1,354,645 +C intergenic (+98/‑15) h16_A1248 → / → h16_A1249 predicted signal transduction protein containingCHASE, PAS, GGDEF and EAL domains/Hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604791,354,6431.C94.3% 125.6 / 0.5 35intergenic (+96/‑17)h16_A1248/h16_A1249predicted signal transduction protein containingCHASE, PAS, GGDEF and EAL domains/Hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (18/15);  total (19/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01

ATGCGCGCCTGACTGCCCGTCAGTCCGGTGCGTGCTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTCCTCAAGCTGGCGGCGAA  >  AM260479/1354536‑1354760
                                                                                                            |                                                                                                                      
aTGCGCGCCTGACTGCCCGTCAGTCCGGTGCGTGCTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGcc                                                                                              <  2:779806/134‑1 (MQ=255)
                      gTCCGGTGCGTGCTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACttt                                                                                                             <  1:1021043/97‑1 (MQ=255)
                      gTCCGGTGCGTGCTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACttt                                                                                                             >  2:1021043/1‑97 (MQ=255)
                                                         gCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAAt                                                      >  2:406756/1‑117 (MQ=255)
                                                         gCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAAt                                                      <  1:406756/117‑1 (MQ=255)
                                                            aGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAg                                                                                           <  2:632810/77‑1 (MQ=255)
                                                            aGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAg                                                                                           >  1:632810/1‑77 (MQ=255)
                                                            aGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGc                                                            >  1:1013270/1‑108 (MQ=255)
                                                            aGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGc                                                            <  2:1013270/108‑1 (MQ=255)
                                                              aTCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCgg                                >  2:158136/1‑134 (MQ=255)
                                                                      ggACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATgcgc                                                                     <  2:1109641/89‑1 (MQ=255)
                                                                      ggACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATgcgc                                                                     >  1:1109641/1‑89 (MQ=255)
                                                                      ggACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCt                                                  >  1:1101723/1‑108 (MQ=255)
                                                                      ggACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCt                                                  <  2:1101723/108‑1 (MQ=255)
                                                                      ggACCGTACTATTGAAATCATGTGGCAGGCCATCGACGccc                                                                                                                     <  1:807784/41‑2 (MQ=255)
                                                                      ggACCGTACTATTGAAATCATGTGGCAGGCCATCGACGccc                                                                                                                     >  2:807784/1‑40 (MQ=255)
                                                                                aTTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAg                                                                                           >  1:924224/1‑57 (MQ=39)
                                                                                aTTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAg                                                                                           <  2:924224/57‑1 (MQ=39)
                                                                                aTTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCg                                                             <  1:535914/87‑1 (MQ=255)
                                                                                aTTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCg                                                             >  2:535914/1‑87 (MQ=255)
                                                                                aTTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCt                                                  >  1:600244/1‑98 (MQ=255)
                                                                                aTTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCt                                                  <  2:600244/98‑1 (MQ=255)
                                                                                aTTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAg                                           >  1:101953/1‑105 (MQ=255)
                                                                                aTTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAg                                           <  2:101953/105‑1 (MQ=255)
                                                                                aTTGAAATCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTCCTc                >  1:254326/1‑132 (MQ=255)
                                                                                       tCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAg                                                                                           <  2:240281/50‑1 (MQ=38)
                                                                                       tCATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAg                                                                                           >  1:240281/1‑50 (MQ=38)
                                                                                        cATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAg                                                                               <  1:949392/61‑1 (MQ=255)
                                                                                        cATGTGGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAg                                                                               >  2:949392/1‑61 (MQ=255)
                                                                                          tgtgGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCGGGCGGAGAGCGGCTACCGGCAATTCCTCAAGCTggcggc     >  2:660544/1‑133 (MQ=255)
                                                                                          tgtgGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTCCTCAAGCTggcggc     >  1:676871/1‑133 (MQ=255)
                                                                                          tgtgGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTCCTCAAGCTggcggc     >  1:360835/1‑133 (MQ=255)
                                                                                            tgGCAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTCCTCAAGCTGGCGGCGaa  <  2:360835/134‑1 (MQ=255)
                                                                                               cAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTCCTCAAGCTGgcggcg    >  2:1146727/1‑129 (MQ=255)
                                                                                               cAGGCCATCGACGCCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTCCTCAAGCTGgcggcg    <  1:1146727/129‑1 (MQ=255)
                                                                                                  gCCATCGACGCCATGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTAcc                           >  1:835331/1‑103 (MQ=255)
                                                                                                  gCCATCGACGCCATGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTAcc                           <  2:835331/103‑1 (MQ=255)
                                                                                                            |                                                                                                                      
ATGCGCGCCTGACTGCCCGTCAGTCCGGTGCGTGCTTAAGCAGTCCAATAACAGGGAGCTAGATCATGCAGGACCGTACTATTGAAATCATGTGGCAGGCCATCGACGCCTGGACTTTTCCAGGATGAAGGCCAAGCTGCTGCATCAGAAGCATGCGCATTGGTCGCCGGAATCGCTGGAGCAGGCGGAGAGCGGCTACCGGCAATTCCTCAAGCTGGCGGCGAA  >  AM260479/1354536‑1354760

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