Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,250,936 G→A A304V (GCC→GTC)  psuG ← pseudouridine‑5'‑phosphate glycosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,250,9360GA94.3% 169.0 / ‑0.4 53A304V (GCC→GTC) psuGpseudouridine‑5'‑phosphate glycosidase
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (25/25);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCATTATATCCATAATGATTTTCCCTTCATGCCGGTAAACCCGGCGTCAGCGCCAGGTTTTGGTATGCTTGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTCCCTGAGCT  >  NZ_CP009273/2250801‑2251067
                                                                                                                                       |                                                                                                                                   
ctcATTATATCCATAATGATTTTCCCTTCATGCCGGTAAACCCGGCGTCAGCGCCAGGTTTTGGTATGCTTGATGAGTACGGGCGACGGCTTTCTTCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGACaa                                                                                                                                  >  2:216086/1‑139 (MQ=255)
ctcATTATATCCATAATGATTTTCCCTTCATGCCGGTAAACCCGGCGTCAGCGCCAGGTTTTGGTATGCTTGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGACaa                                                                                                                                  >  2:796339/1‑139 (MQ=255)
   aTTATATCCATAATGATTTTCCCTTCATGCCGGTAAACCCGGCGTCAGCGCCAGGTTTTGGTATGCTTGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAAttt                                                                                                                               >  2:1148628/1‑139 (MQ=255)
       tatCCATAATGATTTTCCCTTCATGCCGGTAAACCCGGCGTCAGCGCCAGGTTTTGGTATGCTTGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTtcgct                                                                                                                           <  2:1237706/139‑1 (MQ=255)
       tatCCATAATGATTTTCCCTTCATGCCGGTAAACCCGGCGTCAGCGCCAGGTTTTGGTATGCTTGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGACAATTtcgct                                                                                                                           <  1:414650/139‑1 (MQ=255)
        atCCATAATGATTTTCCCTTCATGCCGGTAAACCCGGCGTCAGCGCCAGGTTTTGGTATGCTTGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGACAATTtcgctc                                                                                                                          >  2:738157/1‑139 (MQ=255)
        atCCATAATGATTTTCCCTTCATGCCGGTAAACCCGGCGTCAGCGCCAGGTTTTGGTATGCTTGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGACAATTtcgctc                                                                                                                          >  1:907016/1‑139 (MQ=255)
             tAATGATTTTCCCTTCATGCCGGTAAACCCGGCGTCAGCGCCAGGTTTTGGTATGCTTGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGACAATTtcgctcgc                                                                                                                        >  2:386394/1‑136 (MQ=255)
              aaTGATTTTCCCTTCATGCCGGTAAACCCGGCGTCAGCGCCAGGTTTTGGTATGCTTGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGACAATTTCGCTCGCCAg                                                                                                                     >  2:810770/1‑138 (MQ=255)
                     ttCCCTTCATGCCGGTAAACCCGGCGTCAGCGCCAGGTTTTGGTATGCTTGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGACAATTTCGCTCGCCAGAATGGCgt                                                                                                             >  1:642776/1‑139 (MQ=255)
                     ttCCCTTCATGCCGGTAAACCCGGCGTCAGCGCCAGGTTTTGGTATGCTTGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGACAATTTCGCTCGCCAGAATGGCgt                                                                                                             >  2:1238165/1‑139 (MQ=255)
                     ttCCCTTCATGCCGGTAAACCCGGCGTCAGCGCCAGGTTTTGGTATGCTTGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGACAATTTCGCTCGCCAGAATGGCgt                                                                                                             >  1:759561/1‑139 (MQ=255)
                     ttCCCTTCATGCCGGTAAACCCGGCGTCAGCGCCAGGTTTTGGTATGCTTGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGACAATTTCGCTCGCCAGAATGGCg                                                                                                              <  2:535513/138‑1 (MQ=255)
                     ttCCCTTCATGCCGGTAAACCCGGCGTCAGCGCCAGGTTTTGGTATGCTTGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGACAATTTCGCTCGCCAGAATGGCg                                                                                                              <  2:339988/138‑1 (MQ=255)
                     ttCCCTTCATGCCGGTAAACCCGGCGTCAGCGCCAGGTTTTGGTATGCTTGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGACAATTTCGCTCGCCAGAATGGCg                                                                                                              <  2:728050/138‑1 (MQ=255)
                      tCCCTTCATGCCGGTAAACCCGGCGTCAGCGCCAGGTTTTGGTATGCTTGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGACAATTTCGCTCGCCAGAATGGCgtt                                                                                                            >  1:276700/1‑139 (MQ=255)
                      tCCCTTCATGCCGGTAAACCCGGCGTCAGCGCCAGGTTTTGGTATGCTTGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGACAATTTCGCTCGCCAGAATGGCg                                                                                                              <  1:452729/137‑1 (MQ=255)
                      tCCCTTCATGCCGGTAAACCCGGCGTCAGCGCCAGGTTTTGGTATGCTTGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGACAATTTCGCTCGCCAGAATGGCg                                                                                                              >  2:452729/1‑137 (MQ=255)
                           tCATGCCGGTAAACCCGGCGTCAGCGCCAGGTTTTGGTATGCTTGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGACAATTTCGCTCGCCAGAATGGCGttgttg                                                                                                        >  2:327484/1‑138 (MQ=255)
                            cATGCCGGTAAACCCGGCGTCAGCGCCAGGTTTTGGTATGCTTGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGACAATTTCGCTCGCCAGAATGGCGTTGTTGaa                                                                                                      >  1:722026/1‑139 (MQ=255)
                             aTGCCGGTAAACCCGGCGTCAGCGCCAGGTTTTGGTATGCTTGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGACAATTTCGCTCGCCAGAATGGCGTTGTTGaa                                                                                                      <  2:15360/138‑1 (MQ=255)
                                  ggTAAACCCGGCGTCAGCGCCAGGTTTTGGTATGCTTGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGACAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAg                                                                                                <  1:783039/139‑1 (MQ=255)
                                  ggTAAACCCGGCGTCAGCGCCAGGGTTTGGAATGCTTGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAg                                                                                                >  1:890335/1‑139 (MQ=255)
                                  ccaaaaCCCGGCGTCAGCGCCAGGTTTTGGTATGCTTGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGACAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAg                                                                                                >  1:790489/4‑139 (MQ=255)
                                                     ccAGGTTTTGGTATGCTTGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGACAATTTCGCTCGCCAGAATGGCGTTGATGAACACCAGCTGGATGTTGGATTTCAgg                                                                             <  1:216086/139‑1 (MQ=255)
                                                         gTTTTGTTATGCTTGATGATTACGGGCGACGGCTTTCTGCCCGTCAAAAATTACCCCGCGAGACGCTGATATTCTTTGACAATTTCGCTCGCCAGAAGGGCGTTGTTGAACCCCAGCTGGATGTTGGATTTCAGGCTGt                                                                         <  2:1181905/139‑1 (MQ=255)
                                                         gTTTTGGTATGCTTGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGACAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGt                                                                         <  1:386394/139‑1 (MQ=255)
                                                         gTTTTGGTATGCTTGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGACAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGt                                                                         >  2:497826/1‑139 (MQ=255)
                                                              ggTATGCTTGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGACAATTTCGCTCGCCAGAATGGCg                                                                                                              >  2:165057/1‑97 (MQ=255)
                                                              ggTATGCTTGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGACAATTTCGCTCGCCAGAATGGCg                                                                                                              <  1:165057/97‑1 (MQ=255)
                                                               gTATGCTTGAGGAGTACGGGCGACGGCTTTCGGCCCGTCAGAAATTAACCCGCGAGACGCGGATATTCTTTGACAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCAccg                                                                    <  2:322227/138‑1 (MQ=255)
                                                                    cTTGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGACAATTTCGCTCGCCAGAATGGCgt                                                                                                             <  2:775686/92‑1 (MQ=255)
                                                                    cTTGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGACAATTTCGCTCGCCAGAATGGCgt                                                                                                             >  1:775686/1‑92 (MQ=255)
                                                                    cTTGATGAGTACGGGCGACGGCTTTCTCCCCGTCAGAAATTACCCCGCGAGACGCTGATATTCTTTGACATTTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGt                                                               <  1:397961/138‑1 (MQ=255)
                                                                               cGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGACAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAAc                                                   >  1:287358/1‑139 (MQ=255)
                                                                                       ggCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGACAATTtcg                                                                                                                             >  2:182608/1‑57 (MQ=255)
                                                                                       ggCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGACAATTtcg                                                                                                                             <  1:182608/57‑1 (MQ=255)
                                                                                       ggCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGACAATTTCGCTCGCCAGAATgg                                                                                                                <  1:374982/70‑1 (MQ=255)
                                                                                       ggCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGACAATTTCGCTCGCCAGAATgg                                                                                                                >  2:374982/1‑70 (MQ=255)
                                                                                       ggCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGACAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGt                                                               >  2:752738/1‑119 (MQ=255)
                                                                                       ggCTTTCTGCCCGTCAAAAATTAACCCGCGAGACGCTGATATTCTTTGACAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGt                                                               <  1:752738/119‑1 (MQ=255)
                                                                                        gCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGACAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCcag                                          <  2:1247527/139‑1 (MQ=255)
                                                                                                cccGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGACAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAAt                                   <  1:497826/138‑1 (MQ=255)
                                                                                                   gTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGACAATTTCGCTCGCCAGAATGGCg                                                                                                              <  1:467798/60‑1 (MQ=255)
                                                                                                   gTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGACAATTTCGCTCGCCAGAATGGCg                                                                                                              >  2:467798/1‑60 (MQ=255)
                                                                                                   gTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGACAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCg                               >  2:730237/1‑139 (MQ=255)
                                                                                                     cAGAAATTACCCCGCGAGACGCTGATATTCTTTGACAATTTCGCTCGCCAGAATGGCGTTGTTGAACCCCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTa                             <  2:1119715/139‑1 (MQ=255)
                                                                                                       gAAATTACCCCGCGAGACGCTGATATTCTTTGACAATTTCGCTCGCCAGAATGGCGTTGTTGAACCCCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTAct                           <  1:1204893/139‑1 (MQ=255)
                                                                                                             aaCCCGCGAGACGCTGATATTCTTTGACAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTActttctt                      >  1:854815/1‑138 (MQ=255)
                                                                                                             aaCCCGCGAGACGCTGATATTCTTTGACAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTActttctt                      >  2:101107/1‑138 (MQ=255)
                                                                                                                   cgAGACGCTGATATTCTTTGACAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAAt               <  2:854815/139‑1 (MQ=255)
                                                                                                                    gagcCGCTGTTATTCTTTGACATTTCCGCTCGCCAAAATGGCGTTGTTGAACACCAGCGGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATg              <  2:642776/139‑1 (MQ=255)
                                                                                                                                 tCTTTGACAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGTGGTTGGATTCCAGGCGGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTCCCTGagct  <  2:424684/138‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                   
CTCATTATATCCATAATGATTTTCCCTTCATGCCGGTAAACCCGGCGTCAGCGCCAGGTTTTGGTATGCTTGATGAGTACGGGCGACGGCTTTCTGCCCGTCAGAAATTAACCCGCGAGACGCTGATATTCTTTGGCAATTTCGCTCGCCAGAATGGCGTTGTTGAACACCAGCTGGATGTTGGATTTCAGGCTGTCACCGCCGGTCAGTTCAGCAACGCGCGCCAGCAGGAATGGCGTACTTTCTTTACCAATGACTCCCTGAGCT  >  NZ_CP009273/2250801‑2251067

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: