Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 26,844 A→T 23.3% intergenic (+40/+145) PP_0021 → / ← PP_0022 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294726,8440AT23.3% 86.7 / 26.0 43intergenic (+40/+145)PP_0021/PP_0022hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (14/19);  new base T (6/4);  total (20/23)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.31e-01

TTCGAATCCCTTATCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCACGCCGCGGTGTCCTA  >  NC_002947/26697‑26991
                                                                                                                                                    |                                                                                                                                                   
ttCGAATCCCTTATCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGt                                                                                                                                                   >  2:1032575/1‑150 (MQ=255)
   gAATCCCTTATCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGt                                                                                                                                                <  1:263595/150‑1 (MQ=255)
      tCCCTTATCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTCACCACAGTAGTTTa                                                                                                                                             >  2:1030797/1‑150 (MQ=255)
        ccTTATCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACa                                                                                                                                           >  1:824609/1‑150 (MQ=255)
        ccTTATCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACa                                                                                                                                           <  2:1084701/150‑1 (MQ=255)
         cTTATCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAg                                                                                                                                          >  2:1134650/1‑150 (MQ=255)
                     ccGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTg                                                                                                                              <  1:963906/150‑1 (MQ=255)
                          tGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTa                                                                                                                         <  2:1194867/150‑1 (MQ=255)
                             tcatcCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACaa                                                                                                                      <  1:624301/150‑1 (MQ=255)
                              catcCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAAATGACCGTTTGGTTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTCAACCAGCATGATGCACGGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAAc                                                                                                                     <  1:822238/150‑1 (MQ=255)
                                     aTCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTa                                                                                                              <  1:943220/150‑1 (MQ=255)
                                                       tGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGt                                                                                                                               <  2:488702/115‑1 (MQ=255)
                                                       tGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGt                                                                                                                               >  1:488702/1‑115 (MQ=255)
                                                         aaGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACt                                                                                          <  1:1073495/150‑1 (MQ=255)
                                                          aGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTc                                                                                         <  1:874023/150‑1 (MQ=255)
                                                           gATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGa                                                                                                           <  2:77931/131‑1 (MQ=255)
                                                           gATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGa                                                                                                           >  1:77931/1‑131 (MQ=255)
                                                                 ccGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGc                                                                                  >  2:1080321/1‑150 (MQ=255)
                                                                        gaACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTgg                                                                                               <  2:793244/130‑1 (MQ=255)
                                                                        gaACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTgg                                                                                               >  1:793244/1‑130 (MQ=255)
                                                                            tGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACAACAATGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCCCGctac                                                                       >  2:554530/1‑148 (MQ=255)
                                                                                   tCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACg                                                                                                                    >  1:317562/1‑98 (MQ=255)
                                                                                   tCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACg                                                                                                                    <  2:317562/98‑1 (MQ=255)
                                                                                     tgctgGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGAtt                                                              <  1:1134650/150‑1 (MQ=255)
                                                                                      gctgGGGGTGGCTCAGGCCTGAGGAGGCTTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGTCAGCGCGCTCCTCATTGATTc                                                             <  2:880444/150‑1 (MQ=38)
                                                                                      gctgGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTc                                                             <  1:1032575/150‑1 (MQ=255)
                                                                                         ggAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGAGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGtt                                                          <  2:857513/150‑1 (MQ=255)
                                                                                                           aggaggCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAacac                                        <  2:824609/150‑1 (MQ=255)
                                                                                                                        tAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCgg                           >  2:842593/1‑150 (MQ=255)
                                                                                                                             cACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTgg                                                                                               <  1:816047/77‑1 (MQ=255)
                                                                                                                             cACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTgg                                                                                               >  2:816047/1‑77 (MQ=255)
                                                                                                                                     ttCACTGTAAACCACTGTGGTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCAcgc              >  2:799078/1‑150 (MQ=38)
                                                                                                                                        aCTGTAAACCACTGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGctcc                                                                       <  2:115482/91‑1 (MQ=37)
                                                                                                                                        aCTGTAAACCACTGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGctcc                                                                       >  1:115482/1‑91 (MQ=37)
                                                                                                                                        aCTGTAAACCACTGTGGTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCAcgccgc           >  1:597448/1‑150 (MQ=38)
                                                                                                                                         ctgtaaaCCACTGTGGTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTc                                                         <  1:1028823/98‑1 (MQ=25)
                                                                                                                                         ctgtaaaCCACTGTGGTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTc                                                         >  2:1028823/6‑103 (MQ=25)
                                                                                                                                          tGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCACGCCGCgg         >  1:546575/1‑150 (MQ=255)
                                                                                                                                            taaaCCACTGTGGTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCa                                                                    <  2:594442/87‑1 (MQ=25)
                                                                                                                                            taaaCCACTGTGGTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCa                                                                    >  1:594442/3‑89 (MQ=25)
                                                                                                                                                ccACTGTGGTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCACGCCGCGgtgt      <  2:838502/147‑1 (MQ=255)
                                                                                                                                                ccACTGTGGTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCACGCCGCGgtgt      >  1:838502/1‑147 (MQ=255)
                                                                                                                                                ccACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCACGCCGCGGTGTcc    >  1:211355/1‑150 (MQ=255)
                                                                                                                                                 cacaGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCACGCCGCGGTGTCCt   >  1:177914/1‑150 (MQ=255)
                                                                                                                                                  acaGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGAGCTCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCACGCCGCGGTGTCCta  >  1:210482/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TTCGAATCCCTTATCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCACGCCGCGGTGTCCTA  >  NC_002947/26697‑26991

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: