Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,383,860 T→A 14.2% D69V (GAC→GTC)  PP_2985 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,383,8600TA14.2% 99.2 / 5.2 41D69V (GAC→GTC) PP_2985hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (17/18);  new base A (2/4);  total (19/22)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAACCCCAGCATCTGAGTGATCGTCAGGTGGTCGCGGGTACAAACATGGAGCTTGTGCGTACGAATCCCATCCACATTGCGCCGGTAGAAATGATGGCCCTCTGACAGATCACTTCCCCGTACATA‑CCCCAAAGCCATCAACACCTCGTCCCGCGTCGAGGTATTGCGATGCGCGTGAACCTCAA‑TGAGCACATCCAGTTCGGGCTTC‑GCCGCCAGCCCTGGGACAGACGTGCTTCCAACGTGGTGGATAGCAATCAGTTCAGCGTCGAATGCCGAGGCTATGAGTGGCTTGTCTGT  >  NC_002947/3383712‑3384008
                                                                                                                                                     |                                                                                                                                                      
gAACCCCAGCATCTGAGTGATCGTCAGGTGGTCGCGGGTACAAACATGGAGCTTGTGCGTACGAATCCCATCCACATTGCGCCGGTAGAAATGATGGCCCTCTGACAGATCACTTCCCCGTACATA‑CCCCAAAGCCATCAACACCTCGt                                                                                                                                                        >  5:328507/1‑149 (MQ=255)
   ccccAGCATCTGAGTGATCGTCAGGTGGTCGCGGGTACAAACATGGAGCTTGTGCGTACGAATCCCATCCACATTGCGCCGGTAGAAATGATGGCCCTCTGACAGATCACTTCCCCGTACATA‑CCCCAAAGCCATCAACACCTCGTccc                                                                                                                                                     >  1:219318/1‑149 (MQ=255)
   ccccAGCATCTGAGTGATCGTCAGGTGGTCGCGGGTACAAACATGGAGCTTGTGCGTACGAATCCCATCCACATTGCGCCGGTAGAAATGATGGCCCTCTGACAGATCACTTCCCCGTACATA‑CCCCAAAGCCATCAACACCTCGTccc                                                                                                                                                     <  4:341619/149‑1 (MQ=255)
     ccAGCATCTGAGTGATCGTCAGGTGGTCGCGGGTACAAACATGGAGCTTGTGCGTACGAATCCCATCCACATTGCGCCGGTAGAAATGATGGCCCTCTGACAGATCACTTCCCCGTACATA‑CCCCAAAGCCATCAACACCTCGTCCcgc                                                                                                                                                   >  2:295784/1‑149 (MQ=255)
      cAGCATCTGAGTGATCGTCAGGTGGTCGCGGGTACAAACATGGAGCTTGTGCGTACGAATCCCATCCACATTGCGCCGGTAGAAATGATGGCCCTCTGACAGATCACTTCCCCGTACATA‑CCCCAAAGCCATCAACACCTCGTCCcgcg                                                                                                                                                  >  5:270637/1‑149 (MQ=255)
          aTCTGAGTGATCGTCAGGTGGTCGCGGGTACAAACATGGAGCTTGTGCGTACGAATCCCATCCACATTGCGCCGGTAGAAATGATGGCCCTCTGACAGATCACTTCCCCGTACATA‑CCCCAAAGCCATCAACACCTCGTCCCGCGTc                                                                                                                                                >  4:105314/1‑147 (MQ=255)
          aTCTGAGTGATCGTCAGGTGGTCGCGGGTACAAACATGGAGCTTGTGCGTACGAATCCCATCCACATTGCGCCGGTAGAAATGATGGCCCTCTGACAGATCACTTCCCCGTACATA‑CCCCAAAGCCATCAACACCTCGTCCCGCGTCGa                                                                                                                                              <  1:1577/149‑1 (MQ=255)
                 tGATCGTCAGGTGGTCGCGGGTACAAACATGGAGCTTGTGCGTACGAATCCCATCCACATTGCGCCGGTAGAAATGATGGCCCTCTGACAGATCACTTCCCCGTACATA‑CCCCAAAGCCATCAACACCTCGACGCGGGACGAGGTGTTg                                                                                                                                       <  3:202334/149‑1 (MQ=255)
                    tCGTCAGGTGGTCGCGGGTACAAACATGGAGCTTGTGCGTACGAATCCCATCCACATTGCGCCGGTAGAAATGATGGCCCTCTGACAGATCACTTCCCCGTACATA‑CCCCAAAGCCATCAACACCTCGTCCCGCGTCGAGGTATTGCGa                                                                                                                                    <  3:56156/149‑1 (MQ=255)
                      gTCAGGTGGTCGCGGGTACAAACATGGAGCTTGTGCGTACGAATCCCATCCACATTGCGCCGGTAGAAATGATGGCCCTCTGACAGATCACTTCCCCGTACATA‑CCCCAAAGCCATCAATACCTCGACGCGGGACGAGGTGTTgatggc                                                                                                                                  <  2:373773/149‑6 (MQ=255)
                      gTCAGGTGGTCGCGGGTACAAACATGGAGCTAGTGCGTACGAATCACATCCACATTGCGCAGGTAGAAATGAAGGCCCTCTGACAGATAACTTCCCCGTACATA‑CCCCAAAGCCATCAAAACATCGACACGCGTCGAGGTATTACGAtc                                                                                                                                  >  3:1746/1‑148 (MQ=255)
                         aGGTGGTCGCGGGTACAAACATGGAGCTTGTGCGTACGAATCCCATCCACATTGCGCCGGTAGAAATGATGGCCCTCTGACAGATCACTTCCCCGTACATA‑CCCCAAAGCCATCAACACCTCGTCCCGCGTCGAGGTATTGCGATgcgc                                                                                                                               >  2:197192/1‑149 (MQ=255)
                              gtCGCGGGTACAAACATGGAGCTTGTGCGTACGAATCCCATCCACATTGCGCCGGTAGAAATGATGGCCCTCTGACAGATCACTTCCCCGTACATA‑CCCCAAAGCCATCAACACCTCGTCCCGCGTCGAGGTATTGCGATGCGCGTGaa                                                                                                                          <  8:104840/149‑1 (MQ=255)
                                cgcgGGTACAAACATGGAGCATGTGCGTACGAATCCCATCCACATTGCGCCGGTAGAAATGATGGCCCTCTGACAGATCACTTCCCCGTACATACCCCCAAAGCCATCAACACCTCGTCCCGCGTCGAGGTATTGCGATGCGCGTGAAc                                                                                                                         <  1:295784/149‑1 (MQ=255)
                                                                aTCCCATCCACATTGCGCCGGTAGAAATGATGGCCCTCTGACAGATCACTTCCCCGTACATA‑CCCCAAAGCCATCAACACCTCGTCCCGCGTCGAGGTATTGCGATGCGCGTGAACCTCAA‑TGAGCACATCCAGTTCGGGCTTc‑gccgc                                                                                      >  4:375443/1‑149 (MQ=255)
                                                                    cATCCACATTGCGCCGGTAGAAATGATGGCCCTCTGACAGATCACTTCCCCGTACATA‑CCCCAAAGCCATCAACACCTCGTCCCGCGTCGAGGTATTGCGATGCGCGTGAACCTCAA‑TGAGCACATCCAGTTCGGGCTTC‑GCCGCCAGc                                                                                  >  2:85257/1‑149 (MQ=255)
                                                                       ccACATTGCGCCGGTAGAAATGATGGCCCTCTGACAGATCACTTCCCCGTACATA‑CCCCAAAGCCATCAACACCTCGTCCCGCGTCGAGGTATTGCGATGCGCGTGAACCTCAA‑TGAGCACATCCAGTTCGGGCTTC‑GCCGCCAGCCCt                                                                               >  1:321073/1‑149 (MQ=255)
                                                                           aTTGCGCCGGTAGAAATGATGGCCCTCTGACAGATCACTTCCCCGTACATA‑CCCCAAAGCCATCAACACCTCGTCCCGCGTCGAGGTATTGCGATGCGCGTGAACCTCAA‑TGAGCACATCCAGTTCGGGCTTC‑GCCGCCAGCCCTGGGa                                                                           <  2:219318/149‑1 (MQ=255)
                                                                              gcgcCGGTAGAAATGATGGCCCTCTGACAGATCACTTCCCCGTACATA‑CCCCAAAGCCATCAACACCTCGTCCCGCGTCGAGGTATTGCGATGCGCGTGAACCTCAA‑TGAGCACATCCAGTTCGGGCTTC‑GCCGCCAGCCCTGGGACa                                                                         <  1:197192/148‑1 (MQ=255)
                                                                                    gTAGAAATGATGGCCCTCTGACAGATCACTTCCCCGTACATA‑CCCCAAAGCCATCAACACCTCGTCCCGCGTCGAGGTATTGCGATGCGCGTGAACCTCAA‑TGAGCACATCCAGTTCGGGCTTT‑GCCGCCAGCCCTGGGACAGACGTGc                                                                  >  1:129836/1‑149 (MQ=255)
                                                                                    gTAGAAATGATGGCCCTCTGACAGATCACTTCCCCGTACATA‑CCCCAAAGCCATCAACACCTCGTCCCGCGTCGAGGTATTGCGATGCGCGTGAACCTCAA‑TGAGCACATCCAGTTCGGGCTTC‑GCCGCCAGCCCTGGGACAGACGTGc                                                                  <  3:105314/149‑1 (MQ=255)
                                                                                      aGAAATGATGGCCCTCTGACAGATCACTTCCCCGTACATA‑CCCCAAAGCCATCAACACCTCGTCCCGCGTCGAGGTATTGCGATgcgc                                                                                                                               <  3:209132/88‑1 (MQ=255)
                                                                                      aGAAATGATGGCCCTCTGACAGATCACTTCCCCGTACATA‑CCCCAAAGCCATCAACACCTCGTCCCGCGTCGAGGTATTGCGATgcgc                                                                                                                               >  4:209132/1‑88 (MQ=255)
                                                                                             atgGCCCTCTGACAGATCACTTCCCCGTACATA‑CCCCAAAGCCATCAACACCTCGTCCCGCGTCGAGGTATTGCGATGCGCGTGAACCTCAA‑TGAGCACATCCAGTTCGGGCTTC‑GCCGCCAGCCCTGGGACAGACGTGCTTCCAACgt                                                         <  5:156481/149‑1 (MQ=255)
                                                                                                 cccTCTGACAGATCACTTCCCCGTACATA‑CCCCAAAGCCATCAACACCTCGTCCCGCGTCGAGGTATTGCGATGCGCGTGAACCTCAA‑TGAGCACATCCAGTTCGGGCTTC‑GCCGCCAGCCCTGGGACAGACGTGCTTCCAACgtggtg                                                     <  6:270637/149‑1 (MQ=255)
                                                                                                   ctctGACAGATCACTTCCCCGTACATA‑CCCCAAAGCCATCAACACCTCGTCCCGCGTCGAGGTATTGCGATGCGCGTGAACCTCAA‑TGAGCACATCCAGTTCGGGCTTC‑GCCGCCAGCCCTGGGACAGACGTGCTTCCAACGTGGTTGa                                                   <  7:172709/149‑1 (MQ=255)
                                                                                                          aGATCACTTCCCCGTACATA‑CCCCAAAGCCATCAACACCTCGTCCCGCGTCGAGGTATTGCGATGCGCGTGAACCTCAA‑TGAGCACATCCAGTTCGGGCTTC‑GCCGCCAGCCCTGGGACAGACGTGCTTCCAACGTGGTGGATAGCAAt                                            >  7:219944/1‑149 (MQ=255)
                                                                                                           gATCACTTCCCCGTACATA‑CCCCAAAGCCATCAACACCTCGTCCCGCGTCGAGGTATTGCGATGCGCGTGAACCTCAA‑TGAGCACATCCAGTTCGGGCTTCTGCCGCCAGCCCTGGGACAGACGTGCTTCCAACGTGGTGGATAGCAAt                                            <  1:85257/149‑1 (MQ=255)
                                                                                                               aCTTCCCCGTACATA‑CCCCAAAGCCATCAACACCTCGTCCCGCGTCGAGGTATTGCGATGCGCGTGAACCTCAA‑TGAGCACATCCAGTTCGGGCTTC‑GCCGCCAGCCCTGGGACAGACGTGCTTCCAACGTGGTGGATAg                                                <  8:44309/140‑1 (MQ=255)
                                                                                                               aCTTCCCCGTACATA‑CCCCAAAGCCATCAACACCTCGTCCCGCGTCGAGGTATTGCGATGCGAGTGAACCTCAA‑TGAGCACATCCAGTTCGGGCTTC‑GCCGCCAGCCCTGGGACAGACGTGCTTACAACGTGGTGGATAg                                                >  7:44309/1‑140 (MQ=255)
                                                                                                                cTTCCCCGTACATA‑CGCGAAAGACTTCAACACCTCGTCCCGCGTGGAGGTATTGCGATGCGCGTGAACCTCAA‑TGAGCACATCCAGTTCGGGCTTC‑GCCGCCAGCCGTGGGACAGACGTGGTTCCTACGTGGTGGATAGCAATCAGtt                                       <  4:1746/148‑1 (MQ=255)
                                                                                                                 ttCCCCGTACATA‑CCCCAAAGCCATCAACACCACGTCCCGCGTCGAGGTATTGCGATGCGCGTGAACCTCAA‑TGAGCACATCCAGTTCGGGCTTC‑GCCGCCAGCCCTGGGACAGACGTGCTTCCAACGTGGTGGATAGCAATCAGTTCa                                     >  1:242187/1‑149 (MQ=255)
                                                                                                                     ccGTACATA‑CCCCAAAGCCATCAACACCTCGTCCCGCGTCGAGGTATTGCGATGCGCGTGAACCTCAA‑TGAGCACATCCAGTTCGGGCTTC‑GCCGCCAGCCCTGGGACAGACGTGCTTCCAACGTGGTGGATa                                                 >  4:299000/1‑133 (MQ=255)
                                                                                                                     ccGTACATA‑CCCCAAAGCCATCAACACCTCGTCCCGCGTCGAGGTATTGCGATGCGCGTGAACCTCAA‑TGAGCACATCCAGTTCGGGCTTC‑GCCGCCAGCCCTGGGACAGACGTGCTTCCAACGTGGTGGATa                                                 <  3:299000/133‑1 (MQ=255)
                                                                                                                     ccGTACATA‑CCCCAAAGCCATCAACACCTCGTCCCGCGTCGAGGTATTGCGATGCGCGTGAACCTCAA‑TGAGCACATCCAGTTCGGGCTTC‑GCCGCCAGCCCTGGGACAGACGTGCTTCCAACGTGGTGGATAGCAATCAGTTCAGCGt                                 >  1:344778/1‑149 (MQ=255)
                                                                                                                                      gCCATCAACACCTCGTCCCGCGTCGAGGTATTGCGATGCGCGTGAACCTCAA‑TGAGCACATCCAGTTCGGGCTTC‑GCCGCCAGCCCTGGCACAGACGTGCTTCCAACGTGGTGGATAGCAATCAGTTCAGCGTCGAATGCCGAGGCTAt                 >  7:106034/1‑149 (MQ=255)
                                                                                                                                         aTCAACACCTCGTCCCGCGTCGAGGTATTGCGATGCGCGTGAACCTCAA‑TGAGCACATCCAGTTCGGGCTTA‑GCCGCCAGACCTGGGACAGACGTGCTTCCAACGTGGTGGATAGCAATCAGTTCAGCGTCGAATGCCGAGGCTATGAg              >  7:59001/1‑149 (MQ=255)
                                                                                                                                           cAATACCTCGACGCGGGACGAGGTGTTGCGATGCGCGTGAACCTCAACTGAGCACATCCAGTTCGGGCTTC‑GCCGCCa                                                                                    <  1:34937/78‑1 (MQ=21)
                                                                                                                                              tacCTCGACGCGGGACGAGGTGTTGCGATGCGCGTGAACCTCAA‑TGAGCACATCCAGTTCGGGCTTC‑GCCGCCAGCCCTGGGACAGACGTGCTTCCAACGTGGTGGATAGCAATCAGTTCAGCGTCGAATGCCGAGGCTATGAGTgg           >  2:115341/2‑147 (MQ=255)
                                                                                                                                              tacCTCGACGCGGGACGAGGTGTTGCGATGCGCGTGAACCTCAA‑TGAGCACATCCAGTTCGGGCTTC‑GCCGCCAGCCCTGGGACAGACGTGCTTCCAACGTGGTGGATAGCAATCAGTTCAGCGTCGAATGCCGAGGCTATGAGTgg           <  1:115341/146‑1 (MQ=255)
                                                                                                                                              cacCTCGTCCCGCGTCGAGGTATTGCGATGCGCGTGAACCTCAA‑TGAGCACATCCAGTTCGGGCTTC‑GCCGCCAGCCCTGGGACAGACGTGCTTCCAACGTGGTGGATAGCAATCAGTTCAGCGTCg                               <  7:2837/127‑1 (MQ=255)
                                                                                                                                                    gTCCCGCGTCGAGGTATTGCGATGCGCGTGAACCTCAA‑TGAGCACATCCAGTTCGGGCTTT‑GCCGCCAGCCCTGGGACAGACGTGCTTCCAACGTGGTGGATAGCAATCAGTTCAGCGTCGAATGCCGAGGCTATGAGTGGCTtgtc     <  2:129836/147‑1 (MQ=255)
                                                                                                                                                     tCCCGCGTCGAGGTATTGCGATGCGCGTGAACCTCAA‑TGAGCACATCCAGTTCGGGCTTC‑GCCGCCAGCCCTGGGACAGACGTGCTTCCAACGTGGTGGATAGCAATCAGTTCAGCGTCGAATGCCGAGGCTATGAGTGGCTtgtctgt  <  2:242187/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GAACCCCAGCATCTGAGTGATCGTCAGGTGGTCGCGGGTACAAACATGGAGCTTGTGCGTACGAATCCCATCCACATTGCGCCGGTAGAAATGATGGCCCTCTGACAGATCACTTCCCCGTACATA‑CCCCAAAGCCATCAACACCTCGTCCCGCGTCGAGGTATTGCGATGCGCGTGAACCTCAA‑TGAGCACATCCAGTTCGGGCTTC‑GCCGCCAGCCCTGGGACAGACGTGCTTCCAACGTGGTGGATAGCAATCAGTTCAGCGTCGAATGCCGAGGCTATGAGTGGCTTGTCTGT  >  NC_002947/3383712‑3384008

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: ‑