Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 4,012,374 G→A 9.4% A500T (GCC→ACC)  rbbA → fused ribosome‑associated ATPases

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090484,012,3740GA 9.4% 229.2 / 10.6 74A500T (GCC→ACC) rbbAfused ribosome‑associated ATPases
Reads supporting (aligned to +/- strand):  ref base G (28/39);  new base A (5/2);  total (33/41)
Fisher's exact test for biased strand distribution p-value = 2.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

ATCTCACTGATGCACGCCGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCGGCAGGGCAGAGCAACGAAGCCGA  >  AP009048/4012305‑4012441
                                                                     |                                                                    
aTCTCACTGATGCACGCCGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAgc                                                                     >  1:1100094/1‑71 (MQ=255)
   tcACTGATGCACGCCGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAaccgc                                                                  >  1:2044593/1‑71 (MQ=255)
     aCTGATGCACGCCGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAg                                                                      <  1:1514773/65‑1 (MQ=255)
     aCTGATGCACGCCGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCa                                                                <  1:3865230/71‑1 (MQ=255)
      cTGATGCACGCCGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAACCGCCAg                                                               >  1:2581693/1‑71 (MQ=255)
       tGATGCACGCCGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGt                                                              >  1:4309698/1‑71 (MQ=255)
       tGATGCACGCCGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGt                                                              >  1:2481405/1‑71 (MQ=255)
       tGATGCACGCCGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGt                                                              >  1:2216268/1‑71 (MQ=255)
            cACGCCGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGa                                                         <  1:4140233/71‑1 (MQ=255)
             aCGCCGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAgc                                                                     >  1:3034048/1‑58 (MQ=255)
             aCGCCGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTgg                                                          <  1:4430216/69‑1 (MQ=255)
              cGCCGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTgg                                                          >  1:3887074/1‑68 (MQ=255)
               gCCGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTgg                                                          <  1:4178913/67‑1 (MQ=255)
               gCCGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAACCGCCAGTCTGGaa                                                        >  1:3467875/1‑69 (MQ=255)
                 cGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAgg                                                    >  1:1073578/1‑71 (MQ=255)
                 cGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAgg                                                    <  1:488207/71‑1 (MQ=255)
                  ggAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAgg                                                    >  1:4021303/1‑70 (MQ=255)
                  ggAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGc                                                   <  1:92226/71‑1 (MQ=255)
                   gAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTgg                                                          <  1:1959841/63‑1 (MQ=255)
                   gAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCa                                                  >  1:4088991/1‑71 (MQ=255)
                    aaaaGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAagag                                                     <  1:2385783/67‑1 (MQ=255)
                      aaGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAgg                                                    <  1:1460204/66‑1 (MQ=255)
                      aaGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCAttt                                              <  1:4528557/71‑1 (MQ=255)
                      aaGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCAtt                                               >  1:156090/1‑70 (MQ=255)
                       aGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTa                                             >  1:1700671/1‑71 (MQ=255)
                          gCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCt                                                            >  1:2575656/1‑54 (MQ=255)
                           cTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAaga                                                      >  1:2282359/1‑59 (MQ=255)
                             tGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGcc                                        >  1:3475786/1‑70 (MQ=255)
                             tGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGcc                                        >  1:1043421/1‑70 (MQ=255)
                              gCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGcc                                        <  1:3643597/69‑1 (MQ=255)
                              gCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTa                                      >  1:2372235/1‑71 (MQ=255)
                              gCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAACCGCCAGTCTGGAAGAGGCATTTATCGCCTa                                      >  1:4371731/1‑71 (MQ=255)
                                cAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTAtt                                    <  1:4221126/71‑1 (MQ=255)
                                  gCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTAtt                                    >  1:1444732/1‑69 (MQ=255)
                                  gCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTg                                  <  1:3885545/71‑1 (MQ=255)
                                  gCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTg                                  <  1:2978109/71‑1 (MQ=255)
                                  gCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAACCGCCAGTCTGGAAGAGGCAttt                                              <  1:4536850/59‑1 (MQ=255)
                                   cGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTg                                  >  1:2680751/1‑70 (MQ=255)
                                   aggTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTg                                  >  1:1191591/2‑70 (MQ=255)
                                       acacCGCAGGAACTGGTTGAGAAACGCGGAACCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAgg                              >  1:333886/1‑70 (MQ=255)
                                        cacCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAg                               <  1:3256425/68‑1 (MQ=255)
                                         acCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAg                           >  1:3260163/1‑71 (MQ=255)
                                          ccGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGaa                                                        <  1:268521/42‑1 (MQ=255)
                                           cGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCAt                                                <  1:2933804/48‑1 (MQ=255)
                                             cAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCgg                        >  1:2158961/1‑70 (MQ=255)
                                                gAACTGGTTGAGAAACGCGGAACCGCCAGTCTGGAAGAGGCATTTATc                                           >  1:4416762/1‑48 (MQ=255)
                                                 aaCTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCggcaggg                   >  1:3121323/1‑71 (MQ=255)
                                                  aCTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAaga                                                      >  1:3124612/1‑36 (MQ=255)
                                                  aCTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGc                          >  1:3953749/1‑63 (MQ=255)
                                                  aCTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCggcaggg                   >  1:2688326/1‑70 (MQ=255)
                                                  aCTGGTTGAGAAACGCGGAGCCGCCAGTCTCGAAGAGGCATTTATCGCCTATTTGCAGGAAGCggcagggc                  <  1:4215182/71‑1 (MQ=255)
                                                   cTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTTATCGCCTATTTGCAGGAAGCggcaggg                   >  1:293597/1‑70 (MQ=255)
                                                   cTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCggcagggca                 <  1:2717548/71‑1 (MQ=255)
                                                   cTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCggcagggca                 <  1:1654888/71‑1 (MQ=255)
                                                   cTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCggcagggca                 <  1:2444069/71‑1 (MQ=255)
                                                   cTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCggcagggca                 <  1:3333810/71‑1 (MQ=255)
                                                     ggTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCggcagggca                 >  1:2895825/1‑69 (MQ=255)
                                                     ggTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCggcagggca                 >  1:602267/1‑69 (MQ=255)
                                                      gTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCggcagggcag                <  1:1798560/69‑1 (MQ=255)
                                                      gTTGAGAAACGCGGAACCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCggcagggcag                <  1:743188/69‑1 (MQ=255)
                                                       ttGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTg                                  <  1:345205/50‑1 (MQ=255)
                                                       ttGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCggcaggg                   <  1:4159733/65‑1 (MQ=255)
                                                       ttGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCGGCAGGGCAGAGc             <  1:1874177/71‑1 (MQ=255)
                                                        tGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCGGCAGGGCAGAGc             >  1:3004803/1‑70 (MQ=255)
                                                         gagaAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCGGCAGGGCAGAGCa            <  1:4051510/70‑1 (MQ=255)
                                                         gagaAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCGGCAGGGCAGAGCa            <  1:1832528/70‑1 (MQ=255)
                                                          agaAACGCGGAGCCGCCAGTCTGGAAGAGGCAttt                                              <  1:622479/35‑1 (MQ=255)
                                                           gaAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCt                                       <  1:3145767/41‑1 (MQ=255)
                                                           gaAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCggcaggg                   <  1:1663728/61‑1 (MQ=255)
                                                           gaAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCGGCAGGGCAGAGCAAcg         <  1:2753011/71‑1 (MQ=255)
                                                           gaAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCGGCAGGGCAGAGCAAcg         <  1:4184919/71‑1 (MQ=255)
                                                             aaCGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCGGCAGGGCAGAGCAAcgaa       <  1:3050729/71‑1 (MQ=255)
                                                             aaCGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCGGCAGGGCAGAGCAAc          >  1:3985870/1‑68 (MQ=255)
                                                                gcgGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCGGCAGGGCAGAGCAAcgaagcc    <  1:3019027/71‑1 (MQ=255)
                                                                  ggAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCGGCAGGGCAGAGCAAcgaagccga  <  1:943319/71‑1 (MQ=255)
                                                                  ggAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCGGCAGGGCAGAGCAAcgaagccga  <  1:1410362/71‑1 (MQ=255)
                                                                  ggAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCGGCAGGGCAGAGCAAcgaagccga  <  1:4021655/71‑1 (MQ=255)
                                                                  ggAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCGGCAGGGCAGAGCAAcgaagccga  <  1:2147559/71‑1 (MQ=255)
                                                                     gccgccAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCGGCAGGGCAGAGCAAcgaagccga  >  1:1906260/1‑68 (MQ=255)
                                                                     |                                                                    
ATCTCACTGATGCACGCCGGAAAAGTGCTTGCCAGCGGTACACCGCAGGAACTGGTTGAGAAACGCGGAGCCGCCAGTCTGGAAGAGGCATTTATCGCCTATTTGCAGGAAGCGGCAGGGCAGAGCAACGAAGCCGA  >  AP009048/4012305‑4012441

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