Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,258,312 T→C D19G (GAC→GGC)  psuG ← pseudouridine‑5'‑phosphate glycosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,258,3120TC100.0% 212.2 / NA 57D19G (GAC→GGC) psuGpseudouridine‑5'‑phosphate glycosidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (28/29);  total (28/29)

GTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACTAACGATATAACGTTGGCAATCGATAAATCGGGGTTATTGGTGTATTCACAGGAGAGCGCCATTGAC  >  NC_000913/2258184‑2258448
                                                                                                                                |                                                                                                                                        
gTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTcc                                                                                                                               <  2:815132/140‑1 (MQ=255)
   cAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCggg                                                                                                                            <  2:1148293/140‑1 (MQ=255)
      cgcCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCgggg                                                                                                                           >  2:488924/1‑138 (MQ=255)
      cgcCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCgggg                                                                                                                           <  1:488924/138‑1 (MQ=255)
        cccTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAAt                                                                                                                       <  2:330991/140‑1 (MQ=255)
         ccTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAAtt                                                                                                                      <  1:432653/140‑1 (MQ=255)
           tGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTg                                                                                                                    >  2:36861/1‑140 (MQ=255)
             tttACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGta                                                                                                                  <  1:152943/140‑1 (MQ=255)
             tttACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGta                                                                                                                  >  2:483381/1‑140 (MQ=255)
                       gTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAgg                                                                                                        >  1:252025/1‑140 (MQ=255)
                       gTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAgg                                                                                                        >  1:564035/1‑140 (MQ=255)
                       gTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAgg                                                                                                        >  1:851529/1‑140 (MQ=255)
                       gTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAgg                                                                                                        >  2:522809/1‑140 (MQ=255)
                       gTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAgg                                                                                                        >  2:1154520/1‑140 (MQ=255)
                       gTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAgg                                                                                                        >  2:572900/1‑140 (MQ=255)
                                  ttcaaTTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTa                                                                                             >  2:1113969/1‑140 (MQ=255)
                                     aaTTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAAtt                                                                                          <  1:36861/140‑1 (MQ=255)
                                     aaTTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAAtt                                                                                          <  2:666891/140‑1 (MQ=255)
                                        tGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAg                                                                                       <  1:1157128/140‑1 (MQ=255)
                                          cGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGAc                                                                                     >  1:937108/1‑140 (MQ=255)
                                            gTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACAt                                                                                   >  2:181479/1‑140 (MQ=255)
                                                gggCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTg                                                                               <  1:910916/140‑1 (MQ=255)
                                                 ggCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGc                                                                              >  1:170457/1‑140 (MQ=255)
                                                  gCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCg                                                                             <  2:170457/140‑1 (MQ=255)
                                                       ttGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTtctc                                                                        >  2:869205/1‑140 (MQ=255)
                                                           gggAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACt                                                                    <  1:572900/140‑1 (MQ=255)
                                                           gggAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACt                                                                    <  1:1098365/140‑1 (MQ=255)
                                                           gggAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACt                                                                    >  1:224705/1‑140 (MQ=255)
                                                            ggAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGGCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACTa                                                                   >  2:532941/1‑140 (MQ=255)
                                                            ggAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACTa                                                                   <  1:35576/140‑1 (MQ=255)
                                                            ggAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACTa                                                                   <  1:1113969/140‑1 (MQ=255)
                                                             gAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACTaa                                                                  <  1:522809/140‑1 (MQ=255)
                                                                    aTCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACTAACGAtata                                                           <  2:224705/140‑1 (MQ=255)
                                                                        cGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACTAACGATATAACGt                                                       >  2:1040616/1‑140 (MQ=255)
                                                                        cGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACTAACGATATAACGt                                                       >  1:18987/1‑140 (MQ=255)
                                                                        cGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGG‑TTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACTa                                                                   >  1:1002898/1‑127 (MQ=255)
                                                                        cGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGG‑TTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACTa                                                                   <  2:1002898/127‑1 (MQ=255)
                                                                          tGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACTAACGATATAACGTTg                                                     >  1:189768/1‑140 (MQ=255)
                                                                           gagaAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACTAACGATATAACGTTgg                                                    >  2:442201/1‑140 (MQ=255)
                                                                                    tGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACTAACGATATAACGTTGGCAATCGATa                                           >  1:1102192/1‑140 (MQ=255)
                                                                                      gTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACTAACGATATAACGTTGGCAATCGATaaa                                         >  1:593152/1‑140 (MQ=255)
                                                                                         gATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACTAACGATATAACGTTGGCAATCGATAAATCg                                      >  1:876345/1‑140 (MQ=255)
                                                                                         gATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACTAACGATATAACGTTGGCAATCGATAAATCg                                      >  1:197228/1‑140 (MQ=255)
                                                                                         gATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACTAACGATATAACGTTGGCAATCGATAAATCg                                      >  1:246404/1‑140 (MQ=255)
                                                                                            tCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACTAACGATATAACGTTGGCAATCGATAAATCgggg                                   <  1:869205/140‑1 (MQ=255)
                                                                                               aGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACTAACGATATAACGTTGGCAATCGATAAATCGGGGTTa                                <  2:252025/140‑1 (MQ=255)
                                                                                               aGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACTAACGATATAACGTTGGCAATCGATAAATCGGGGTTa                                <  1:857030/140‑1 (MQ=255)
                                                                                               aGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACTAACGATATAACGTTGGCAATCGATAAATCGGGGTTa                                <  2:18987/140‑1 (MQ=255)
                                                                                                 cgcCACAACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACTAACGATATAACGTTGGCAATCGATAAATCGGGGTTAtt                              <  2:876345/140‑1 (MQ=255)
                                                                                                    cacaACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACTAACGATATAACGTTGGCAATCGATAAATCGGGGTTATTGgt                           <  1:532941/140‑1 (MQ=255)
                                                                                                     acaACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACTAACGATATAACGTTGGCAATCGATAAATCGGGGTTATTGgtg                          >  1:64011/1‑140 (MQ=255)
                                                                                                      caACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACTAACGATATAACGTTGGCAATCGATAAATCGGGGTTATTGgtgt                         <  1:483381/140‑1 (MQ=255)
                                                                                                      caACCGGTTTTTTGTTTTTTAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACTAACGATATAACGTTGGCAATCGATAAATCGGGGTTATTGgtgt                         <  1:1154520/140‑1 (MQ=255)
                                                                                                                    ttttttAAAGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACTAACGATATAACGTTGGCAATCGATAAATCGGGGTTATTGGTGTATTCACAGGAGAgc           >  2:326825/1‑140 (MQ=255)
                                                                                                                           aaGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACTAACGATATAACGTTGGCAATCGATAAATCGGGGTTATTGGTGTATTCACAGGAGAGCGCCATTg    <  2:94616/140‑1 (MQ=255)
                                                                                                                           aaGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACTAACGATATAACGTTGGCAATCGATAAATCGGGGTTATTGGTGTATTCACAGGAGAGCGCCATTg    <  2:189768/140‑1 (MQ=255)
                                                                                                                           aaGCGCCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACTAACGATATAACGTTGGCAATCGATAAATCGGGGTTATTGGTGTATTCACAGGAGAGCGCCATTg    <  1:1040616/140‑1 (MQ=255)
                                                                                                                             gcgcCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACTAACGATATAACGTTGGCAATCGATAAATCGGGGTTATTGGTGTATTCACAGGAGAGCGCCATTGAc  <  2:937108/140‑1 (MQ=255)
                                                                                                                                |                                                                                                                                        
GTACAGCGCCCTGTTTACGAATAGTTTCTTCAACTTCAATTGCGGTCTGGGCATTTTGTGGGAACGGCATCCCGTGAGAAATAATGGTCGATTCCAGCGCCACAACCGGTTTTTTGTTTTTTAAAGCGTCCTGCACTTCCGGGGAAATTTGTAATAATTCAGGGGAAATTTTTAATTCAGACATTCTGCGTTCTCCACTAACGATATAACGTTGGCAATCGATAAATCGGGGTTATTGGTGTATTCACAGGAGAGCGCCATTGAC  >  NC_000913/2258184‑2258448

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑