Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 172,649 C→T 8.2% intergenic (+3613/+74) PP_RS08750 → / ← fhuA_1 tRNA‑Val/Ferrichrome outer membrane transporter/phagereceptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB122172,6490CT 8.2% 146.0 / 3.8 49intergenic (+3613/+74)PP_RS08750/fhuA_1tRNA‑Val/Ferrichrome outer membrane transporter/phagereceptor
Reads supporting (aligned to +/- strand):  ref base C (25/20);  new base T (2/2);  total (27/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.02e-01

GCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTT  >  ACB122/172524‑172782
                                                                                                                             |                                                                                                                                      
gCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGAtt                                                                                                                               >  5:203398/1‑135 (MQ=39)
 ccAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTc                                                                                                                              >  4:66451/1‑135 (MQ=38)
  cAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATTGTCAAGATTCa                                                                                                                             >  2:38443/1‑135 (MQ=37)
  cAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATTGTCAAGATTCa                                                                                                                             >  2:38439/1‑135 (MQ=37)
  cAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCa                                                                                                                             >  8:91272/1‑135 (MQ=37)
       cacaGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCa                                                                                                                       <  8:12439/136‑1 (MQ=38)
             aTTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCa                                                                                                                       >  1:156420/1‑130 (MQ=38)
             aTTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCa                                                                                                                       <  2:156420/130‑1 (MQ=38)
                 cTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTg                                                                                                              <  1:163599/135‑1 (MQ=255)
                    gacgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGt                                                                                                             <  4:133228/133‑1 (MQ=255)
                         ccATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATg                                                                                                    >  2:120011/1‑136 (MQ=255)
                                    tGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCACCCGTCGGCGGTGCAGACAGGGGc                                                                                         >  5:11261/1‑136 (MQ=255)
                                       aaCCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTAATCGTCAAGATTCAAACCCAAAGCCCCAGTCTGCGATGCAGACAGGGGCttt                                                                                      >  7:6666/1‑136 (MQ=255)
                                       aaCCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCACCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCttt                                                                                      >  5:107864/1‑136 (MQ=255)
                                              tGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCtttt                                                                               >  1:141945/1‑136 (MQ=255)
                                              tGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCtttt                                                                               >  4:17306/1‑136 (MQ=255)
                                              tGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCtttt                                                                               >  1:40663/1‑136 (MQ=255)
                                              tGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCAGTGATGCAAGCACACAGGGGCATTGttgttt                                                                               >  1:187397/1‑131 (MQ=14)
                                              ggagcccatcccGAACTCAGCAGTGAAACGAGGCATCGCCGATGGTCGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCtttt                                                                               <  4:200543/132‑1 (MQ=255)
                                                atcccatcccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATTGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCtttt                                                                               <  1:38439/134‑1 (MQ=255)
                                                atcccatcccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATTGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCtttt                                                                               <  1:38443/134‑1 (MQ=255)
                                                        ccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTc                                                                      >  4:47663/1‑135 (MQ=255)
                                                         cGAACTCAGCAGTGGAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCaa                                                                    >  6:174103/1‑136 (MQ=255)
                                                                   agTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGt                                                             >  7:70713/1‑133 (MQ=255)
                                                                   agTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTACGGGTTTCAATATCGGt                                                             >  5:179958/1‑133 (MQ=255)
                                                                           gATGCATCCCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATa                                                  <  6:107864/136‑1 (MQ=255)
                                                                                   gccttgggTTGTGTGGGGTTGCCCCTTGTGCGATTAGGGCATCGCCAAGATTCAACCCCAACCCCCCTGTCTGCGAGCCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGAc                                            <  7:190491/128‑1 (MQ=255)
                                                                                      gATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAACGCCCCTGTCTGCTATGCAGACTGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTgg                                       <  1:43272/136‑1 (MQ=255)
                                                                                              tgtgtgGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTcaca                               <  4:48457/136‑1 (MQ=255)
                                                                                              tgtgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTcaca                               <  8:189731/136‑1 (MQ=255)
                                                                                              tgtgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTcaca                               <  3:17306/136‑1 (MQ=255)
                                                                                                 ggggTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTg                                                                                                              <  2:184100/55‑1 (MQ=255)
                                                                                                 ggggTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTg                                                                                                              >  1:184100/1‑55 (MQ=255)
                                                                                                 ggggTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTg                            >  5:163986/1‑136 (MQ=255)
                                                                                                        ccccATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGctct                      >  6:158403/1‑135 (MQ=255)
                                                                                                        ccccATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCGGTCTGCGATGCAGACAGGGGGGTTGTGGTTTGCGGGGTTGAATATCGGTTAGAAAGGATAGCTGCCCGTGGCGGTCACATTgggcgatt                     >  8:76376/1‑132 (MQ=25)
                                                                                                          ccATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCa                   >  6:185832/1‑136 (MQ=255)
                                                                                                           cATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCAc                  <  8:42281/136‑1 (MQ=255)
                                                                                                                gagaGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATg             >  2:33549/1‑136 (MQ=255)
                                                                                                                  gagTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTa           <  2:40663/136‑1 (MQ=255)
                                                                                                                   agTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAg          <  3:47663/136‑1 (MQ=255)
                                                                                                                   agTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAg          <  2:141945/136‑1 (MQ=255)
                                                                                                                      aGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAg          >  7:53032/1‑133 (MQ=255)
                                                                                                                      aGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAg          <  8:53032/133‑1 (MQ=255)
                                                                                                                        ggaaTCGTCAAGATTCAAACCCAAAGCCCCTGGCTGCGATCCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGGGGCGCTCACATTGCGCTCTTCACCCATGTAGCagta     <  3:82750/133‑1 (MQ=255)
                                                                                                                        gTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCagta     <  8:70713/136‑1 (MQ=255)
                                                                                                                          cATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTa           >  3:70336/1‑128 (MQ=255)
                                                                                                                          cATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTa           <  4:70336/128‑1 (MQ=255)
                                                                                                                           aTCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGtt  >  4:183001/1‑136 (MQ=255)
                                                                                                                             |                                                                                                                                      
GCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTT  >  ACB122/172524‑172782

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