Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,016,329 A→T 9.7% I345N (ATC→AAC) ‡ PP_4421 ← aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,016,3290AT 9.7% 111.6 / 2.3 41I345N (ATC→AAC) ‡PP_4421aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (21/16);  new base T (2/2);  total (23/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGATCCGCTACGAACTCAATGGCGGCCATCATGCCCACACCACGCACTTCGCCGACGATTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACCAACGCGATCAGCGCCTTCCTCCATCACCTTATACACCTTCTCACCAACAATCGAAGCCGACAGCGGTGCGTA  >  NC_002947/5016181‑5016474
                                                                                                                                                    |                                                                                                                                                 
cGGATCCGCTACGAACTCAATGGCGGCCATCATGCCCACACCACGCACTTCGCCGACGATTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGCCGTGCGCTGGAATGTTCTCACCCTcgacga                                                                                                                                                   >  7:138685/1‑149 (MQ=255)
       gCTACGAACTCAATGGCGGCCATCATGCCCACACCACGCACTTCGCCGACGATTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACATCGTCGAg                                                                                                                                            >  2:179379/1‑149 (MQ=255)
        cTACGAACTCAATGGCGGCCATCATGCCCACACCACGCACTTCGCCGACGATTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACATCGTCGAgg                                                                                                                                           <  1:179379/149‑2 (MQ=255)
            gAACTCAATGGCGGCCATCATGCCCACACCACGCACTTCGCCGACGATTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAg                                                                                                                                       >  1:221714/1‑149 (MQ=255)
                tCAATGGCGGCCATCATGCCCACACCACGCACTTCGCCGACGATTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACATCGTCGAgggtga                                                                                                                                       <  4:140392/145‑6 (MQ=255)
                tCAATGGCGGCCATCATGCCCACACCACGCACTTCGCCGACGATTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACATCGTCGAgggtga                                                                                                                                       >  3:140392/1‑140 (MQ=255)
                   aTGGCGGCCATCATGCCCACACCACGCACTTCGCCGACGATTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTtgc                                                                                                                                >  4:209117/1‑149 (MQ=255)
                           catcatGCCCACACCACGCACTTCGCCGACGATTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGa                                                                                                                        <  4:171778/149‑1 (MQ=255)
                                   ccACACCACGCACTTCGCCGACGATTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGa                                                                                                                  >  1:49344/1‑147 (MQ=255)
                                      caccacGCACTTCGCCGACGATTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAggg                                                                                                             <  2:221714/149‑1 (MQ=255)
                                         cacGCACTTCGCCGACGATTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCcgc                                                                                                                          <  1:243339/133‑1 (MQ=255)
                                         cacGCACTTCGCCGACGATTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCcgc                                                                                                                          >  2:243339/1‑133 (MQ=255)
                                               cTTCGCCGACGATTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAg                                                                                                    >  5:299306/1‑149 (MQ=255)
                                                ttCGCCGACGATTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAg                                                                                                    >  6:257883/1‑148 (MQ=255)
                                                  cGCCGACGATTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGt                                                                                                   <  3:209117/147‑1 (MQ=255)
                                                   gCCGACGATTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAg                                                                                                 <  6:49202/148‑1 (MQ=255)
                                                     cgacgaTTGGCAGCGAGGAGAAGGTCTCCTGGATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCAGAGCGTTTGCTGCCGGGACACCGATAGGGTGGCCCGAGTAGgtg                                                                                              <  2:49344/149‑1 (MQ=255)
                                                          aTTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGcc                                                                                         <  6:299306/149‑1 (MQ=255)
                                                              gCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTgg                                                                                     >  5:355963/1‑149 (MQ=255)
                                                              gCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTgg                                                                                     >  3:162462/1‑149 (MQ=255)
                                                              gCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTgg                                                                                     <  5:257883/149‑1 (MQ=255)
                                                              gCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTgg                                                                                     >  1:346183/1‑149 (MQ=255)
                                                                aGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGga                                                                                   <  4:162462/149‑1 (MQ=255)
                                                                    aggagAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAg                                                                               <  2:346183/149‑1 (MQ=255)
                                                                             tctcCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCcgc                                                                                                                          >  3:202416/1‑97 (MQ=255)
                                                                             tctcCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCcgc                                                                                                                          <  4:202416/97‑1 (MQ=255)
                                                                             tctcCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACCAAcg                                                                      <  6:355963/149‑1 (MQ=255)
                                                                                 cTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACCAACGCGAt                                                                  >  4:186969/1‑149 (MQ=255)
                                                                                 cTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACCAACGCGAt                                                                  <  3:186969/149‑1 (MQ=255)
                                                                                   gCATCGACTTTTGGAAATACTGACAAACCTTGTCGAAGTGCTCTGGAATGTTCTCACCCTCGACGATGTCGAGCACATCGTTTGCTGCCTCGAAAAAGATAGGGTGGCACGAGTAGgtg                                                                                              >  5:355758/1‑119 (MQ=255)
                                                                                      tCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACCAACGCGATCAgcg                                                             >  7:302491/1‑149 (MQ=255)
                                                                                         aCTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGga                                                                                   >  3:157818/1‑124 (MQ=255)
                                                                                         aCTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGga                                                                                   <  4:157818/124‑1 (MQ=255)
                                                                                         aCTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACCAACGCGATCAGCGCCt                                                          >  7:314005/1‑149 (MQ=255)
                                                                                           ttttGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACCAACGCGATCAGCGCCTtc                                                        >  4:111724/1‑149 (MQ=255)
                                                                                              tGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACCAACGCGATCAGCGCCTtcct                                                      <  8:138685/148‑1 (MQ=255)
                                                                                                              ccTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACCAACGCGATCAGCGCCTTCCTCCATCACCTTATACAcc                                     >  5:257128/1‑149 (MQ=255)
                                                                                                                ttGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACCAACGCGATCAGCGCCTTCCTCCATCACCTTATACACCtt                                   >  6:19055/1‑149 (MQ=255)
                                                                                                                              tGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACCAACGCGATCAGCGCCTTCCTCCATCACCTTATACACCTTCTCACCAACAATc                      >  5:308113/1‑148 (MQ=255)
                                                                                                                                            ccTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACCAACGCGATCAGCGCCTTCCTCCATCACCTTATACACCTTCTCa                               >  4:133353/1‑125 (MQ=255)
                                                                                                                                            ccTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACCAACGCGATCAGCGCCTTCCTCCATCACCTTATACACCTTCTCa                               <  3:133353/125‑1 (MQ=255)
                                                                                                                                            ccTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACCAACGCGATCAGCGCCTTCCTCCATCACCTTATACACCTTCTCACCAACAATCGAAGCCGACAGCGGt       >  5:334914/1‑149 (MQ=255)
                                                                                                                                                gacgaTGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACCAACGCGATCAGCGCCTTCCTCCATCACCTTATACACCTTCTCACCAACAATCGAAGCCGACAGCGGTGCGt   <  6:334914/149‑1 (MQ=255)
                                                                                                                                                 acgaTGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACCAACGCGATCAGCGCCTTCCTCCATCACCTTATACACCTTCTCACCAACAATCGAAGCCGACAGCGGTGCGTa  <  5:19055/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CGGATCCGCTACGAACTCAATGGCGGCCATCATGCCCACACCACGCACTTCGCCGACGATTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACCAACGCGATCAGCGCCTTCCTCCATCACCTTATACACCTTCTCACCAACAATCGAAGCCGACAGCGGTGCGTA  >  NC_002947/5016181‑5016474

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