Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,470,272 A→C 9.3% V321G (GTG→GGG)  PP_3961 ← transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,470,2720AC 9.3% 111.8 / 4.7 43V321G (GTG→GGG) PP_3961transporter
Reads supporting (aligned to +/- strand):  ref base A (19/20);  new base C (2/2);  total (21/22)
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.88e-01

TGAATCCGGCCTGGGCACCCCACACATCGACCACCTGCCCAGATGCTGCGGCACCGCCGGCAACTCCAACGTTGAGACCGGCGAGCAACCAGGTCATCCCTTCGGTCAGCTTGCTTTCGGGCACGATGCGCTCTACCAGCGACATCGCCACGATCATGGTCGGGGCGAAGAACAGGCCCGCAACCAGCACGGCGATGGCCAGCCCGATGATGTTTCCCACCACCAGCAGCGGCAAGGTTGTCGCCGCCGTCGCCAACCCACCCAGCAGTAGCAGCTTGTGTAGCGGGGTATC  >  NC_002947/4470123‑4470414
                                                                                                                                                     |                                                                                                                                              
tGAATCCGGCCTGGGCACCCCACACATCGACCACCTGCCCAGATGCTGCGGCACCGCCGGCAACTCCAACGTTGAGACCGGCGAGCAACCAGGTCATCCCTTCGGTCAGCTTGCTTTCGGGCACGATGCGCTCTACCAGCGACATCGCCa                                                                                                                                                <  1:523220/150‑1 (MQ=255)
        gCCTGGGCACCCCACACATCGACCACCTGCCCAGATGCTGCGGCACCGCCGGCAACTCCAACGTTGAGACCGGCGAGCAACCAGGTCATCCCTTCGGTCAGCTTGCTTTCGGGCACGATGCGCTCTACCAGCGACATCGCCACGATCATg                                                                                                                                        >  1:295149/1‑150 (MQ=255)
                  cccACACATCGACCACCTGCCCAGATGCTGCGGCACCGCCGGCAACTCCAACGTTGAGACCGGCGAGCAACCAGGTCATCCCTTCGGTCAGCTTGCTTTCGGGCACGATGCGCTCTACCAGCGACATCGCCACGATCATGGTCGGGGCga                                                                                                                              >  2:569273/1‑150 (MQ=255)
                   ccACACATCGACCACCTGCCCAGATGCTGCGGCACCGCCGGCAACTCCAACGTTGAGACCGGCGAGCAACCAGGTCATCCCTTCGGTCAGCTTGCTTTCGGGCACGATGCGCTCTACCAGCGACATCGCCCCGACCATGATCGTGGCgat                                                                                                                             <  2:287414/150‑2 (MQ=255)
                   ccACACATCGACCACCTGCCCAGATGCTGCGGCACCGCCGGCAACTCCAACGTTGAGACCGGCGAGCAACCAGGTCATCCCTTCGGTCAGCTTGCTTTCGGGCACGATGCGCTCTACCAGCGACATCGCCACGATCATGGTCGGGGCgaa                                                                                                                             <  2:48526/150‑1 (MQ=255)
                          tCGACCACCTGCCCAGATGCTGCGGCACCGCCGGCAACTCCAACGTTGAGACCGGCGAGCAACCAGGTCATCCCTTCGGTCAGCTTGCTTTCGGGCACGATGCGCTCTACCAGCGACATCGCCACGATCATGGTCGGGGCGAAGAACAgg                                                                                                                      >  2:166060/1‑150 (MQ=255)
                                 ccTCCCCAGTTGCTGCGGCACCGCCGGCAACTCCAACGTTGAGACCGGCGAGCAACCAGGTCATCCCTTCGGTCAGCTTGCTTTCGGGCACGATGCGCTCTACCAGCGACATCGCCACGATCATGGTCGGGGCGAAGAACAGGCCCGCaa                                                                                                               <  1:407457/150‑1 (MQ=255)
                                              tgcGGCACCGCCGGCAACTCCAACGTTGAGACCGGCGAGCAACCCGGACATCCCTATGGACAGTTTGCGGTCTGGCGCGATGCGCTCTACCAGAGGCATGGCCACTAAGATGGTCGGGGCGAAGAACAGGCCCGCAACCAGCACGGCGAt                                                                                                  <  1:390660/150‑1 (MQ=255)
                                               gcGGCACCGCCGGCAACTCCAACGTTGAGACCGGCGAGCAACCAGGTCATCCCTTCGGTCAGCTTGCTTTCGGGCACGATGCGCTCTACCAGCGACATCGCCACGATCATGGTCGGGGCGAAGAACAGGCCCGCAACCAGCACGGCGATg                                                                                                 <  1:497486/150‑1 (MQ=255)
                                               gcGGCACCGCCGGCAACTCCAACGTTGAGACCGGCGAGCAACCAGGTCATCCCTTCGGTCAGCTTGCTTTCGGGCACGATGCGCTCTACCAGCGACATCGCCACGATCATGGTCGGGGCGAAGAACAGGCCCGCAACCAGCACGGCGATg                                                                                                 <  1:45268/150‑1 (MQ=255)
                                               gcGGCACCGCCGGCAACTCCAACGTTGAGACCGGCGAGCAACCAGGTCATCCCTTCGGTCAGCTTGCTTTCGGGCACGATGCGCTCTACCAGCGACATCGCCACGATCATGGTCGGGGCGAAGAACAGGCCCGCAACCAGCACGGCGATg                                                                                                 <  1:45264/150‑1 (MQ=255)
                                                                tCCAACGTTGAGACCGGCGAGCAACCAGGTCATCCCTTCGGTCAGCTTGCTTTCGGGCACGATGCGCTCTACCAGCGACATCGCCACGATCATGGTCGGGGCGAAGAACAGGCCCGCAAc                                                                                                              <  1:207001/120‑1 (MQ=255)
                                                                tCCAACGTTGAGACCGGCGAGCAACCAGGTCATCCCTTCGGTCAGCTTGCTTTCGGGCACGATGCGCTCTACCAGCGACATCGCCACGATCATGGTCGGGGCGAAGAACAGGCCCGCAAc                                                                                                              >  2:207001/1‑120 (MQ=255)
                                                                           gaCCGGCGAGCAACCAGGTCATCCCTTCGGTCAGCTTGCTTTCGGGCACGATGCGCTCTACCAGCGACATCGCCACGATCATGGTCGGGGCGAAGAACAGGCCCGCAACCAGCACGGCGATGGCCAGCCCGATGATGTTTCccaccacca                                                                     <  2:295149/150‑1 (MQ=255)
                                                                                        ccAGGTCATCCCTTCGGTCAGCTTGCTTTCGGGCACGATGCGCTCTACCAGCGACATCGCCACGATCATGGTCGGGGCGAAGAACAGGCCCGCAAc                                                                                                              <  2:407767/96‑1 (MQ=255)
                                                                                        ccAGGTCATCCCTTCGGTCAGCTTGCTTTCGGGCACGATGCGCTCTACCAGCGACATCGCCACGATCATGGTCGGGGCGAAGAACAGGCCCGCAAc                                                                                                              >  1:407767/1‑96 (MQ=255)
                                                                                        ccAGGTCATCCCTTCGGTCAGCTTGCTTTCGGGCACGATGCGCTCTACCAGCGACATCGCCACGATCATGGGCGGGGCGAAGAACAGGCCCGCAACCAGCACGGCGATGGCCAGCCCGATGATGTTTCCCACCACCAGCAGCGGCAAGGt                                                        >  1:703976/1‑150 (MQ=255)
                                                                                           ggTCATCCCTTCGGTCAGCTTGCTTTCGGGCACGATGCGCTCTACCAGCGACATCGCCACGATCATGGTCGGGGCGAAGAACAGGCCCGCAACCAGCACGGCGATGGCCAGCCCAATGATGTTTCCCACCACCAGCAGCGGCAAGGTTGt                                                     >  1:271019/1‑150 (MQ=255)
                                                                                               aTCCCTTCGGTAAGCTTGCTTTCGGGCACGATGCGCTCTACCAGCGACATCGCCACGATCATGGTCGGGGCGAAGCGGATTACTGCGATCAGCACGGCGATGGACAGCCCGATGATGTTTGTAACGAAGAGGAGCGGGAAGGTTGTcgcc                                                 >  1:317624/1‑150 (MQ=255)
                                                                                                   cTTCGGTCAGCTTGCTTTCGGGCACGATGCGCTCTACCAGCGACATCGCCACGATCATGGTCGGGGCGAAGAACAGGCCCGCAACCAGCACGGCGATGGCCa                                                                                             <  1:489712/102‑1 (MQ=255)
                                                                                                   cTTCGGTCAGCTTGCTTTCGGGCACGATGCGCTCTACCAGCGACATCGCCACGATCATGGTCGGGGCGAAGAACAGGCCCGCAACCAGCACGGCGATGGCCa                                                                                             >  2:489712/1‑102 (MQ=255)
                                                                                                                gcttTCGGGCACGATGCGCTCTACCAGCGACATCGCCACGATCATGGTCGGGGCGAAGAACAGGCCCGCAACCAGCACGGCGATGGCCAGCCCGATGATGTTTCCCACCACCAGCAGCGGCAAGGTTGTCGCCGCCGTCGCCAacccacc                                >  2:473690/1‑150 (MQ=255)
                                                                                                                 cttTCGGGCACGATGCGCTCTACCAGCGACATCGCCACGATCATGGTCGGGGCGAAGAACAGGCCCGCAACCAGCACGGCGATGGCCAGCCCGATGATGTTTCCCACCACCAGCAGCGGCAAGGTTGTCGCCGCCGTCGCCAacccac                                 >  1:400413/1‑148 (MQ=255)
                                                                                                                 cttTCGGGCACGATGCGCTCTACCAGCGACATCGCCACGATCATGGTCGGGGCGAAGAACAGGCCCGCAACCAGCACGGCGATGGCCAGCCCGATGATGTTTCCCACCACCAGCAGCGGCAAGGTTGTCGCCGCCGTCGCCAacccac                                 <  2:400413/148‑1 (MQ=255)
                                                                                                                  tttCGGGCACGATGCGCTCTACCAGCGACATCGCCACGATCATGGTCGGGGCGAAGAACAGGCCCGCAAc                                                                                                              <  2:98875/70‑1 (MQ=255)
                                                                                                                  tttCGGGCACGATGCGCTCTACCAGCGACATCGCCACGATCATGGTCGGGGCGAAGAACAGGCCCGCAAc                                                                                                              >  1:98875/1‑70 (MQ=255)
                                                                                                                       ggCAAGATCCGCTCTACCAGCGACATCGCCACGATCATGGTCGGGGCGAAGAACAGGCCCGCAACCAGCACGGCGATGGCCAGCCCGATGATGTTTCCCACCACCAGCAGCGGCAAGGTTGTCGCCGCCGTCGCCAACCCACCCAGCAGt                         <  1:596045/150‑1 (MQ=255)
                                                                                                                        gCACGATGCGCTCTAACCGAAACATCGCCACACACTCTGTCGTGCCGGAGGATATTCCCGCAACCAGCACGGCGATGGCCAGCCCAATGATGTTTCCCACCACCAGCAGCGGCAAGGTTGTCGCCGCCGTCGCCAACCCACCCAGCAGTa                        <  2:271019/150‑1 (MQ=255)
                                                                                                                           cGATGCGCTCTACCAGCGACATCGCCACGATCATGGTCGGGGCGAAGAACAGGCCCGCAACCAGCACGGCGATGGCCAGCCCGATGATGTTTCCCACCACCAGCAGCGGCAAGGTTGTCGCCGCCGTCGCCAACCCACCCAGCAGTagca                     >  2:224187/1‑150 (MQ=255)
                                                                                                                            gATGCGCTCTACCAGCGACATCGCCACGATCATGGTCGGGGCGAAGAACAGGCCCGCAACCAGCACGGCGATGGCCAGCCCGATGATGTTTCCCACCACcag                                                                    >  2:399446/1‑102 (MQ=255)
                                                                                                                            gATGCGCTCTACCAGCGACATCGCCACGATCATGGTCGGGGCGAAGAACAGGCCCGCAACCAGCACGGCGATGGCCAGCCCGATGATGTTTCCCACCACcag                                                                    <  1:399446/102‑1 (MQ=255)
                                                                                                                            gATGCGCTCTACCAGCGACATCGCCACGATCATGGTCGGGGCGAAGAACAGGCCCGCAACCAGCACGGCGATGGCCAGCCCGATGATGTTTCCCACCACCAGCAGCGGGAAGGTTGTCGCCGCCGTCGCCAACCCACCCAGCAGTagcag                    >  1:303380/1‑150 (MQ=255)
                                                                                                                            gATGCGCTCTACCAGCGACATCGCCACGATCATGGTCGGGGCGAAGAACAGGCCCGCAACCAGCACGGCGATGGCCAGCCCGATGATGTTTCCCACCACCAGCAGCGGCAAGGTTGTCGCCGCCGTCGCCAACCCACCCAGCAGTagcag                    >  2:244690/1‑150 (MQ=255)
                                                                                                                             aTGCGCTCTACCAGCGACATCGCCACGATCATGGTCGGGGCGAAGAACAGGCCCGCAACCAGCACGGCGATGGCCAGCCCGATGATGTTTCCCACCACCAGCAGCGGCAAGGTTGTCGCCGCCGTCGCCAACCCACCCAGCAGTagcagc                   <  2:455419/150‑1 (MQ=255)
                                                                                                                               gcgcTCTACCAGCGACATCGCCACGATCATGGTCGGGGCGAAGAACAGGCCCGCAACCAGCACGGCGATGGCCAGCCCGATGATGTTTCCCACCACCAGCAGCGGCAAGGTTGTCGCCGCCGTCGCCAACCCACCCAGCAGTAGCAGCtt                 >  2:388036/1‑150 (MQ=255)
                                                                                                                               gcgcTCTACCAGCGACATCGCCACGATCATGGTCGGGGCGAAGAACAGGCCCGCAACCAGCACGGCGATGGCCAGCCCGATGATGTTTCCCACCACCAGCAGCGGCAAGGTTGTCGCCGCCGTCGCCAACCCACCCAGCAGTAGCAGCtt                 <  1:166060/150‑1 (MQ=255)
                                                                                                                                  ctctACCAGCGACATCGCCACGATCATGGTCGGGGCGAAGAACAGGCCCGCAAc                                                                                                              <  1:749809/54‑1 (MQ=255)
                                                                                                                                  ctctACCAGCGACATCGCCACGATCATGGTCGGGGCGAAGAACAGGCCCGCAAc                                                                                                              >  2:749809/1‑54 (MQ=255)
                                                                                                                                       ccAGCGACATCGCCACGATCATGGTCGGGGCGAAGAACAGGCCCGCAACCAGCACGGCGATGGCCAGCCCGATGATGTTTCCCACCACCAGCAGCGGCAAGGTTGTCGCCGCCGTCGCCAACCCACCCAGCAGTAGCAGCTTGTGTAGCg         >  1:86818/1‑150 (MQ=255)
                                                                                                                                       ccAGCGACATCGCCACGATCATGGTCGGGGCGAAGAACAGGCCCGCAACCAGCACGGCGATGGCCAGCCCGATGATGTTTCCCACCACCAGCAGCGGCAAGGTTGTCGCCGCCGTCGCCAACCCACCCAGCAGTAGCAGCTTGTGTAGCg         <  1:224187/150‑1 (MQ=255)
                                                                                                                                        cAGCGACATCGCCACGATCATGGTCGGGGCGAAGAACAGGCCCGCAACCAGCACGGCGATGGCCAGCCCGATGATGTTTCCCACCACCAGCAGCGGCAAGGTTGTCGCCGCCGTCGCCAACCCACCCAGCAGTAGCAGCTTGTGTAGCgg        <  1:569273/150‑1 (MQ=255)
                                                                                                                                              cttCGCCCCGACCATGATCGTGGCGAAGAACAGGCCCGCAACCAGCACGGCGATGGCCAGCCCGATGATGTTTCCCACCACCAGCAGCGGCAAGGTTGTCGCCGCCGTCGCCAACCCACCCAGCAGTAGCAGCTTGTGTAGCGGGGTATc  >  1:319423/3‑150 (MQ=255)
                                                                                                                                               ttCGCCCCGACCATGATCGTGGCGAAGAACAGGCCCGCAACCAGCACGGCGATGGCCAGCCCGATGATGTTTCCCACCACCAGCAGCGGCAAg                                                          >  2:44180/2‑93 (MQ=38)
                                                                                                                                               ttCGCCCCGACCATGATCGTGGCGAAGAACAGGCCCGCAACCAGCACGGCGATGGCCAGCCCGATGATGTTTCCCACCACCAGCAGCGGCAAg                                                          <  1:44180/92‑1 (MQ=38)
                                                                                                                                                     |                                                                                                                                              
TGAATCCGGCCTGGGCACCCCACACATCGACCACCTGCCCAGATGCTGCGGCACCGCCGGCAACTCCAACGTTGAGACCGGCGAGCAACCAGGTCATCCCTTCGGTCAGCTTGCTTTCGGGCACGATGCGCTCTACCAGCGACATCGCCACGATCATGGTCGGGGCGAAGAACAGGCCCGCAACCAGCACGGCGATGGCCAGCCCGATGATGTTTCCCACCACCAGCAGCGGCAAGGTTGTCGCCGCCGTCGCCAACCCACCCAGCAGTAGCAGCTTGTGTAGCGGGGTATC  >  NC_002947/4470123‑4470414

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑