Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,256,841 G→T I165I (ATC→ATA psuT ← putative pseudouridine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,256,8410GT96.0% 163.2 / ‑0.4 50I165I (ATC→ATApsuTputative pseudouridine transporter
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base T (30/18);  total (32/18)
Fisher's exact test for biased strand distribution p-value = 5.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.87e-01

GGATCGCCATCAGCGATGCCGCTAACAGGTAGTCAATTGGTACGCCCATTCCGGCATAACCAATCATCATCGAACCAGCAATGGACGCCATCCCGCTACAAATTGCGGTAAACAACTCGTTGCGATTCATGCGATCGATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAATAGTAGTAACCGCAACAAAAGATTCGATTTTGCTGATGTTGAGGGCTTTCTGGAAAATGCTGCCAAGGATGCGAATCAGCAGCCCCATCACGC  >  NC_000913/2256705‑2256980
                                                                                                                                        |                                                                                                                                           
ggATCGCCATCAGCGATGCCGCTAACAGGTAGTCAATTGGTACGCCCATTCCGGCATAACCAATCATCATCGAACCAGCAATGGACGCCATCCCGCTACAAATTGCGGTAAACAACTCGTTGCGATTCATGCGATCtata                                                                                                                                          <  2:1117693/140‑1 (MQ=255)
  aTCTCCATCAGCGATGCCGCTAACAGGTAGTCAATTGGTACGCCCATTCCGGCATAACCAATCATCATCGAACCAGCAATGGACGCCATCCCGCTACAAATTGCGGTAAACAACTCGTTGCGATTCATGCGATCTATaaa                                                                                                                                        <  2:24301/140‑1 (MQ=255)
  aTCGCCATCAGCGATGCCGCTAACAGGTAGTCAATTGGTACGCCCATTCCGGCATAACCAATCATCATCGAACCAGCAATGGACGCCATCCCGCTACAAATTGCGGTAAACAACTCGTTGCGATTCATGCGATCTATaaa                                                                                                                                        <  1:100206/140‑1 (MQ=255)
  aTCGCCATCAGCGATGCCGCTAACAGGTAGTCAATTGGTACGCCCATTCCGGCATAACCAATCATCATCGAACCAGCAATGGACGCCATCCCGCTACAAATTGCGGTAAACAACTCGTTGCGATTCATGCGATCTATaaa                                                                                                                                        <  2:1044968/140‑1 (MQ=255)
  aTCGCCATCAGCGATGCCGCTAACAGGTAGTCAATTGGTACGCCCATTCCGGCATAACCAATCATCATCGAACCAGCAATGGACGCCATCCCGCTACAAATTGCGGTAAACAACTCGTTGCGATTCATGCGATCTATaaa                                                                                                                                        <  1:447838/140‑1 (MQ=255)
     gCCATCAGCGATGCCGCTAACAGGTAGTCAATTGGTACGCCCATTCCGGCATAACCAATCATCATCGAACCAGCAATGGACGCCATCCCGCTACAAATTGCGGTAAACAACTCGTTGCGATTCATGCGATCTATAAACgg                                                                                                                                     >  2:733080/1‑140 (MQ=255)
           aGCGATGCCGCTAACAGGTAGTCAATTGGAACGCCCATTCCGGCATAACCAATCATCATCGAACCAGCAATGGACGCCATCCCGCTACAAATTGCGGTAAACAACTCGTTGCGATTCATGCGATCTATAAACGGTTTAAc                                                                                                                               >  2:1208814/1‑140 (MQ=255)
                                  aaTTGGTACGCCCATTCCGGCATAACCAATCATCATCGAACCAGCAATGGACGCCATCCCGCTACAAATTGCGGTAAACAACTCGTTGCGATTCATGCGATCTATAAACGGTTTAACGATCGCCGGGATCTCATTTTGcc                                                                                                        >  1:136978/1‑140 (MQ=255)
                                   aTTGGTACGCCCATTCCGGCATAACCAATCATCATCGAACCAGCAATGGACGCCATCCCGCTACAAATTGCGGTAAACAACTCGTTGCGATTCATGCGATCTATAAACGGTTTAACGATCGCCGGGATCTCATTTTGccc                                                                                                       <  1:807574/140‑1 (MQ=255)
                                        tACGCCCATTCCGGCATAACCAATCATCATCGAACCAGCAATGGACGCCATCCCGCTACAAATTGCGGTAAACAACTCGTTGCGATTCATGCGATCTATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGa                                                                                                  <  2:422352/140‑1 (MQ=255)
                                          cGCCCATTCCGGCATAACCAATCATCATCGAACCAGCAATGGACGCCATCCCGCTACAAATTGCGGTAAACAACTCGTTGCGATTCATGCGATCTATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGaaa                                                                                                >  2:513867/1‑140 (MQ=255)
                                             ccATTCCGGCATAACCAATCATCATCGAACCAGCAATGGACGCCATCCCGCTACAAATTGCGGTAAACAACTCGTTGCGATTCATGCGATCTATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAAta                                                                                             >  2:99894/1‑140 (MQ=255)
                                              cATTCCGGCATAACCAATCATCATCGAACCAGCAATGGATGCCATCCCGCTACAAATTGCGGTAAACAACTCGTTGCGATTCATGCGATCTATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAAtag                                                                                            >  1:188776/1‑140 (MQ=255)
                                                 tCCGGCATAACCAATCATCATCGAACCAGCAATGGACGCCATCCCGCTACAAATTGCGGTAAACAACTCGTTGCGATTCATGCGATCTATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGa                                                                                                  <  1:1003589/131‑1 (MQ=255)
                                                 tCCGGCATAACCAATCATCATCGAACCAGCAATGGACGCCATCCCGCTACAAATTGCGGTAAACAACTCGTTGCGATTCATGCGATCTATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGa                                                                                                  >  2:1003589/1‑131 (MQ=255)
                                                 tCCGGCATAACCAATCATCATCGAACCAGCAATGGACGCCATCCCGCTACAAATTGCGGTAAACAACTCGTTGCGATTCATGCGATCTATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAAtagtag                                                                                         >  2:205422/1‑140 (MQ=255)
                                                      cATAACCAATCATCATCGAACCAGCAATGGACGCCATCCCGCTACAAATTGCGGTAAACAACTCGTTGCGATTCATGCGATCGATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAATAGTAGTAAcc                                                                                    >  2:3624/1‑140 (MQ=255)
                                                        taaacAATCATCATCGAACCAGCAATGGACGCCATCCCGCTACAAATTGCGGTAAACAACTCGTTGCGATTCATGCGATCTATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAATAGTAGTAACCGc                                                                                  >  2:131516/5‑140 (MQ=255)
                                                               tcatcatcGAACCAGCAATGGACGCCATCCCGCTACAAATTGCGGTAAACAACTCGTTGCGATTCATGCGATCTATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAATAGTAGTAACCGCAACaaaa                                                                           >  2:652695/1‑140 (MQ=255)
                                                                   catcGAACCAGCAATGGACGCCATCCCGCTACAAATTGCGGTAAACAACTCGTTGCGATTCATGCGATCTATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAATAGTAGTAAc                                                                                     <  2:94581/126‑1 (MQ=255)
                                                                   catcGAACCAGCAATGGACGCCATCCCGCTACAAATTGCGGTAAACAACTCGTTGCGATTCATGCGATCTATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAATAGTAGTAAc                                                                                     >  1:94581/1‑126 (MQ=255)
                                                                   catcGAACCAGCAATGGACGCCATCCCGCTACAAATTGCGGTAAACAACTCGTTGCGATTCATGCGATCTATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAATAGTAGTAACCGCAACAAAAGAtt                                                                       >  2:207553/1‑140 (MQ=255)
                                                                    atcGAACCAGCAATGGACGCCATCCCGCTACAAATTGCGGTAAACAACTCGTTGCGATTCATGCGATCTATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAATAGTAGTAACCGCAACAAAAGATTc                                                                      >  2:1162380/1‑140 (MQ=255)
                                                                    atcGAACCAGCAATGGACGCCATCCCGCTACAAATTGCGGTAAACAACTCGTTGCGATTCATGCGATCTATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAATAGTAGTAACCGCAACAAAAGATTc                                                                      >  2:15218/1‑140 (MQ=255)
                                                                    atcGAACCAGCAATGGACGCCATCCCGCTACAAATTGCGGTAAACAACTCGTTGCGATTCATGCGATCTATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAATAGAAGTAACCGCAACAAAAGATTc                                                                      >  2:220170/1‑140 (MQ=255)
                                                                      cGAACCAGCAATGGACGCCATCCCGCTACAAATTGCGGTAAACAACTCGTTGCGATTCATGCGATCTATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAATAGTAGTAACCGCAACAAAAGATTCGa                                                                    <  1:652695/140‑1 (MQ=255)
                                                                      cGAACCAGCAATGGACGCCATCCCGCTACAAATTGCGGTAAACAACTCGTTGCGATTCATGCGATCTATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAATAGTAGTAACCGCAACAAAAGATTCGa                                                                    >  1:1104525/1‑140 (MQ=255)
                                                                        aaCCAGCAATGGACGCCATCCCGCTACAAATTGCGGTAAACAACTCGTTGCGATTCATGCGATCTATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAATAGTAGTAACCGCAACAAAAGATTCGAtt                                                                  >  2:1099575/1‑140 (MQ=255)
                                                                             gCAATGGACGCCATCCCGCTACAAATTGCGGTAAACAACTCGTTGCGATTCATGCGATCTATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAATAGTAGTAACCGCAACAAAAGATTCGATTTTGCt                                                             >  2:376395/1‑140 (MQ=255)
                                                                             gCAATGGACGCCATCCCGCTACAAATTGCGGTAAACAACTCGTTGCGATTCATGCGATCTATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAATAGTAGTAACCGCAACAAAAGATTCGATTTTGCt                                                             <  1:1208814/140‑1 (MQ=255)
                                                                              cAATGGACGCCATCCCGCTACAAATTGCGGTAAACAACTCGTTGCGATTCATGCGATCTATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAATAGTAGTAACCGCAACAAAAGATTCGATTTTGCTg                                                            <  1:99894/140‑1 (MQ=255)
                                                                              cAATGGACGCCATCCCGCTACAAATTGCGGTAAACAACTCGTTGCGATTCATGCGATCTATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAATAGTAGTAACCGCAACAAAAGATTCGATTTTGCTg                                                            <  1:376395/140‑1 (MQ=255)
                                                                                 tGGACGCCATCCCGCTACAAATTGCGGTAAACAACTCGTTGCGATTCATGCGATCTATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAATAGTAGTAACCGCAACAAAAGATTCGATTTTGCTGATg                                                         <  1:682237/140‑1 (MQ=255)
                                                                                   gACGCCATCCCGCTACAAATTGCGGTAAACAACTCGTTGCGATTCATGCGATCTATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAATAGTAGTAACCGCAACAAAAGAtt                                                                       >  1:1037580/1‑124 (MQ=255)
                                                                                   gACGCCATCCCGCTACAAATTGCGGTAAACAACTCGTTGCGATTCATGCGATCTATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAATAGTAGTAACCGCAACAAAAGAtt                                                                       <  2:1037580/124‑1 (MQ=255)
                                                                                      gCCATCCCGCTACAAATTGCGGTAAACAACTCGTTGCGATTCATGCGATCTATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAATAGTAGTAACCGCAACAAAAGATTCGATTTTGCTGATGTTGAg                                                    >  2:579019/1‑140 (MQ=255)
                                                                                        cATCCCGCTACAAATTGCGGAAAACAACTCGTTGCGATTCATGCGATCTATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAATAGTAGTAACCGCAACAAAAGATTCGATTTTGCTGATGTTGa                                                     >  2:657481/1‑137 (MQ=255)
                                                                                        cATCCCGCTACAAATTGCGGAAAACAACTCGTTGCGATTCATGCGATCTATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAATAGTAGTAACCGCAACAAAAGATTCGATTTTGCTGATGTTGa                                                     <  1:657481/137‑1 (MQ=255)
                                                                                          tCCCGCTACAAATTGCGGTAAACAACTCGTTGCGATTCATGCGATCTATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAATAGTAGTAACCGCAACAAAAGATTCGATTTTGCTGATGTTGAGGGCt                                                <  2:1104525/140‑1 (MQ=255)
                                                                                           cccGCTACAAATTGCGGTNAACAACTCGTTNCGATTCATGCGATCTATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAATAGTAGTAACCGCAACAAAAGATTCGATTTTGCTGATGTTGAGGGCtt                                               >  2:728285/1‑140 (MQ=255)
                                                                                             cGCTACAAATTGCGGTAAACAACTCGTTGCGATTCATGCGATCTATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAATAGTAGTAACCGCAACAAAAGATTCGATTTTGCTGATGTTGAGGGCTTTc                                             >  1:1111368/1‑140 (MQ=255)
                                                                                                     aTTGCGGTAAACAACTCGTTGCGATTCATGCGATCTATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAATAGTAGTAACCGCAACAAAAGATTCGATTTTGCTGATGTTGAGGGGTTTCTGGAAAAt                                     <  2:188776/140‑1 (MQ=255)
                                                                                                         cGGTAAACAACTCGTTGCGATTCATGCGATCTATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAATAGTAGTAACCGCAACAAAAGATTCGATTTTGCTGATGTTGAGGGCTTTCTGGAAAAtgctg                                 >  2:884997/1‑140 (MQ=255)
                                                                                                          ggTAAACAACTCGTTGCGATTCATGCGATCTATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAATAGTAGTAACCGCAACAAAAGATTCGATTTTGCTGATGTTGAGGGCTTTCTGGAAAAtgctgc                                >  1:1042040/1‑140 (MQ=255)
                                                                                                            tAAACAACTCGTTGCGATTCATGCGATCTATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAATAGTAGTAACCGCAACAAAAGATTCGATTTTGCTGATGTTGAGGGCTTTCTGGAAAATGCTGCCa                              >  1:1104921/1‑140 (MQ=255)
                                                                                                            tAAACAACTCGTTGCGATTCATGCGATCTATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAATAGTAGTAACCGCAACAAAAGATTCGATTTTGCTGATGTTGAGGGCTTTCTGGAAAATGCTGCCa                              <  1:5898/140‑1 (MQ=255)
                                                                                                               acaacTCGTTGCGATTCATGCGATCGATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAATAGTAGTAACCGCAACAAAAGATTCGATTTTGCTGATGTTGAGGGCTTTCTGGAAAATGCTGCCAAgg                           >  2:813558/1‑140 (MQ=255)
                                                                                                                 aacTCGTTGCGATTCATGCGATCTATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAATAGTAGTAACCGCAACAAAAGATTCGATTTTGCTGATGTTGAGGGCTTTCTGGAAAATGCTGCCAAGGAt                         >  1:631428/1‑140 (MQ=255)
                                                                                                                        tGCGATTCATGCGATCTATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAATAGTAGTAACCGCAACAAAAGATTCGATTTTGCTGATGTTGAGGGCTTTCTGGAAAATGCTGCCAAGGATGCGAATc                  >  1:1177660/1‑140 (MQ=255)
                                                                                                                        tGCGATTCATGCGATCTATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAATAGTAGTAACCGCAACAAAAGATTCGATTTTGCTGATGTTGAGGGCTTTCTGGAAAATGCTGCCAAGGATGCGAATc                  >  1:596484/1‑140 (MQ=255)
                                                                                                                        tGCGATTCATGCGATCTATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAATAGTAGTAACCGCAACAAAAGATTCGATTTTGCTGATGTTGAGGGCTTTCTGGAAAATGCTGCCAAGGATGCGAATc                  >  1:1209169/1‑140 (MQ=255)
                                                                                                                         gCGATTCATGCGATCTATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAATAGTAGTAACCGCAACAAAAGATTCGATTTTGCTGATGTTGAGGGCTTTCTGGAAAATGCTGCCAAGGATGCGAATca                 <  1:579019/140‑1 (MQ=255)
                                                                                                                          cGATTCATGCGATCTATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAATAGTAGTAACCGCAACAAAAGATTCGATTTTGCTGATGTTGAGGGCTTTCTGGAAAATGCTGCCAAGGATGCGAATcag                >  2:259041/1‑140 (MQ=255)
                                                                                                                                   cgatctatAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAATAGTAGTAACCGCAACAAAAGATTCGATTTTGCTGATGTTGAGGGCTTTCTGGAAAATGCTGCCAAGGATGCGAATCAGCAGCCCCAt       <  1:259041/140‑1 (MQ=255)
                                                                                                                                    gatctatAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAATAGTAGTAACCGCAACAAAAGATTCGATTTTGCTGATGTTGAGGGCTTTCTGGAAAATGCTGCCAAGGATGCGAATCAGCAGCCCCATc      >  1:820250/1‑140 (MQ=255)
                                                                                                                                     atctatAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAATAGTAGTAACCGCAACAAAAGATTCGATTTTGCTGATGTTGAGGGCTTTCTGGAAAATGCTGCCAAGGATGCGAATCAGCAGCCCCATCa     <  2:820250/140‑1 (MQ=255)
                                                                                                                                        gatAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAATAGTAGTAACCGCAACAAAAGATTCGATTTTGCTGATGTTGAGGGCTTTCTGGAAAATGCTGCCAAGGATGCGAATCAGCAGCCCCATCACGc  <  1:3624/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                           
GGATCGCCATCAGCGATGCCGCTAACAGGTAGTCAATTGGTACGCCCATTCCGGCATAACCAATCATCATCGAACCAGCAATGGACGCCATCCCGCTACAAATTGCGGTAAACAACTCGTTGCGATTCATGCGATCGATAAACGGTTTAACGATCGCCGGGATCTCATTTTGCCCGAGGAAAATAGTAGTAACCGCAACAAAAGATTCGATTTTGCTGATGTTGAGGGCTTTCTGGAAAATGCTGCCAAGGATGCGAATCAGCAGCCCCATCACGC  >  NC_000913/2256705‑2256980

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: