Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,151,401 C→T 9.3% I14I (ATC→ATT rluC → ribosomal large subunit pseudouridine synthase C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,151,4010CT 9.3% 113.6 / 2.5 43I14I (ATC→ATTrluCribosomal large subunit pseudouridine synthase C
Reads supporting (aligned to +/- strand):  ref base C (20/19);  new base T (2/2);  total (22/21)
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.81e-01

CCGGCAAAGAGCCTGGTGCCGCTGGCCTGCCGCTTTTGGCGGCGTTCGCGACTATAGCAGCAATGATTAAGTGCTTCAATTCCATAAAAAATTGTTATGATCCCGTAATGACGACCAATACCCCTCCGACTTCCGGCGTTCAGCTGATCGAAGTCGCGCCGGAGCTTGCCGGCCAACGCATCGATAATTTCCTCATCACGGCCCTCAAGGGCGTGCCCAAGACACTGGTCTACCGCATCCTGCGCAAGGGTGAGGTGAGGGTCAACAAGGGCCGCGTCAAGCCGGAG  >  NC_002947/2151253‑2151539
                                                                                                                                                    |                                                                                                                                          
ccGGCAAAGAGCCTGGTGCCGCTGGCCTGCCGCTTTTGGCGGCGTTCGCGACTATAGCAGCAATGATTAAGTGCTTCAATTCCATAAAAAATTGTTATGATCCCGTAATGACGACCAATACCCCTCCGACTTCCGGCGTTCAGCTGATCg                                                                                                                                           <  1:745391/150‑1 (MQ=255)
  ggCAAAGAGCCTGGTGCCGCTGGCCTGCCGCTTTTGGCGGCGTTCGCGACTATAGCAGCAATGATTAAGTGCTTCAATTCCATAAAAAATTGTTATGATCCCGTAATGACGACCAATACCCCTCCGACTTCCGGCGTTCAGCTGATCGaa                                                                                                                                         <  1:817579/150‑1 (MQ=255)
       agagCCTGGTGCCTTTGGCCTGCCGATTTTGGTGGCGTTCGCGACTATAGCAGCAATGATTAAGTGATTCAATTCCATAAAAAATTGTTATGAGCCCGTAATGACTACGAATGCCCCTCCGACTTCCGTCGTTCAGCTGATCGAACTcgc                                                                                                                                    <  2:572817/150‑1 (MQ=255)
       agagCCTGGTGCCGCTGGCCTGCCGCTTTTGGCGGCGTTCGCGACTATAGCAGCAATGATTAAGTGCTTCAATTCCATAAAAAATTGTTATGATCCCGTAATGACGACCAATACCCCTCCGACTTCCGGCGTTCAGCTGATCGAAGTcgc                                                                                                                                    >  1:1019696/1‑150 (MQ=255)
           ccTGGTGCCGCTGGCCTGCCGCTTTTGGCGGCGTTCGCGACTATAGCAGCAATGAATACTTGCTTCACTTCCATAAAAAATTGTTATGATCCCGTAATGACGACCAATACCCCTCCGACTTCCGGCGTTCAGCTGATCGAAGTCGCGCCg                                                                                                                                <  1:291677/150‑1 (MQ=255)
                          cTGCCGCTTTTGGCGGCGTTCGCGACTATAGCAGCAATGATTAAGTGCTTCAATTCCATAAAAAATTGTTATGATCCCGTAATGACGACCAATACCCCTCCGACTTCCGGCGTTCAGCTGATCGAAGTCGCGCCGGAGCTTGCCGGCCaa                                                                                                                 >  2:1126175/1‑150 (MQ=255)
                           tGCCGCTTTTGGCGGCGTTCGCGACTATAGCAGCAATGATTAAGTGCTTCAATTCCATAAAAAATTGTTATGATCCCGTAATGACGACCAATACCCCTCCGACTTCCGGCGTTCAGCTGATCGAAGTCGCGCCGGAGCTTGCCGGCCAAc                                                                                                                <  2:1286764/150‑1 (MQ=255)
                           tGCCGCTTTTGGCGGCGTTCGCGACTATAGCAGCAATGATTAAGTGCTTCAATTCCATAAAAAATTGTTATGATCCCGTAATGACGACCAATACCCCTCCGACTTCCGGCGTTCAGCTGATCGAAGTCGCGCCGGAGCTTGCCGGCCAAc                                                                                                                <  1:868900/150‑1 (MQ=255)
                              cGCTTTTGGCGGCGTTCGCGACTATAGCAGCAATGATTAAGTGCTTCAATTCCATAAAAAATTGTTATGATCCCGTAATGACGACCAATACCCCTCCGACTTCCGGCGTTCAGCTGATCGAAGTCGCGCCGGAGCTTGCCGGCCAACGCa                                                                                                             >  2:282394/1‑150 (MQ=255)
                                cTTTTGGCGGCGTTCGCGACTATAGCAGCAATGATTAAGTGCTTCAATTCCATAAAAAATTGTTATGATCCCGTAATGACGACCAATACCCCTCCGACTTCCGGCGTTCAGCTGATCGAAGTCGCGCCGGAGCTTGCCGGCCAACGCATc                                                                                                           >  2:192455/1‑150 (MQ=255)
                                 ttttGGCGGCGTTCGCGACTATAGCAGCAATGATTAAGTGCTTCAATTCCATAAAAAATTGTTATGATCCCGTAATGACGACCAATACCCCTCCGACTTCCGGCGTTCAGCTGATCGAAGTCGCGCCGGAGCTTGCCGGc                                                                                                                    >  2:1082070/1‑140 (MQ=255)
                                 ttttGGCGGCGTTCGCGACTATAGCAGCAATGATTAAGTGCTTCAATTCCATAAAAAATTGTTATGATCCCGTAATGACGACCAATACCCCTCCGACTTCCGGCGTTCAGCTGATCGAAGTCGCGCCGGAGCTTGCCGGc                                                                                                                    <  1:1082070/140‑1 (MQ=255)
                                    tGGCGGCGTTCGCGACTATAGCAGCAATGATTAAGTGCTTCAATTCCATAAAAAATTGTTATGATCCCGTAATGACGACCAATACCCCTCCGACTTCCGGCGTTCAGCTGATCGAAGTCGCGCCGGAGCTTGCCGGCCAACGCATCGATa                                                                                                       <  2:894684/150‑1 (MQ=255)
                                      gcggcgTTCGCGACTATAGCAGCAATGATTAAGTGCTTCAATTCCATAAAAAATTGTTATGATCCCGTAATGACGACCAATACCCCTCCGACTTCCGGCGTTCAGCTGATCGAAGTCGCGCCGGAGCTTGCCGGCCAACGCATCGATAAt                                                                                                     >  2:506584/1‑150 (MQ=255)
                                         gcgTTCGCGACTATAGCAGCAATGATTAAGTGCTTCAATTCCATAAAAAATTGTTATGATCCCGTAATGACGACCAATACCCCTCCGACTTCCGGCGTTCAGCTGATCGAAGTCGCGCCGGAGCTTGCCGGCCAACGCATCGATAATTTc                                                                                                  <  1:1126175/150‑1 (MQ=255)
                                               gcgACTATAGCAGCAATGATTAAGTGCTTCAATTCCATAAAAAATTGTTATGATCCCGTAATGACGACCAATACCCCTCCGACTTCCGGCGTTCAGCTGATCGAAGTCGCGCCGGAGCTTGCCGGCCAACGCCTCACTAATTTCCtcatc                                                                                            >  2:415575/1‑150 (MQ=255)
                                                 gACTATAGCAGCAATGATTAAGTGCTTCAATTCCATAAAAAATTGTTATGATCCCGTAATGACGACCAATACCCCTCCGACTTCCGGCGTTCAGCTGATCGAAGTCGCGCCGGAGCTTGCCGGCCAACGCATCGATAATTTCCTCATCAc                                                                                          >  2:498471/1‑150 (MQ=255)
                                                    tataGCAGCAATGATTAAGTGCTTCAATTCCATAAAAAATTGTTATGATCCCGTAATGACGACCAATACCCCTCCGACTTCCGGCGTTCAGCTGATCGAAGTCGCGCCGGAGCTTGCCGGCCAACGCATCGATAATTTCCTCATCACGGc                                                                                       <  1:192455/150‑1 (MQ=255)
                                                      taGCAGCAATGATTAAGTGCTTCAATTCCATAAAAAATTGTTATGATCCCGTAATGACGACCAATACCCCTCCGACTTCCGGCGTTCAGCTGATCGAAGTCGCGCCGGAGCTTGCCGGCCAACGCATCGATAATTTCCTCATCACGGccc                                                                                     >  1:465082/1‑150 (MQ=255)
                                                                        gCTTCAATTCCATAAAAAATTGTTATGATCCCGTAATGACGACCAATACCCCTCCGACTTCCGGCGTTCAGCTGATCGAAGTCGCGCCGGAGCTTGCCGGCCAACGCATCGATAATTTCCTCATCACGGCCCTCAAGGGCGTGCCCAAGa                                                                   >  2:778044/1‑150 (MQ=255)
                                                                        gCTTCAATTCCATAAAAAATTGTTATGATCCCGTAATGACGACCAATACCCCTCCGACTTCCGGCGTTCAGCTGATCGAAGTCGCGCCGGAGCTTGCCGGCCAACGCATCGATAATTTCCTCATCACGGCCCTCAAGGGCGTGCCCAAGa                                                                   <  2:1019696/150‑1 (MQ=255)
                                                                        gCTTCAATTCCATAAAAAATTGTTATGATCCCGTAATCACGACCAATACCCCTCCGACTTCCGGCGTTCAGCGGATCAAAGTCGCGCCGGAGCTTTCCGGCCAACGCATCGATAATTTCCTCATACCCGGCCTCAAGGGCGTGCCCAAGa                                                                   <  1:38913/150‑1 (MQ=255)
                                                                         cTTCAATTCCATAAAAAATTGTTATGATCCCGTAATGACGACCAATACCCCTCCGACTTCCGGCGTTCAGCTGATCGAAGTCGCGCCGGa                                                                                                                              <  1:1029929/90‑1 (MQ=255)
                                                                         cTTCAATTCCATAAAAAATTGTTATGATCCCGTAATGACGACCAATACCCCTCCGACTTCCGGCGTTCAGCTGATCGAAGTCGCGCCGGa                                                                                                                              >  2:1029929/1‑90 (MQ=255)
                                                                           tCAATTCCATAAAAAATTGTTATGATCCCGTAATGACGACCAATACCCCTCCGACTTCCGGCGTTCAGCTGATCGAAGTCGCGCCGGAGCTTGCCGGCCAACGCATCGATAATTTCCTCATCACGGCCCTCAAGGGCGTGCCCAAGacac                                                                <  1:506584/150‑1 (MQ=255)
                                                                             aaTTCCATAAAAAATTGTTATGATCCCGTAATGACGACCAATACCCCTCCGACTTCCGGCGTTCAGCTGATCGAAGTCGCGCCGGAGCTTGCCGGCCAACGCATCGATAATTTCCTCATCACGGCCCTCAAGGGCGTGCCCAAGACACTg                                                              >  1:335848/1‑150 (MQ=255)
                                                                               ttCCATAAAAAATTTTTATGATCCCGTAACGACGACCATTCCACCTACCATTTCCGGCGTTCAGCTGATCAAAGTCAGGCCGGAGCTTGCCGGCCAACGCATCGATAATTTCCTCATCACGGCCCTCAAGGGCGTGCCCAAGACACTGGt                                                            <  1:926443/150‑1 (MQ=255)
                                                                                   aTAAAAAATTGTTATGATCCCGTAATGACGACCAATACCCCTCCGACTTCCGGCGTTCAGCTGATCGAAGTCGCGCCGGAGCTTGCCGGCCAACGAATCGATAATTTCCTCATCACGGCCCTCAAGGGCGTGCCCAAGACACTGGTCTAc                                                        >  2:305434/1‑150 (MQ=255)
                                                                                                  gATCCCGTAATGACGACCAATACCCCTCCGACTTCCGGCGTTCAGCTGATCGAAGTCGCGCCGGAGCTTGCCGGCCAACGCATCGATAATTTCCTCATCACGGCCCTCAAGGGCGTGCCCAAGACACTGGTCTACCGCATCCTGCGCAAg                                         >  2:643/1‑150 (MQ=255)
                                                                                                             gacgacCAATACCCCTCCGACTTCCGGCGTTCAGCTGATCGAAGTCGCGCCGGAGCTTGCCGGCCAACGCATCGATAATTTCCTCATCACGGCCCTCAAGGGCGTGCCCAAGACACTGGTCTACCGCATCCTGCGCAAGGgtgaggtgag                              >  1:1086146/1‑150 (MQ=255)
                                                                                                             gacgacCAATACCCCTCCGACTTCCGGCGTTCAGCTGATCGAAGTCGCGCCGGAGCTTGCCGGCCAACGCATCGATAATTTCCTCATCACGGCCCTCAAGGGCGTGCCCAAGACACTGGTCTACCGCATCCTGCGCAAGGgtgaggtgag                              <  1:344342/150‑1 (MQ=255)
                                                                                                                     aTACCCCTCCGACTTCCGGCGTTCAGCTGATCTAAGTCGCGCCGGAGCGTGCCGGCAAACGCATCGATAATTTCCTCATCAAGGCCATCAAGGGCGTGCCCAAGACACTGGTCTACCGCATCCTGCGCAAGGGTGAGGTGAGGGTcaaca                      >  2:77717/1‑150 (MQ=255)
                                                                                                                       aCCCCTCCGACTTCCGGCGTTCAGCTGATCGAAGTCGCGCCGGAGCTTGCCGGCCAACGCATCGATAATTTCCTCATCACGGCCCTCAAGGGCGTGCCCAAGACACTGGTCTACCGCATCCTGCGCAAGGGTGAGGTGAGGGTCAACAAg                    <  1:498471/150‑1 (MQ=255)
                                                                                                                           cTCCGACTTCCGGCGTTCAGCTGATCGAAGTCGCGCCGGAGCTTGCCGGCCAACGCATCGATAATTTCCTCATCACGGCCCTCAAGGGCGTGCCCAAGACACTGGTCTACCGCATCCTGCGCAAGGGTGAGGTGAGGGTCAACAAGGGcc                >  2:945482/1‑150 (MQ=255)
                                                                                                                                     cGGCGTTCAGCTGATCGAAGTCGCGCCGGAGCTTGCCGGCCAACGCATCGATAATTTCCTCATCACGGCCCTCAAGGGCGTGCCCAAGACACTGGTCTACCGCATCCTGCGCAAGGGTGAGGTGAGGGCCAACAAGGGCCGCGGCAAGcc      >  1:145801/1‑150 (MQ=255)
                                                                                                                                        cGTTCAGCTGATCGAAGTCGCGCCGGAGCTTGCCGGCCAACGCATCGATAATTTCCTCATCACGGCCCTCAAGGGCGTGCCCAAGACACTGGTCTACCGCATCCTGCGCAAGGGTGAGGTGAGGGTCAACAAGGGCCGCGTCAAGCCGGa   <  2:1086146/150‑1 (MQ=255)
                                                                                                                                         gTTCAGCTGATCGAAGTCGCGCCGGAGCTTGCCGGCCAACGCATCGATAATTTCCTCATCACGGCCCTCAAGGGCGTGCCCAAGACACTGGTCTACCGCATCCTGCGCAAGGCTGAGGTGAGGGTCAACAAGGGCCGCGTCAAGCCGGAg  >  2:941954/1‑150 (MQ=255)
                                                                                                                                              gCTGATTGAAGTCGCGCCGGAGCTTGCCGGCCAACGCATCGATAATTTCCTCATCACGGCCCTCAAGGGCGTGCCCAAGACACTGGTCTACCGCATCCTGCGCAAGggtgag                                   >  2:532228/1‑112 (MQ=255)
                                                                                                                                              gCTGATTGAAGTCGCGCCGGAGCTTGCCGGCCAACGCATCGATAATTTCCTCATCACGGCCCTCAAGGGCGTGCCCAAGACACTGGTCTACCGCATCCTGCGCAAGggtgag                                   >  2:532382/1‑112 (MQ=255)
                                                                                                                                              gCTGATTGAAGTCGCGCCGGAGCTTGCCGGCCAACGCATCGATAATTTCCTCATCACGGCCCTCAAGGGCGTGCCCAAGACACTGGTCTACCGCATCCTGCGCAAGggtgag                                   <  1:532228/112‑1 (MQ=255)
                                                                                                                                              gCTGATTGAAGTCGCGCCGGAGCTTGCCGGCCAACGCATCGATAATTTCCTCATCACGGCCCTCAAGGGCGTGCCCAAGACACTGGTCTACCGCATCCTGCGCAAGggtgag                                   <  1:532382/112‑1 (MQ=255)
                                                                                                                                              gCTGATCGAAGTCGCGCCGGAGCTTGCCGGCCAACGCATCGATAATTTCCTCATCACGGCCCTCAAGGGCGTGCCCAAGACACTGGTCTACCGCATCCTGCGCAAGGGTGAGGTGAGGGTCAACAAGGGCCGCGTCAAGCCGGa   <  2:983/144‑1 (MQ=255)
                                                                                                                                              gCTGATCGAAGTCGCGCCGGAGCTTGCCGGCCAACGCATCGATAATTTCCTCATCACGGCCCTCAAGCGCGTGCCCAAGACACTGGTCTACCGCATCCTGCGCAAGGGTGAGGTGAGGGTCAACAAGGGCCGCGTCAAGCCGGa   >  1:983/1‑144 (MQ=255)
                                                                                                                                                    |                                                                                                                                          
CCGGCAAAGAGCCTGGTGCCGCTGGCCTGCCGCTTTTGGCGGCGTTCGCGACTATAGCAGCAATGATTAAGTGCTTCAATTCCATAAAAAATTGTTATGATCCCGTAATGACGACCAATACCCCTCCGACTTCCGGCGTTCAGCTGATCGAAGTCGCGCCGGAGCTTGCCGGCCAACGCATCGATAATTTCCTCATCACGGCCCTCAAGGGCGTGCCCAAGACACTGGTCTACCGCATCCTGCGCAAGGGTGAGGTGAGGGTCAACAAGGGCCGCGTCAAGCCGGAG  >  NC_002947/2151253‑2151539

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: