Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,076,727 (T)5→4 intergenic (+458/‑85) putP → / → BW25113_RS05295 sodium/proline symporter PutP/FTR1 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,076,7230T.89.5% 56.0 / 4.0 19intergenic (+454/‑89)putP/BW25113_RS05295sodium/proline symporter PutP/FTR1 family protein
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base . (7/10);  total (8/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.67e-01

CCCTGGAATAATCGTTTATATCCCTTGGCATTACCTCTCTTTGTTTACATTCCAACATCATTTTATAAACATTCCGCTTGTGTTTTTCTTTGCCGTAATGATAATCGCTATCACTGCGATTTACTTTTCTTTGCATAGATTGACTCAGAAAAA  >  NZ_CP009273/1076641‑1076793
                                                                                  |                                                                      
cccTGGAATAATCGTTTATATCCCTTGGCATTACCTCTCTTTGTTTACATTCCAACATCATTTTATAAACATTCCGCTTGTGttttcttt                                                                 <  1:356278/90‑5 (MQ=255)
  cTGGAATAATCGTTTATATCCCTTGGCATTACCTCTCTTTGTTTACATTCCAACATCATTTTATAAACATTCCGCTTGTGTTTTCTTTGc                                                              >  2:647932/1‑90 (MQ=255)
  cTGGAATAATCGTTTATATCCCTTGGCATTACCTCTCTTTGTTTACATTCCAACATCATTTTATAAACATTCCGCTTGTGTTTTCTTTGc                                                              >  2:542897/1‑90 (MQ=255)
      aataatCGTTTATATCCCTTGGCATTACCTCTCTTTGTTTACATTCCAACATCATTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTa                                                          <  2:520591/90‑1 (MQ=255)
         aatCGTTTATATCCCTTGGCATTACCTCTCTTTGTTTACATTCCAACATCATTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATg                                                       >  1:397488/1‑90 (MQ=255)
               ttATATCCCTTGGCATTACCTCCCTTTGTTTACATTCCAACATCATTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGATAATc                                                 >  2:249393/1‑90 (MQ=255)
                         tGGCATTACCTCTCTTTGTTTACATTCCAACATCATTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGATAATCGCTATCACTg                                       >  2:378843/1‑90 (MQ=255)
                          ggCATTACCTCTCTTTGTTTACATTCCAACATCATTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGATAATCGCTATCACTGc                                      <  1:156049/90‑1 (MQ=255)
                            cATTACCTCTCTTTGTTTACATTCCAACATCATTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGATAATCGCTATCACTGCGa                                    <  2:387291/90‑1 (MQ=255)
                                       tttGTTTACATTCCAACATCATTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGATAATCGCTATCACTGCGATTTACTTTTCt                         <  2:512887/90‑1 (MQ=255)
                                        ttGTTTACATTCCAACATCATTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGATAATCGCTATCACTGCGATTTACtt                             >  2:79206/1‑85 (MQ=255)
                                        ttGTTTACATTCCAACATCATTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGATAATCGCTATCACTGCGATTTACtt                             <  1:79206/85‑1 (MQ=255)
                                                 ttCCAACATCATTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGATAATCGCTATCACTGCGATTTACTTTTCTTTGCATAGAt               <  1:378843/90‑1 (MQ=255)
                                                      aCATCATTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGATAATCGCTATCACTGCGATTTACTTTTCTTTGCATAGATTGACt          <  2:132560/90‑1 (MQ=255)
                                                           atTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGATAATCGCTATCACTGCGATTTACtt                             <  2:45999/66‑1 (MQ=255)
                                                           atTTTATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGATAATCGCTATCACTGCGATTTACtt                             >  1:45999/1‑66 (MQ=255)
                                                              ttATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGATAATCGCTGTCACTGCGATTTACTTTTCTTTGCATAGATTGACTCAGaaaaa  <  2:323541/90‑1 (MQ=255)
                                                              ttATAAACATTCCGCTTGTGTTTTCTTTGCCGTAATGATAATCGCTATCACTGCGATTTACTTTTCTTTGCATAGATTGACTCAGaaaaa  <  2:127715/90‑1 (MQ=255)
                                                                ataAACATTCCGCTTGTGTTTTTCTTTGCCGTAATGATAATCGCTATCACTGCGATTTACTTTTCtt                        <  2:185776/67‑1 (MQ=255)
                                                                ataAACATTCCGCTTGTGTTTTTCTTTGCCGTAATGATAATCGCTATCACTGCGATTTACTTTTCtt                        >  1:185776/1‑67 (MQ=255)
                                                                                  |                                                                      
CCCTGGAATAATCGTTTATATCCCTTGGCATTACCTCTCTTTGTTTACATTCCAACATCATTTTATAAACATTCCGCTTGTGTTTTTCTTTGCCGTAATGATAATCGCTATCACTGCGATTTACTTTTCTTTGCATAGATTGACTCAGAAAAA  >  NZ_CP009273/1076641‑1076793

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