Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260480 274,470 (G)10→9 intergenic (+18/‑626) h16_B0246 → / → fhaB Hypothetical protein/filamentous hemagglutinin / adhesin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM260480274,4610G.93.3% 115.8 / 1.8 30intergenic (+9/‑635)h16_B0246/fhaBHypothetical protein/filamentous hemagglutinin / adhesin
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base . (15/13);  total (16/14)
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.42e-01

GGTGGATCTGGAGCAGGGGATGTGGGAAAGCTTGCTTTCCCTCGCGGCCTGGCTGGCGCATGAGACGCAGGTGTTGAGTGAGCTGGCCGTGCTGCAGGGTCAGGTATTGCAGACTGACTGACGTTTAGCGGGGGGGGGGCAGCCGTGCAATCGGCTGCCGCGCGCCTCCTCGGCGCGCAGGTCAACTTCCCCGATTTCAAAGCGGTTGAACCACACGGCAAAGACGCCATCACGCCGAGTCGGGCGCAGGCGAGAC  >  AM260480/274332‑274587
                                                                                                                                 |                                                                                                                              
ggTGGATCTGGAGGAGGGGATGGGGGGAATCTTGCTTTCCCTCGCGGGCTGGCTGGCGGATGAGAGGCATGGGTTGAGTGAGCTGGCCGTGGTGCAGGGTCGGGGATTGCAGACTGACTGACGTTTAGCggggg                                                                                                                            <  2:296016/134‑1 (MQ=255)
ggTGGATCTGGAGCAGGGGATGTGGGAAAGCTTGCTTTCCCTCGCGGCCTGGCTGGCGCATGAGACGCAGGTGTTGAGTGAGCTGGCCGTGCTGCAGGGTCAGGTATTGCAGACTGACTGACGTTTAGCggggg                                                                                                                            <  2:627391/134‑1 (MQ=255)
ggTGGATCTGGAGCAGGGGATGTGGGAAAGCTTGCTTTCCCTCGCGGCCTGGCTGGCGCATGAGACGCAGGTGTTGAGTGAGCTGGCCGTGCTGCAGGGTCAGGTATTGCAGACTGACTGACGTTTAGCggggg                                                                                                                            <  1:729648/134‑1 (MQ=255)
              aGGGGATGTGGGAAAGCTTGCTTTCCCTCGCGGCCTGGCTGGCGCATGAGACGCAGGTGTTGAGTGAGCTGGCCGTGCTGCAGGGTCAGGTATTGCAGACTGACTGACGTTTAGCGGGGGGGGGCAGCCGTGc                                                                                                              >  1:19162/1‑133 (MQ=255)
                    tgtgGGAAAGCTTGCTTTCCCTCGCGGCCTGGCTGGCGCATGAGACGCAGGTGTTGAGTGAGCTGGCCGTGCTGCAGGGTCAGGTATTGCAGACTGACTGACGTTTAGCGGGGGGGGGCAGCCGTGCAATCGGc                                                                                                       >  1:823529/1‑134 (MQ=255)
                              cttgcttTCCCTCGCGGCCTGGCTGGCGCATGAGACGCAGGTGTTGAGTGAGCTGGCCGTGCTGCAGGGTCAGGTATTGCAGACTGACTGACGTTTAGCGGGGGGGGGCAGCCGTGCAATCGGCTGcccccccc                                                                                             <  2:892154/134‑7 (MQ=255)
                                                                 cGCAGGTGTTGAGTGAGCTGGCCGTGCTGCAGGGTCAGGTATTGCAGACTGACTGACGTTTAGCGGGGGGGGGCAGCCGTGCAATCGGCTGCCGCGCGCCTCCTCGGCGCGCAGGTCAACTTCCCCGATTTCa                                                           >  1:35620/1‑133 (MQ=255)
                                                                  gCAGGTGTTGAGTGAGCTGGCCGTGCTGCAGGGTCAGGTATTGCAGACTGACTGACGTTTAGCGGGGGGGGGCAGCCGTGCAATCGGCTGCCGCGCGCCTCCTCGGCGCGCAGGTCAACTTCCCCGATTTCaaa                                                         <  1:1078588/134‑1 (MQ=255)
                                                                  gCAGGTGTTGAGTGAGCTGGCCGTGCTGCAGGGTCAGGTATTGCAGACTGACTGACGTTTAGCGGGGGGGGGCAGCCGTGCAATCGGCTGCCGCGCGCCTCCTCGGCGCGCAGGTCAACTTCCCCGATTTCaaa                                                         <  2:35620/134‑1 (MQ=255)
                                                                             gtgagCTGGCCGTGCTGCAGGGTCAGGTATTGCAGACTGACTGACGTTTAGCGGGGGGGGGCAGCCGTGCAATCGGCTGCCGCGCGCCTCCTCGGCGCGCAGGTCAACTTCCCCGATTTCAAAGCGGTTGAAc                                               >  1:631079/1‑133 (MQ=255)
                                                                                  cTGGCCGTGCTGCAGGGTCAGGTATTGCATACTGACTGACGTTTAGCGGGGGGGGGCAGCCGTGCAATCGGCTGCcgcgc                                                                                               >  1:161762/1‑80 (MQ=255)
                                                                                  cTGGCCGTGCTGCAGGGTCAGGTATTGCATACTGACTGACGTTTAGCGGGGGGGGGCAGCCGTGCAATCGGCTGCcgcgc                                                                                               <  2:161762/80‑1 (MQ=255)
                                                                                    ggCCGTGCTGCAGGGTCAGGTATTGCAGACTGACTGACGTTTAGCGGGGGGGGGCAGCCGTGCAATCGGCTGCCGCGCGCctcctc                                                                                       >  1:361161/1‑86 (MQ=255)
                                                                                    ggCCGTGCTGCAGGGTCAGGTATTGCAGACTGACTGACGTTTAGCGGGGGGGGGCAGCCGTGCAATCGGCTGCCGCGCGCctcctc                                                                                       <  2:361161/86‑1 (MQ=255)
                                                                                    ggCCGTGCTGCAGGGTCAGGTATTGCAGACTGACTGACGTTTAGCGGGGGGGGGCAGCCGTGCAATCGGCTGCCGCGCGCCTCCTCGGCGCGCAGGTCAACTTCCCCGATTTCAAAGCGGTTGAACCACACgg                                        >  2:211407/1‑133 (MQ=255)
                                                                                     gCCGTGCTGCAGGGTCAGGTATTGCAGACTGACTGACGTTTAGCGGGGGGGGGCAGCCGTGCAATCGGCTGCCGCGCGCCTCCTCGGCGCGCAGGTCAACTTCCCCGATTTCAAAGCGGTTGAACCACACGGc                                       >  2:879255/1‑133 (MQ=255)
                                                                                     gCCGTGCTGCAGGGTCAGGTATTGCAGACTGACTGACGTTTAGCGGGGGGGGGCAGCCGTGCAATCGGCTGCCGCGCGCCTCCTCGGCGCGCAGGTCAACTTCCCCGATTTCAAAGCGGTTGAACCACACGGCa                                      <  1:741256/134‑1 (MQ=255)
                                                                                           ctgcAGGGTCAGGTATTGCAGACTGACTGACGTTTAGCGGGGGGGGGCAGCCGTGCAATCGGCTGCCGCGCGCCTCCTCGGCGCGCAGGTCAACTTCCCCGATTTCAAAGCGGTTGAACCACACGGCAAAGAc                                 <  1:879255/133‑1 (MQ=255)
                                                                                            tgcAGGGTCAGGTATTGCAGACTGACTGACGTTTAGCGGGGGGGGGCAGCCGTGCAATCGGCTGCCGCGCGCCTCCTCGGCGCGCAGGTCAACTTCCCCGATTTCAAAGCGGTTGAACCACACGGCAAAGACGc                               >  1:101567/1‑134 (MQ=255)
                                                                                                         aTTGCAGACTGACTGACGTTTAGCGGGGGGGGGCAGCCGTGCAATCGGCt                                                                                                      >  1:768623/1‑50 (MQ=255)
                                                                                                         aTTGCAGACTGACTGACGTTTAGCGGGGGGGGGCAGCCGTGCAATCGGCt                                                                                                      <  2:768623/50‑1 (MQ=255)
                                                                                                          ttGCAGACTGACTGACGTTTAGCGGGGGGGGGCAGCCGTGCAATCGGCTGCCGCGCGCCTCCTCGGCGCGCAGGTCAACTTCCCCGATTTCAAAGCGGTTGAACCACACGGCAAAGACGCCATCACGCCGAGTc                 <  2:756382/134‑1 (MQ=255)
                                                                                                            gCAGACTGACTGACGTTTAGCGGGGGGGGGCAGCCGTGCAATCGGCTGCCGCGCGCCTCCTCGGCGCGCAGGTCAACTTCCCCGATTTCAAAGCGGTTGAACCACACGGCAAAGACGCCATCACGCCGAGTCg                <  1:211407/133‑1 (MQ=255)
                                                                                                               gactgactGACGTTTAGCGGGGGGGGGCAGCCGTGCAATCGGCTGCCGCGCGCCTCCTCGGCGCGCAGGTCAACTTCCCCGATTTCAAAGCGGTTGAACCACACGGCaa                                     <  1:864804/109‑1 (MQ=255)
                                                                                                               gactgactGACGTTTAGCGGGGGGGGGCAGCCGTGCAATCGGCTGCCGCGCGCCTCCTCGGCGCGCAGGTCAACTTCCCCGATTTCAAAGCGGTTGAACCACACGGCaa                                     >  2:864804/1‑109 (MQ=255)
                                                                                                                 ctgactgaCGTTTAGCGGGGGGGGGGCAGCCGTGCAATCGGCTGCCGCGCGCCTCCTCGGCGCGCAGGTCAACTTCCCCGATTTCAAAGCGGTTGAACCACACGGCa                                      <  2:1080622/107‑1 (MQ=255)
                                                                                                                 ctgactgaCGTTTAGCGGGGGGGGGGCAGCCGTGCAATCGGCTGCCGCGCGCCTCCTCGGCGCGCAGGTCAACTTCCCCGATTTCAAAGCGGTTGAACCACACGGCa                                      >  1:1080622/1‑107 (MQ=255)
                                                                                                                    actgacGTTTAGCGGGGGGGGGCAGCCGTGCAATCGGCt                                                                                                      <  2:859725/39‑1 (MQ=255)
                                                                                                                    actgacGTTTAGCGGGGGGGGGCAGCCGTGCAATCGGCt                                                                                                      >  1:859725/1‑39 (MQ=255)
                                                                                                                      tgacGTTTAGCGGGGGGGGGCAGCCGTGCAATCGGCTGCCGCGCGCCTCCTCGGCGCGCAGGTCAACTTCCCCGATTTCAAAGCGGTTGAACCACACGGCAAAGACGCCATCACGCCGAGTCGGGCGCAGGc       >  2:648900/1‑132 (MQ=255)
                                                                                                                        acGTTTAGCGGGGGGGGGCAGCCGTGCAATCGGCTGCCGCGCGCCTCCTCGGCGCGCAGGTCAACTTCCCCGATTTCAAAGCGGTTGAACCACACGGCAAAGACGCCATCACGCCGAGTCGGGCGCAGGCgaga   >  2:655674/1‑134 (MQ=255)
                                                                                                                        acGTTTAGCGGGGGGGGGCAGCCGTGCAATCGGCTGCCGCGCGCCTCCTCGGCGCGCAGGTCAACTTCCCCGATTTCAAAGCGGTTGAACCACACGGCAAAGACGCCATCACGCCGAGTCGGGCGCAGGCgaga   >  2:49516/1‑134 (MQ=255)
                                                                                                                         cGTTTAGCGGGGGGGGGCAGCCGTGCAATCGGCTGCCGCGCGCCTCCTCGGCGCGCAGGTCAACTTCCCCGATTTCAAAGCGGTTGAACCACACGGCAAAGACGCCATCACGCCGAGTCGGGCGCAGGCGAGAc  <  2:991647/134‑1 (MQ=255)
                                                                                                                                 |                                                                                                                              
GGTGGATCTGGAGCAGGGGATGTGGGAAAGCTTGCTTTCCCTCGCGGCCTGGCTGGCGCATGAGACGCAGGTGTTGAGTGAGCTGGCCGTGCTGCAGGGTCAGGTATTGCAGACTGACTGACGTTTAGCGGGGGGGGGGCAGCCGTGCAATCGGCTGCCGCGCGCCTCCTCGGCGCGCAGGTCAACTTCCCCGATTTCAAAGCGGTTGAACCACACGGCAAAGACGCCATCACGCCGAGTCGGGCGCAGGCGAGAC  >  AM260480/274332‑274587

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