Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,380,027 C→T 7.8% A167T (GCA→ACA)  PP_3856 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,380,0270CT 7.8% 144.5 / 2.3 51A167T (GCA→ACA) PP_3856hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (22/25);  new base T (2/2);  total (24/27)
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.95e-01

GTGACCAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCATGCGGTGCGGCTCAACAGAGGGCCAAGGTAGTCGGCAAACTGGGTGGACAGGCTGGTGGTACCAATGCGGAATTTGCCAGCCGAGGACGCCAAGAACCAGTAGCTTTGCCCAGATCCAGGCAACACCGAACTGCTGGCAGG  >  NC_002947/4379882‑4380165
                                                                                                                                                 |                                                                                                                                          
gTGACCAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCAtgcggt                                                                                                                                         <  2:153387/149‑1 (MQ=255)
     cAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGTATCAAGCGAAGCCAtgcggtgcgg                                                                                                                                     <  7:295922/148‑1 (MQ=255)
     cAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGAAGCCAtgcggtgcgg                                                                                                                                     >  8:295922/1‑148 (MQ=255)
         ccACCCATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCATGCGGTGCGGCTCaa                                                                                                                                >  8:319479/1‑149 (MQ=255)
         ccACCCATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCATGCGGTGCGGCTCaa                                                                                                                                >  3:276918/1‑149 (MQ=255)
             ccATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCATGCGGTGCGGCTCa                                                                                                                                 >  4:206249/1‑144 (MQ=255)
             ccATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCATGCGGTGCGGCTCa                                                                                                                                 <  3:206249/144‑1 (MQ=255)
             ccATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCATGCGGTGCGGCTCAACaga                                                                                                                            <  8:200990/149‑1 (MQ=255)
             ccATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCATGCGGTGCGGCTCAACaga                                                                                                                            <  6:292444/149‑1 (MQ=255)
                 aaCCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCATGCGGTGCGGCTCAACAGAGGGc                                                                                                                        >  8:43163/1‑149 (MQ=255)
                  aCCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCATGCGGTGCGGCTCAACAGAGGGc                                                                                                                        <  7:43163/148‑1 (MQ=255)
                  aCCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCATGCGGTGCGGCTCAACAGAGGGc                                                                                                                        <  3:213128/148‑1 (MQ=255)
                   ccAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCATGCGGTGCGGCTCAACAGAggg                                                                                                                         <  2:169361/146‑1 (MQ=255)
                       gCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGCCCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCATGCGGTGCGGCTCAACAGAGGGCCAAGGt                                                                                                                  >  8:265100/1‑149 (MQ=255)
                        ccTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCATGCGGTGCGGCTCAACAGAGGGCCAAGGTa                                                                                                                 <  4:100875/149‑1 (MQ=255)
                                   gCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCATGCGGTGCGGCTCAACAGAGGGCCAAGGTAGTCGGCAAACt                                                                                                      >  5:75337/1‑149 (MQ=255)
                                       agagTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCATGCGGTGCGGCTCAACAGAGGGCCAAGGTAGTCGGCAAACTgg                                                                                                    <  6:75337/147‑1 (MQ=255)
                                        gagTCTTCAACGCATCATTGGTCAGTTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCATGCGGTGCGGCTCAACAGAGGGCCAAGGTAGTCGGCAAACTGGGTg                                                                                                 <  7:204979/149‑1 (MQ=255)
                                          gTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCATGCGGTGCGGCTCAACAGAGGGCCAAGGTAGTCGGCAAACTGGGTGGa                                                                                               >  6:15396/1‑149 (MQ=255)
                                              tCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCATGCGGTGCGGc                                                                                                                                    <  3:143025/108‑1 (MQ=255)
                                              tCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCATGCGGTGCGGc                                                                                                                                    >  4:143025/1‑108 (MQ=255)
                                              tCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCATCCGGTGCGGCTCAACAGAGTGCCAAGGTAGTCGGCAACCTGGGTGGACAgg                                                                                           >  7:98578/1‑149 (MQ=255)
                                                aaCGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCATGCGGTGCGGCTCAACAGAGGGCCAAGGTAGTCGGCAAACTGGGTGGACAGGCt                                                                                         >  7:203851/1‑149 (MQ=255)
                                                        atTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCATGCGGTGCGGCTCAACAGAGGGCCAAGGTAGTCGGCAAACTGGGTGGACAGGCTGGTGGTAc                                                                                 <  5:15396/149‑1 (MQ=255)
                                                            gTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCATGTGGTGCGGCTCAACAGAGGGCCAAGGTAGTCGGCAAACTGGGTGGACAGGCTGGTGGTACCAAt                                                                             >  5:151635/1‑149 (MQ=255)
                                                            gTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCATGTGGTGCGGCTCAACAGAGGGCCAAGGTAGTCGGCAAACTGGGTGGACAGGCTGGTGGTACCAAt                                                                             >  5:79701/1‑149 (MQ=255)
                                                            gTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCATGCGGTGCGGCTCAACAGAGGGcc                                                                                                                       >  4:100005/1‑107 (MQ=255)
                                                            gTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCATGCGGTGCGGCTCAACAGAGGGcc                                                                                                                       <  3:100005/107‑1 (MQ=255)
                                                                 cttggcttgCTGCACGCCATTTACGTACAGGGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCATGCGGTGCGGCTCAACAGAGGGCCAAGGTAGTCGGCAAACTGGGTGGACAGGCTGGTGGTACCAATGCgg                                                                         >  1:29849/1‑148 (MQ=255)
                                                                     gcttgCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCATGCGGTGCGGCTCAACAGAGGGCCAAGGTAGTCGGCAAACTGGGTGGACAGGCTGGTGGTACCAATGCGGAAttt                                                                    <  2:29849/149‑1 (MQ=255)
                                                                         gctgcACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCATGCGGTGCGGCTCAACAGAGGGCCAAg                                                                                                                    <  5:32147/97‑1 (MQ=255)
                                                                         gctgcACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCATGCGGTGCGGCTCAACAGAGGGCCAAg                                                                                                                    >  6:32147/1‑97 (MQ=255)
                                                                          ctgcACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCATGCGGTGCGTCTCAACAGAGGGCCAAGGTAGTCGGCAAACTGGGTGGACAGGCTGGTGGTACCAATGCGGAATTTGCCa                                                                <  6:44147/148‑1 (MQ=255)
                                                                           tgcACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCATGCGGTGCGGCTCAACAGAGGGCCAAGGt                                                                                                                  <  3:106411/97‑1 (MQ=255)
                                                                              aCGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCATGCGGTGCGGCGCAACAGAGGGCCAAGGTAGTCGGCAAACTGGGTGGACAGGCTGGTGTTACCAATGCGGAATTTGCCAGCCGa                                                           <  4:276918/149‑1 (MQ=255)
                                                                               cGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCATGCGGTGCGGCTCAACAGAGGGCCAAGGTAGTCGGCAAACTGGGTGGACAGGCTGGTGGTACCAATGCGGAATTTGCCAGCCGAg                                                          <  8:204452/149‑1 (MQ=255)
                                                                                              gCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCATGTGGTGCGGCTCAACAGAGGGCCAAGGTAGTCGGCAAACTGGGTGGACAGGCTGGTGGTACCAATGCGGAATTTGCCAGCCGAGGACGCCAAGAACCAg                                           <  6:151635/149‑1 (MQ=255)
                                                                                              gCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCATGTGGTGCGGCTCAACAGAGGGCCAAGGTAGTCGGCAAACTGGGTGGACAGGCTGGTGGTACCAATGCGGAATTTGCCAGCCGAGGACGCCAAGAACCAg                                           <  6:79701/149‑1 (MQ=255)
                                                                                                  cATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCATGCGGTGCGGCTCAACAGAGGGCCAAGGTAGTCGGCAAACTGGGTGGACAGGCTGGTGGTACCAATGc                                                                           <  1:313006/113‑1 (MQ=255)
                                                                                                  cATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCATGCGGTGCGGCTCAACAGAGGGCCAAGGTAGTCGGCAAACTGGGTGGACAGGCTGGTGGTACCAATGc                                                                           >  2:313006/1‑113 (MQ=255)
                                                                                                     ggTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCATGCGGTGCGGCTCAACAGAGGGCCAAGGTAGTCGGCAAACTGGGTGGACAGGCTGGTGGTACCAATGCGGAATTTGCCAGCCGAGGACGCCAAGAACCAGTAGCttt                                    >  6:210641/1‑149 (MQ=255)
                                                                                                       tACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCATGCGGTGCGGCTCAACAGAGGGCCAAGGTAGTCGGCAAACTGGGTGGACAGGCTGGTGGTACCAATGCGGAATTTGCCAGCCGAGGACGCCAAGAACCAGTAGCTTTGc                                  >  1:237124/1‑149 (MQ=255)
                                                                                                                 gcgGTCGTATTTGAGGATCAAGCGCAGCCATGCGGTGCGGCTCAACAGAGGGCCAAGGTAGTCGGCAAACTGGGTGGACAGGCTGGTGGTACCAATGCGGAATTTGCCAGCCGAGGACGCCAAGAACCAGTAGCTTTGc                                  <  1:169602/139‑1 (MQ=255)
                                                                                                                 gcgGTCGTATTTGAGGATCAAGCGCAGCCATGCGGTGCGGCTCAACAGAGGGCCAAGGTAGTCGGCAAACTGGGTGGACAGGCTGGTGGTACCAATGCGGAATTTGCCAGCCGAGGACGCCAAGAACCAGTAGCTTTGc                                  >  2:169602/1‑139 (MQ=255)
                                                                                                                 gcgGTCGTATTTGAGGATCAAGCGCAGCCATGCGGTGCGGCTCAACAGAGGGCCAAGGTAGTCGGCAAACTGGGTGGACAGGCTGGTGGTACCAATGCGGAATTTGCCAGCCGAGGACGCCAAGAACCAGTAGCTTTGCCCAGATCCAg                        >  2:280458/1‑149 (MQ=255)
                                                                                                                       gTATTTGAGGATCAAGCGCAGCCATGCGGTGCGGCTCAACAGAGGGCCAAGGTAGTCGGCAAACTGGGTGGACAGGCTGGTGGTACCAATGCGGAATTTGCCAGCCGAGGACGCCAAGAACCAGTAGCTTTGCCCAGATCCAgg                       >  3:180653/1‑144 (MQ=255)
                                                                                                                       gTATTTGAGGATCAAGCGCAGCCATGCGGTGCGGCTCAACAGAGGGCCAAGGTAGTCGGCAAACTGGGTGGACAGGCTGGTGGTACCAATGCGGAATTTGCCAGCCGAGGACGCCAAGAACCAGTAGCTTTGCCCAGATCCAgg                       <  4:180653/144‑1 (MQ=255)
                                                                                                                         aTTTGAGGATCAAGCGCAGCCATGCGGTGCGGCTCAACAGAGGGCCAAGGTAGTCGGCAAACTGGGTGGACAGGCTGGTGGTACCAATGCGGAATTTGCCAGCCGAGGACGCCAAGAACCAGTAGCTTTGCCCAGATCCAGGCAACAcc                >  5:99738/1‑149 (MQ=255)
                                                                                                                            tGAGGATCAAGCGCAGCCATGCGGTGCGGCTCAACAGAGGGCCAAGGTAGTCGGCAAACTGGGTGGACAGGCTGGTGGTACCAATGCGGAATTTGCCAGCCGAGGACGCCAAGAACCAGTAGCTTTGCCCAGATCCAGGCAACACCGaa             >  6:121771/1‑149 (MQ=255)
                                                                                                                                gATCAAGCGCAGCCATGCGGTGCGGCTCAACAGAGGGCCAAGGTAGTCGGCAAACTGGGTGGACAGGCTGGTGGTACCAATGCGGAATTTGCCAGCCGAGGACGCCAAGAACCAGTAGCTTTGCCCAGATCCAGGCAACACCGAActgc         >  1:204399/1‑149 (MQ=255)
                                                                                                                                     aGCGCAGCCATGCGGTGCGGCTCAACAGAGGGCCAAGGTAGTCGGCAAACTGGGTGGACAGGCTGGTGGTACCAATGCGGAATTTGCCAGCCGAGGACGCCAAGAACCAGTAGCTTTGCCCAGATCCAGGCAACACCGAACTGCTGGCa    <  6:99738/149‑1 (MQ=255)
                                                                                                                                      gcgcAGCCATGCGGTGCGGCTCAACAGAGGGCCAAGGTAGTCGGCAAACTGGGTGGACAGGCTGGTGGTACCAATTCGGAATTTGCCAGCCGAGGACGCCAAGAACCAGTAGCTTTGCCCAGATCCAGGCAACACCGAACTGCTGGCAg   <  1:280458/149‑1 (MQ=255)
                                                                                                                                      gcgcAGCCATGCGGTGCGGCTCAACAGAGGGCCAAGGTAGTCGGCAAACTGGGTGGACAGGCTGGTGGTACCAATGCGGAATTTGCCAGCCGAGGACGCCAAGAACCAGTAGCTTTGCCCAGATCCAGGCAACACCGAACTGCTGGCAg   <  7:319479/149‑1 (MQ=255)
                                                                                                                                       cgcAGCCATGCGGTGCGGCTCAACAGAGGGCCAAGGTAGTCGGCAAACTGGGTGGACAGGCTGGTGGTACCAATGCGGAATTTGCCAGCCGAGGACGCCAAGAACCAGTAGCTTTGCCCAGATCCAGGCAACACCGAACTGCTGGCAgg  <  3:330671/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                          
GTGACCAGGCCACCCATAACCAGGCCTGGGGCGAGGCTTAGAGTCTTCAACGCATCATTGGTCAGCTTGGCTTGCTGCACGCCATTTACGTACAGCGTCATGGTACCGACCGCGCGGTCGTATTTGAGGATCAAGCGCAGCCATGCGGTGCGGCTCAACAGAGGGCCAAGGTAGTCGGCAAACTGGGTGGACAGGCTGGTGGTACCAATGCGGAATTTGCCAGCCGAGGACGCCAAGAACCAGTAGCTTTGCCCAGATCCAGGCAACACCGAACTGCTGGCAGG  >  NC_002947/4379882‑4380165

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