Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 956,586 G→T 13.6% intergenic (+16/‑98) alaC → / → PP_0818 aminotransferase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947956,5860GT13.6% 108.2 / 7.5 44intergenic (+16/‑98)alaC/PP_0818aminotransferase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (21/17);  new base T (2/4);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 4.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

GGGGATCGGCTTTGGTGACTATGGGGACGACCATGTGCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCACCGCCTACTTCTGCGGCACTGGCTCGCACCGTTTC  >  NC_002947/956441‑956731
                                                                                                                                                 |                                                                                                                                                 
ggggATCGGCTTTGGTGACTATGGGGACGACCATGTGCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAagc                                                                                                                                               <  1:64297/150‑1 (MQ=255)
ggggATCGGCTTTGGTGACTATGGGGACGACCATGTGCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAagc                                                                                                                                               <  1:642370/150‑1 (MQ=255)
   gATCGGCTTTGGTGACTATGGGGACGACCATGTGCGCTTTGCGTTGATCGAGAACCAGGACGGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAagcagc                                                                                                                                            >  1:147400/1‑150 (MQ=255)
         cTTTGGTGACTATGGGGACGACCATGTGCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCTTCGCAGCCCCatt                                                                                                                                      <  1:637080/150‑3 (MQ=255)
         cTTTGGTGACTAAGGGGACGACCATGTGCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGACCATCCAATGCGGCTGCTTCGCAGCCCCatt                                                                                                                                      <  1:337015/150‑3 (MQ=255)
              gTGACTATGGGGACGACCATGTGCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGcat                                                                                                                                 >  1:476754/1‑148 (MQ=255)
                   tATGGGGACGACCATGTGCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTgg                                                                                                                            <  1:929/150‑1 (MQ=255)
                   tATGGGGACGACCATGTGCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTgg                                                                                                                            >  1:203075/1‑150 (MQ=255)
                       gggACGACAGTGTGCGCGTTGCGTTGATCGAGAAGCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTCCCGGGCTGATGGTTGGGGTCGCAAGCAATCCTAATGGGGCTGCGAAGCAGCCCCAACGGCCTCTTGTGTGa                                                                                                                        <  2:1073976/150‑1 (MQ=255)
                                   tgCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCgcctg                                                                                                            >  1:1091515/1‑150 (MQ=255)
                                          tGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTg                                                                                                     <  2:1091515/150‑1 (MQ=255)
                                             gTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCTTCGCAGCCCCatt                                                                                                                                      <  2:277349/114‑3 (MQ=255)
                                             gTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCTTCGCAGCCCCatt                                                                                                                                      >  1:277349/1‑112 (MQ=255)
                                             gTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGtt                                                                                                  >  1:704766/1‑150 (MQ=255)
                                                 aTCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCTTCGCAGCCCCattgg                                                                                                                                    >  2:114935/1‑108 (MQ=255)
                                                 aTCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCTTCGCAGCCCCattgg                                                                                                                                    <  1:114935/112‑5 (MQ=255)
                                                 aTCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTAccc                                                                                              <  2:203075/150‑1 (MQ=255)
                                                  tCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGa                                                                                                                        <  1:276040/123‑1 (MQ=255)
                                                  tCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGa                                                                                                                        >  2:276040/1‑123 (MQ=255)
                                                         ccAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTg                                                                                                     <  2:1002229/135‑1 (MQ=255)
                                                         ccAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTg                                                                                                     >  1:1002229/1‑135 (MQ=255)
                                                          cAGGAACGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTa                                                                                     >  1:224911/1‑150 (MQ=255)
                                                             gACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCg                                                                                  >  1:1064470/1‑150 (MQ=255)
                                                                  tATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGc                                                                             >  1:987877/1‑150 (MQ=255)
                                                                   aTTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGaa                                                                                                                                                 <  2:831255/81‑1 (MQ=255)
                                                                   aTTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGaa                                                                                                                                                 >  1:831255/1‑81 (MQ=255)
                                                                       gCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCa                                                                        >  2:557333/1‑150 (MQ=255)
                                                                         cAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCGTCTGGTTACCCGTACGCCTACCGTCAGCACTCAtt                                                                      >  1:264157/1‑150 (MQ=255)
                                                                                   gcgGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCaaa                                                            <  2:147400/150‑1 (MQ=255)
                                                                                      ggATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCTCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTAc                                                                                    <  2:793552/123‑1 (MQ=255)
                                                                                      ggATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTAc                                                                                    >  1:793552/1‑123 (MQ=255)
                                                                                         tCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATc                                                      >  1:1082562/1‑150 (MQ=255)
                                                                                             gCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATccccc                                                  >  1:824409/1‑150 (MQ=255)
                                                                                                   gTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAgt                                            <  2:224911/150‑1 (MQ=255)
                                                                                                   gTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAgt                                            <  2:704766/150‑1 (MQ=255)
                                                                                                      ccGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTc                                         <  2:1064470/150‑1 (MQ=255)
                                                                                                      ccGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTc                                         >  1:725086/1‑150 (MQ=255)
                                                                                                        gggCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGc                                                                                        >  1:926317/1‑101 (MQ=255)
                                                                                                        gggCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGc                                                                                        <  2:926317/101‑1 (MQ=255)
                                                                                                          gCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCa                                     >  2:121109/1‑150 (MQ=255)
                                                                                                                 gTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCACCGCCTa                              >  1:64909/1‑150 (MQ=255)
                                                                                                                        tCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCACCGCCTACTTCTGc                       <  1:557333/150‑1 (MQ=255)
                                                                                                                                  gCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCACCGCCTACTTCTGCGGCACTGGCt             >  1:51204/1‑150 (MQ=255)
                                                                                                                                             cTGCGAAGTAGCCCAAACTGCCGCTTGGGTGCGCTGTCAGCCGGAGTTGGGGGGAGCCGTAGGCCGACCGTGAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCACCGCCTACTTCTGCGGCACTGGCTCGCACCGTTTc  <  2:264157/150‑1 (MQ=255)
                                                                                                                                                 gAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAgt                                            >  1:785849/1‑104 (MQ=255)
                                                                                                                                                 gAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAgt                                            <  2:785849/104‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GGGGATCGGCTTTGGTGACTATGGGGACGACCATGTGCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCACCGCCTACTTCTGCGGCACTGGCTCGCACCGTTTC  >  NC_002947/956441‑956731

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