Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 226,300 T→C 5.1% noncoding (542/2904 nt) rrlH → 23S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP009048226,3000TC 5.1% 392.5 / 7.4 137noncoding (542/2904 nt)rrlH23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (68/62);  new base C (3/4);  total (71/66)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 143 total aligned reads displayed.

GGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCG  >  AP009048/226231‑226370
                                                                     |                                                                      
ggCGAAAAGAACCCCGGCCAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATg                                                                       <  1:2069943/71‑1 (MQ=17)
 gCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGc                                                                      <  1:5805552/71‑1 (MQ=34)
  cGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCt                                                                     <  1:6343241/71‑1 (MQ=34)
   gAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCt                                                                     >  1:1933621/1‑70 (MQ=35)
     aaaGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAg                                                                  >  1:7896166/1‑71 (MQ=34)
      aaGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCt                                                                     >  1:4174558/1‑67 (MQ=35)
      aaGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAgg                                                                 >  1:2621903/1‑71 (MQ=34)
           ccccGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCgtg                                                             <  1:363367/70‑1 (MQ=35)
            cccGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCAt                                                                        <  1:5004856/58‑1 (MQ=31)
               ggCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACt                                                        >  1:6392815/1‑71 (MQ=34)
                 cGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGc                                                      >  1:7715246/1‑71 (MQ=34)
                      ggAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTAcc                                                 <  1:7795138/71‑1 (MQ=34)
                       gAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCt                                                >  1:234955/1‑71 (MQ=34)
                       gAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCt                                                >  1:7462528/1‑71 (MQ=34)
                        aGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCtt                                               <  1:664744/71‑1 (MQ=34)
                         gTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCttt                                              >  1:3389556/1‑71 (MQ=34)
                          tGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCtttt                                             >  1:1016332/1‑71 (MQ=34)
                          tGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCtt                                               <  1:4857287/69‑1 (MQ=35)
                          tGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCtt                                               <  1:1040391/69‑1 (MQ=35)
                          tGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCG                                                               >  1:6714369‑M2/1‑53 (MQ=255)
                            aaaaaGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTa                                                                   >  1:1010816/1‑47 (MQ=38)
                            aaaaaGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTa                                                                   >  1:5203070/1‑47 (MQ=38)
                            aaaaaGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTg                                            >  1:6551940/1‑70 (MQ=35)
                             aaaaGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAg                                                                  <  1:2317061/47‑1 (MQ=38)
                             aaaaGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAg                                                                  <  1:5585895/47‑1 (MQ=38)
                             aaaaGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAg                                                                  <  1:1815121/47‑1 (MQ=38)
                              aaaGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGta                                          >  1:2795673/1‑70 (MQ=35)
                              aaaGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGta                                          >  1:6989932/1‑70 (MQ=35)
                               aaGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGa                                                          <  1:4311353/53‑1 (MQ=35)
                               aaGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGa                                                          >  1:7526771/1‑53 (MQ=35)
                               aaGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGC                                                                <  1:2074707‑M2/47‑1 (MQ=255)
                                aGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTa                                                                   >  1:4083065/1‑43 (MQ=38)
                                aGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGa                                                          >  1:4533717/1‑52 (MQ=35)
                                aGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGa                                                          >  1:3382218/1‑52 (MQ=35)
                                aGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATaa                                       >  1:3864628/1‑71 (MQ=34)
                                aGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATaa                                       >  1:6364523/1‑71 (MQ=34)
                                aGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATaa                                       <  1:6998833/71‑1 (MQ=34)
                                   aCCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCt                                                                     >  1:4978433/1‑38 (MQ=40)
                                   aCCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATgg                                    >  1:2271276/1‑71 (MQ=34)
                                   aCCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACATTTTGta                                          >  1:5708179/1‑65 (MQ=25)
                                    ccTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTa                                                                   <  1:3760066/39‑1 (MQ=40)
                                    ccTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTa                                                                   <  1:3088307/39‑1 (MQ=40)
                                     cTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGC                                                                >  1:528305‑M2/1‑41 (MQ=255)
                                     cTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGC                                                                >  1:4549344‑M2/1‑41 (MQ=255)
                                     cTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGC                                                                >  1:7152365‑M2/1‑41 (MQ=255)
                                     cTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGC                                                                >  1:844986‑M2/1‑41 (MQ=255)
                                     cTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGC                                                                >  1:1949896‑M2/1‑41 (MQ=255)
                                      tGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACt                                                        <  1:2829785/48‑1 (MQ=38)
                                      tGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTg                                            >  1:5018490/1‑60 (MQ=35)
                                      tGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATggg                                   >  1:6809811/1‑69 (MQ=35)
                                        aaaCCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCttt                                              >  1:4179325/1‑56 (MQ=35)
                                        aaaCCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTc                                 <  1:513008/69‑1 (MQ=35)
                                         aaCCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGc                              <  1:1310610/71‑1 (MQ=34)
                                         aaCCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGc                              <  1:2629430/71‑1 (MQ=34)
                                         aaCCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGc                              <  1:3635753/71‑1 (MQ=34)
                                         aaCCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGc                              >  1:740508/1‑71 (MQ=34)
                                           ccGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTc                                 <  1:1763792/66‑1 (MQ=35)
                                           ccGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTc                                 <  1:4107163/66‑1 (MQ=35)
                                               gtacgtacAAGCAGTGGGAGCACGCTTAGGCGTacacggtttcaggttctttttcactcccctcgccgg                          >  1:5017502‑M2/1‑33 (MQ=255)
                                               gtacgtacAAGCAGTGGGAGCACGCTTAGGCGTacacggtttcaggttctttttcactccactcgcc                            <  1:409489‑M2/67‑35 (MQ=255)
                                                 acgtacAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTtata                      >  1:2151637/1‑71 (MQ=34)
                                                  cgtacAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTtata                      <  1:3163905/70‑1 (MQ=35)
                                                  cgtacAAGCAGTGGGAGCATGCTTAGCCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTtata                      <  1:794234/70‑1 (MQ=25)
                                                   gtacAAGCAGTGGGAGCATGCTTAGGCGTacacggtttcaggttctttttcactcccctcgccggggttct                    <  1:7452879‑M2/71‑43 (MQ=255)
                                                   gtacAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCttt                                              >  1:5324795/1‑45 (MQ=38)
                                                   gtacAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATAtt                    >  1:1106597/1‑71 (MQ=34)
                                                   gtacAAGCAGTGGGAGCACGCTTAGGCGTacacggtttcaggttctttttcactcccctcgcc                            >  1:4352168‑M2/1‑29 (MQ=255)
                                                    tacAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAg                               <  1:7358497/59‑1 (MQ=35)
                                                    tacAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAg                               <  1:7852256/59‑1 (MQ=35)
                                                    tacAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAg                               <  1:6529289/59‑1 (MQ=35)
                                                     acAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGta                                          >  1:6609759/1‑47 (MQ=38)
                                                     acAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCt                  >  1:425481/1‑71 (MQ=34)
                                                     acAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCt                  >  1:434832/1‑71 (MQ=34)
                                                     acAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCt                  >  1:2545282/1‑71 (MQ=34)
                                                      cAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCt                  <  1:661792/70‑1 (MQ=35)
                                                      cAAGCAGTGGGAGCACGCTTGTACGTacacggtttcaggttctttttcactcccctcgccggggttc                     <  1:2394356‑M2/67‑42 (MQ=255)
                                                       aaGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTg                 >  1:748124/1‑70 (MQ=35)
                                                       aaGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTg                 >  1:6346096/1‑70 (MQ=35)
                                                       aaGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGt                <  1:3492039/71‑1 (MQ=34)
                                                         gCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATaa                                       >  1:696232/1‑46 (MQ=38)
                                                         gCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAg              >  1:566552/1‑71 (MQ=34)
                                                           aGTGGGAGCATGCTTAGGCGTacacggtttcaggttctttttcactcccctcgccggg                         <  1:3922802‑M2/58‑38 (MQ=255)
                                                           aGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCa            <  1:6707024/71‑1 (MQ=34)
                                                             tGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCaa           <  1:6133785/70‑1 (MQ=35)
                                                             tGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAg          >  1:3568052/1‑71 (MQ=34)
                                                               ggAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTta                        <  1:4392620/55‑1 (MQ=35)
                                                               ggAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGt        >  1:5912354/1‑71 (MQ=34)
                                                               ggAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGt        >  1:1472728/1‑71 (MQ=34)
                                                                 aGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTa      <  1:2159529/71‑1 (MQ=34)
                                                                  gCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTaa     <  1:2464687/71‑1 (MQ=34)
                                                                   cATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGc                              <  1:3777794/45‑1 (MQ=38)
                                                                   cATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGc                              <  1:1725804/45‑1 (MQ=38)
                                                                   cATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCt                  <  1:2376826/57‑1 (MQ=35)
                                                                    aTGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTtata                      <  1:5961458/52‑1 (MQ=33)
                                                                    aTGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTtata                      <  1:3326394/52‑1 (MQ=33)
                                                                     tGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGc                              <  1:764172/43‑1 (MQ=32)
                                                                     tGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGc                              <  1:1464785/43‑1 (MQ=32)
                                                                     tGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAAcc   >  1:3407065/1‑70 (MQ=31)
                                                                     tGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCg  <  1:2513160/71‑1 (MQ=31)
                                                                     tGCTTAGGCCTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAAcc   >  1:6518601/1‑70 (MQ=14)
                                                                     |                                                                      
GGCGAAAAGAACCCCGGCGAGGGGAGTGAAAAAGAACCTGAAACCGTGTACGTACAAGCAGTGGGAGCATGCTTAGGCGTGTGACTGCGTACCTTTTGTATAATGGGTCAGCGACTTATATTCTGTAGCAAGGTTAACCG  >  AP009048/226231‑226370

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: