Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,229,010 C→A 9.3% G400C (GGC→TGC)  PP_3707 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,229,0100CA 9.3% 113.8 / 2.3 40G400C (GGC→TGC) PP_3707hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (17/19);  new base A (2/2);  total (19/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

TACATTCTCAGCGGCATCCGGCGTACCTTCGGCCTCAAAAAAATACGGCCCCCATCTCCCGTCTACCAAATCGCCGTGATACTCAGAGATGGTATC‑TCCACCAAGCATGGCGTCCATATAATAAAGATTGTTGTTCTTCATGCCGCCAAACCTCAAATGAGCAAGCAGCGACTGCATAAAAAATGCGACCTGGCTCTTCAGGTTTGTATGGCTCCTCCTAGAGTACAACCCAGGAATCATCATCCCCTTTAAGCTTTTCTCATACAACGTGTATACATAGTGCTTTTC  >  NC_002947/4228864‑4229151
                                                                                                                                                   |                                                                                                                                             
tACATTCTCAGCGGCATCCGGCGTACCTTCGGCCTCAAAAAAATACGGCCCCCATCTCCCGTCTACCAAATCGCCGTGATACACAGAGATGTTATC‑TCCACCAAGCATGGCGTCCATATAATAAAGATTGTTGTTCTTCATGCCGCCa                                                                                                                                              >  7:216538/1‑148 (MQ=255)
         aGCGGCATCCGGCGTACCTTCGGCCTCAAAAAAATACGGCCCCCATCTCCCGTCTACCAAATCGCCGTGATACTCAGAGATGGTATC‑TCCACCAAGCATGGCGTCCATATAATAAAGATTGTTGTTCTTCATGCCGCCAAACCTCAAAt                                                                                                                                    <  4:127819/149‑1 (MQ=255)
            ggCATCCGGCGTACCTTCGGCCTCAAAAAAATACGGCCCCCATCACCCGTCTACCAAATCGCCGTGATACTCAGAGATGGTATC‑TCCACCAAGCATGGCGTCCATATAATAAAGATTGTTGTTCTTCATGCCGCCAAACCTCAAATGag                                                                                                                                 <  3:244808/149‑1 (MQ=255)
                tCCGGCGTACCTTCGGCCTCAAAAAAATACGGCCCCCATCTCCCGTCTACCAAATCGCCGTGATACTCAGAGATGGTATC‑TCCACCAAGCATGGCGTCCATATAATAAAGATTGTTGTTCTTCATGCCGCCAAACCTCAAATGagcaag                                                                                                                             <  4:6309/149‑1 (MQ=255)
                                    aaaaaaaTACGGCCCCCATCTCCCGTCTACCAAATCGCCGTGATACTCAGAGATGGTATC‑TCCACCAAGCATGGCGTCCATATAATAAAGATTGTTGTTCTTCATGCCGCAAAACCTCAAATGag                                                                                                                                 >  6:151721/1‑125 (MQ=255)
                                    aaaaaaaTACGGCCCCCATCTCCCGTCTACCAAATCGCCGTGATACTCAGAGATGGTATC‑TCCACCAAGCATGGAGTCCATATAATAAAGATTGTTGTTCTTCATGCCGCAAAACCTCAAATGag                                                                                                                                 <  5:151721/125‑1 (MQ=255)
                                              ggCCCCCATCTCCCGTCTACCAAATCGCCGTGATACTCAGAGATGGTATC‑TCCACCAAGCATGGCGTCCATATAATAAAGATTGTTGTTCTTCATGCCGCCAAACCTCAAATGAGCAAGCAGCGACTGCATAAAAAATGCGACATGGct                                                                                               >  3:272987/1‑149 (MQ=255)
                                                          ccGTCTACCAAATCGCCGTGATACTCAGAGATGGTATC‑TCCACCAAGCATGGCGTCCATATAATAAAGATTGTTGTTCTTCATGCCGCCAAACCTCAAATGAGCAAGCAGCGACTGCATAAAAAATGCGACCTGGCTCTTc                                                                                           >  4:101748/1‑141 (MQ=255)
                                                          ccGTCTACCAAATCGCCGTGATACTCAGAGATGGTATC‑TCCACCAAGCATGGCGTCCATATAATAAAGATTGTTGTTCTTCATGCCGCCAAACCTCAAATGAGCAAGCAGCGACTGCATAAAAAATGCGACCTGGCTCTTc                                                                                           <  3:101748/141‑1 (MQ=255)
                                                          ccGTCTACCAAATCGCCGTGATACTCAGAGATGGTATC‑TCCACCAAGCATGGCGTCCATATAATAAAGATTGTTGTTCTTCATGCCGCCAAACCTCAAATGAGCAAGCAGCGACTGCATAAAAAATGCGACCTGGCTCTTCAGGt                                                                                       <  6:279137/145‑1 (MQ=255)
                                                          ccGTCTACCAAATCGCCGTGATACTCAGAGATGGTATC‑TCCACCAAGCATGGCGTCCATATAATAAAGATTGTTGTTCTTCATGCCGCCAAACCTCAAATGAGCAAGCAGCGACTGCATAAAAAATGCGACCTGGCTCTTCAGGTTTGt                                                                                   <  2:131471/149‑1 (MQ=255)
                                                           cGTCTACCAAATCGCCGTGATACTCAGAGATGGTATC‑TCCACCAAGCATGGCGTCCATATAATAAAGATTGTTGTTCTTCATGCCGCCAAACCTCAAATGAGCAAGCAGCGACTGCATAAAAAATGCGACCTGGCTCTTCAGGTTTGTa                                                                                  >  3:219359/1‑149 (MQ=255)
                                                                 ccAAATCGCCGTGATACTCAGAGATGGTATCGTCCACCAAGCATGGCGTCCATATAATAAAGATTGTTGTTCTTCATGCCGCCAAACCTCAAATGAGCAAGCAGCGACTGCATAAAAAATGCGACCTGGCTCTTCAGGTTTGt                                                                                   >  1:33469/1‑143 (MQ=255)
                                                                 ccAAATCGCCGTGATACTCAGAGATGGTATC‑TCCACCAAGCATGGCGTCCATATAATAAAGATTGTTGTTCTTCATGCCGCCAAACCTCAAATGAGCAAGCAGCGACTGCATAAAAAATGCGACCTGGCTCTTCAGGTTTGt                                                                                   <  2:33469/142‑1 (MQ=255)
                                                                   aaaTCGCCGTGATACTCAGAGATGGTATC‑TCCACCAAGCATGGCGTCCATATAATAAAGATTGTTGTTCTTCATGCCGCCAAACCTCAAATGAGCAAGCAGCGACTGCATAAAAAATGCGACCTGGCTCTTCAGGTTTGTATGGctcct                                                                          >  2:84320/1‑149 (MQ=255)
                                                                     aTCGCCGTGATACTCAGAGATGGTATC‑TCCACCAAGCATGGCGTCCATATAATAAAGATTGTTGTTCTTCATGCCGCCAAACCTCAAATGAGCAAGCAGCGACTGCATAAAAAATGCGACCTGGCTCTTCAGGTTTGTATGGCtcctcc                                                                        >  4:69178/1‑149 (MQ=255)
                                                                        gCCGTGATACTCAGAGCTGGTATC‑TCCACCAAGCATGGCGTCCATATAATAAAGATTGTTGTTCTTCATGCCGCCAAACCTCAAATGAGAAAGCAGCGACTGCATAAAAAATGCGACCTGGCTCTTCAGGTTTGTATGGCTACTCCTa                                                                      <  4:272987/148‑1 (MQ=255)
                                                                        gCCGTGATACTCAGAGATTGTATC‑TCCACCAACCATGGCTTCCATATAATAAAGATTGTTGTTCTTCTTGCCGCCAAACCTCAAATGAGCTAGCAGCGACTGCATAAAAAATGCGACCTGGCTCTTCAGGTTTGTATGGCTCCTCCTa                                                                      <  1:218082/148‑1 (MQ=255)
                                                                           gagctACTCAGAGATGGTATC‑TCCACCAAGCATGGCGTCCATATAATAAAGATTGTTGTTCTTCATGCCGCCAAACCTCAAATGAGCAAGCAGCGACTGCATAAAAAATGCGACCTGGCTCTTCAGGTTTGTATGGCTCCTCCTAGAGt                                                                  <  8:216538/145‑1 (MQ=255)
                                                                              aTACTCAGAGATGGTATC‑TCCACCAAGCATGGCGTCCATATAATAAAGATTGTTGTTCTTCATGCCGCCAAACCTCAAATGAGCAAGCAGCGACTGCATAAAAAATGCGACCTg                                                                                                  >  8:85829/1‑114 (MQ=255)
                                                                              aTACTCAGAGATGGTATC‑TCCACCAAGCATGGCGTCCATATAATAAAGATTGTTGATCTTCATGCCGCCAAACCTCAAATGAGCAAGCAGCGACTGCATAAAAAATGCGACCTg                                                                                                  <  7:85829/114‑1 (MQ=255)
                                                                                  tCAGAGATGGTATC‑TCCACCAAGCATGGCGTCCATATAATAAAGATTGTTGTTCTTCATGCCGCCAAACCTCAAATGAGCAAGCAGCGACTGCATAAAAAATGCGACCTGGCTCTTCAGGTTTGTATGGCTCCTCCTagagatcgga                                                             >  7:5503/1‑141 (MQ=255)
                                                                                  tCAGAGATGGTATC‑TCCACCAAGCATGGCGTCCATATAATAAAGATTGTTGTTCTTCATGCCGCCAAACCTCAAATGAGCAAGCAGCGACTGCATAAAAAATGCGACCTGGCTCTTCAGGTTTGTATGGCTCCTCCTa                                                                      <  8:5503/138‑1 (MQ=255)
                                                                                  tCAGAGATGGTATC‑TCCACCAAGCATGGCGTCCATATAATAAAGATTGTTGTTCTTCATGCCGCCAAACCTCAAATGAGCAAGCAGCGACTGCATAAAAAATGCGACCTGGCTCTTCAGGTTTGTATGGCTCCTCCTAGAGTACAAcc                                                            >  6:112549/1‑148 (MQ=255)
                                                                                            tATC‑TCCACCAAGCATGGCGTCCATATAATAAAGATTGTTGTTCTTCATGCCGCCAAACCTCAAATGAGCAAGCAGCGACTGCATAAAAAATGCGACCTGGCTCTTCAGGTTTGTATGGCTCCTCCTAGAGTACAACCCAGGAatcat                                                  <  4:219359/148‑1 (MQ=255)
                                                                                                  acaccaAGCCTGGCGTTAATATAATAAACATTGGTTGTCTTTAATCCGCCAAACCTCAAATGAGCAAGCAGCGACTGCATAAAAAATGCGACCTGGCTCTTCAGGTTTGTATGGCTCCTCCTAGAGTACAACCCAGGAATCATCATcc                                             <  8:51985/147‑1 (MQ=255)
                                                                                                   caccaAGCATGGCGTCCATATAATAAAGATTGTTGTTCTTCATGCCGCCAAACCTCAAATGAGCAAGCAGCGACTGCATAAAAAATGCGACCTGGCTCTTCAGGTTTGTATGGCTCCTCCTAGAGTACAACCCAGGAATCATCATcccc                                           <  1:84320/149‑1 (MQ=255)
                                                                                                      caAGCATGGCGTCCATATAATAAAGATTGTTGTTCTTCATGCCGCCAAACCTCAAATGAGCAAGCAGCGACTGCATAAAAAATGCGACCTGGCTCTTCAGGTTTGTATGGCTcctcct                                                                       >  6:306380/1‑118 (MQ=255)
                                                                                                      caAGCATGGCGTCCATATAATAAAGATTGTTGTTCTTCATGCCGCCAAACCTCAAATGAGCAAGCAGCGACTGCATAAAAAATGCGACCTGGCTCTTCAGGTTTGTATGGCTcctcct                                                                       <  5:306380/118‑1 (MQ=255)
                                                                                                                 tCCATATAATAAAGATTGTTGTTCTTCATGCCGCCAAACCTCAAATGAGCAAGCAGCGACTGCATAAAAAATGCGACCTGGCTCTTCAGGTATGTATGGCTCCTCCTAGAGTACAACCCAGGAATCATCATCCCCTTTAAGCTTTtct                              >  1:315287/1‑148 (MQ=255)
                                                                                                                        aataaAGATTGTTGTTCTTCATGCCGCCAAACCTCAAATGAGCAAGCAGCGACTGCATAAAAAATGCGACCTGGCTCTTCAGGTTTGTATGGCTCCTCCTAGAGTACAACCCAGGAATCATCATCCCCTTTAAGCTTTTCTCATACAAc                      >  7:216799/1‑149 (MQ=255)
                                                                                                                          taaAGATTGTTGTTCTTCATGCCGCCAAACCTCAAATGAGCAAGCAGCGACTGCATAAAAAATGCGACCTGGCTCTTCAGGTTTGTATGGCTCCTCCTAGAGTACAACCCAGGAATCATCATCCCCTTTAAGCTTTTCTCATACAACgt                    >  7:124643/1‑149 (MQ=255)
                                                                                                                           aaaGATTGTTGTTCTTCATGCCGCCAAACCTCAAATGAGCAAGCAGCGACTGCATAAAAAATGCGACCTGGCTCTTCAGGTTTGTATGGCTCCTCCTAGAGTACAACCCAGGaa                                                      >  1:285485/1‑114 (MQ=255)
                                                                                                                           aaaGATTGTTGTTCTTCATGCCGCCAAACCTCAAATGAGCAAGCAGCGACTGCATAAAAAATGCGACCTGGCTCTTCAGGTTTGTATGGCTCCTCCTAGAGTACAACCCAGGaa                                                      <  2:285485/114‑1 (MQ=255)
                                                                                                                              gATTGTTGTTCTTCATGCCGCCAAACCTCAAATGAGCAAGCAGCGACTGCATAAAAAATGCGACCTGGCTCTTCAGGTTTGTATGGCTCCTCCTAGAGTACAACCCAGGAATCATCATCCCCTTTAAGCTTTTCTCATACAACGTGta                 >  2:107199/1‑148 (MQ=255)
                                                                                                                              gATTGTTGTTCTTCATGCCGCAAAACCTCAAATGAGCAAGCAGCGACTGCATAAAAAATGCGACCTGGCTCTTCAGGTTTGTATGGCTCCTCCTAGAGTACAACCCAGGAATCATCATCCCCTTTAAGCTTTTCTCAt                           >  2:28470/1‑138 (MQ=255)
                                                                                                                              gATTGTTGTTCTTCATGCCGCAAAACCTCAAATGAGCAAGCAGCGACTGCATAAAAAATGCGACCTGGCTCTTCAGGTTTGTATGGCTCCTCCTAGAGTACAACCCAGGAATCATCATCCCCTTTAAGCTTTTCTCAt                           <  1:28470/138‑1 (MQ=255)
                                                                                                                                   ttgttCTTCATGCCGCCAAACCTCAAATGAGCAAGCAGCGACTGCATAAAAAATGCGACCTGGCTCTTCAGGTTTGTATGGCTCCTCCTAGAGTACAACCCAGGAATCATCATCCCCTTTAAGCTTTTCTCATACAACGTGTATACAt            <  1:107199/148‑1 (MQ=255)
                                                                                                                                         ttcATGCCGCCAAACCTCAAATGAGCAAGCAGCGACTGCATAAAAAATGCGACCTGGCTCTTCAGGTTTGTATGGCTCCTCCTAGAGTACAACCCAGGAATCATCATCCCCTTTAAGCTTTTCTCATACAACGTGTATACATAGTGCtt     >  3:154750/1‑149 (MQ=255)
                                                                                                                                            aTGCCGCCAAACCTCAAATGAGCAAGCAGCGACTGCATAAAAAATGCGACCTGGCTCTTCAGGTTTGTATGGCTCCTCCTAGAGTACAACCCAGGAATCATCATCCCCTTTAAGCTTTTCTCATACAACGTGTATACATAGTGCTTTTc  <  8:216799/149‑1 (MQ=255)
                                                                                                                                                   cAAACCTCAAATGAGCAAGCAGCGACTGCATAAAAAATGCGACCTGGCTCTTCAGGTTTGTATGGCTCCTCCTAGAGTACAACCCAg                                                         >  3:302822/1‑87 (MQ=255)
                                                                                                                                                   cAAACCTCAAATGAGCAAGCAGCGACTGCATAAAAAATGCGACCTGGCTCTTCAGGTTTGTATGGCTCCTCCTAGAGTACAACCCAg                                                         <  4:302822/87‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
TACATTCTCAGCGGCATCCGGCGTACCTTCGGCCTCAAAAAAATACGGCCCCCATCTCCCGTCTACCAAATCGCCGTGATACTCAGAGATGGTATC‑TCCACCAAGCATGGCGTCCATATAATAAAGATTGTTGTTCTTCATGCCGCCAAACCTCAAATGAGCAAGCAGCGACTGCATAAAAAATGCGACCTGGCTCTTCAGGTTTGTATGGCTCCTCCTAGAGTACAACCCAGGAATCATCATCCCCTTTAAGCTTTTCTCATACAACGTGTATACATAGTGCTTTTC  >  NC_002947/4228864‑4229151

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