Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 3,432,092 C→G 17.8% intergenic (+51/‑44) CJ019_03257 → / → gpFI hypothetical protein/Putative prophage major tail sheath protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1223,432,0920CG17.8% 105.6 / 2.6 52intergenic (+51/‑44)CJ019_03257/gpFIhypothetical protein/Putative prophage major tail sheath protein
Reads supporting (aligned to +/- strand):  ref base C (10/31);  new base G (7/2);  total (19/33)
Fisher's exact test for biased strand distribution p-value = 4.43e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTCATGACGAGGCCGAAGCTATTGCCGACGAGGCCATTGCCGTGATCAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTT‑TCGTTT‑CTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCAGCATTGCGCTGCCGTCCT  >  ACB122/3431961‑3432208
                                                                                                                                     |                                                                                                                    
aCTCATGACGAGGCCGAAGCTATTGCCGACGAGGCCATTGCCGTGATCAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCACCCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTT‑TCGTTT‑GTGc                                                                                                                   >  1:7550/1‑135 (MQ=255)
              cGAAGCTATTGCCGACGAGGCCATTGCCGTGATCAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTT‑TCGTTT‑CTGCGGGGCCGCCCcgcg                                                                                                     >  7:71527/1‑135 (MQ=255)
               gAAGCTATTGCCGACGAGGCCATTGCCGTGATCAACGCACGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTGT‑TCGGTT‑CGGCGGGGCCGCCCcccgc                                                                                                    >  8:95542/1‑135 (MQ=255)
               gAAGCTATTGCCGACGAGGCCATTGCCGTGATCAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCACCCCGCCTTCCTTGAGCGCCCCGTCCAGCGGGGCGGGG‑TGGGTT‑GTGGGGGGGCGCCCcccgc                                                                                                    >  7:414694/1‑135 (MQ=255)
               gAAGCTATTGCCGACGAGGCCATTGCCGTGATCAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTGTAGTTT‑TGGCGGGGCccccccccc                                                                                                     >  1:180577/1‑132 (MQ=255)
               gAAGCTATTGCCGACGAGGCCATTGCCGTGATCAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTGT‑TGGGTT‑GGGCGGGGGCGCCCcccgc                                                                                                    >  5:364443/1‑135 (MQ=255)
               gAAGCTATTGCCGACGAGGCCATTGCCGTGATCAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGGGT‑TGGGTT‑GTGGGGGGCCCCCCcccgc                                                                                                    >  5:464180/1‑135 (MQ=255)
               gAAGCTATTGACGACGAGGCAATTGACGTGATCAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCCACCCGCCTTCCTTGAGCGCCCCGTCAATCGGGGCGTTT‑TGCGTT‑CTGCGGGGCCGCCCcccgc                                                                                                    >  8:314061/1‑135 (MQ=255)
                aaGCTATTGCCGACGAGGCCATTGCCGTGATCAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTT‑TCGTTT‑CTGCGGGGCCGCCCcgcgcg                                                                                                   <  3:151900/135‑1 (MQ=255)
                 aGCTATTGCCGACGAGGCCATTGCCGTGATCCACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGACCGCCCCGGCCATCGGGGCGTTT‑TCGTTT‑CTGCGGGGCCGCCCCGCGCgg                                                                                                  <  7:218178/135‑1 (MQ=255)
                 aGCTATTGCCGACGAGGCCATTGCCGTGATCAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTGT‑TCGGTT‑GTGCGGGGCCGCCCCGCGCgg                                                                                                  >  3:280647/1‑135 (MQ=255)
                 aGCTATTGCCGACGAGGCCATTGCCGTGATCAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCAGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTT‑TCTTTT‑CTGCGGGGCCGCCCCGCGCgg                                                                                                  <  8:101485/135‑1 (MQ=255)
                  gCTATTGCCGACGAGGCCATTGCCGTGATCAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCGTTCCTTGAGCGCCCCGTCCATCGGGGCGTTT‑TCGTTT‑CTGCGGGGCCGCCCCGCGCgg                                                                                                  <  1:317365/134‑1 (MQ=255)
                   cTATTGCCGCCGAGGCCATTGCCGTGATCAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTT‑TCGTTT‑CTGCGGGGCCGCCCCGCGCGGcc                                                                                                <  4:162116/135‑1 (MQ=255)
                   cTATTGCCGACGAGGCCATTGCCGTGATCAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTT‑TGGGTT‑GTGGGGGGCCGGCCCGCGCGGcc                                                                                                >  1:233308/1‑135 (MQ=255)
                   cTATTGCCGACGAGGCCATTGCCGTGATCAACGCAAGCGTGGTCGGTGTCGGAAACGCCTGACCCGCAACCCGCCTTCCTTGCGCGCCCCGTCCATCGGGGGGTTT‑TCGTTT‑CTGCGGGGCCGCCCCGCGCGGcc                                                                                                <  8:25910/135‑1 (MQ=255)
                   cTATTCCCGACGAGGCCATTGCCGTGATCACCGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAAGCCGCCTTCCTTGAGCGCCCCGGCCATCGGGGCGTTT‑TCGTTT‑CTGCGGGGCCGCCCCGCGCGGcc                                                                                                <  8:453139/135‑1 (MQ=255)
                     tttGCCGACGAGGCCATTGCCGTGATCTACGCATGCGGGGTCGGTGTCGGCAACGCCTGACCGGCAACCCGGGTTCCTTGAGCGCCCCGTGCTGCGGGGCGTTT‑TCGTTT‑GTGCGGGGCGGCCCCGCGCGGCCCt                                                                                              <  7:38478/134‑1 (MQ=255)
                       tGCCGACGAGGCCATTGCCGTGATCAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTT‑TCGTTT‑CTGCGGGGCCGCCCCGCGCGGCCCt                                                                                              <  7:93862/133‑1 (MQ=255)
                       tGCCGACGAGGCCATTGCCGTGATCAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCGTCCTTGAGCGCCCCGTCCATCGGGGCGTTT‑TCGTTT‑CTGCGGGACCGCCCCGCGCGGCCCt                                                                                              <  6:451388/133‑1 (MQ=255)
                         ccGACGAGGCCATTGCCGTGATCAACGCAAGCGTGTGCGGTGTCGGCAACGCCTGACCCGCAACCAGCCTTCCTTGAGCGCCCCGTCCATCGGGGAGTTT‑TCGTTT‑CTGCGGGGCCGCCCCGCGCGGCCCTTGt                                                                                           <  6:74705/134‑1 (MQ=255)
                         ccGACGAGGCCATTGCCGTGATAAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGCCCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCAGCGGGGCGTTT‑TGGGTT‑AGTGGGGGGCGCCCCCCGCGGgccggggt                                                                                          >  3:282230/1‑130 (MQ=255)
                               aGGCCATTGCCGTGATCAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTT‑TCGTTT‑CTTCGGGGCCGCCCCGCGCGGCCCGTGTCTTAtgt                                                                                    >  8:135143/1‑133 (MQ=255)
                               aGGCCATTGCCGTGATCAACGCAAGCGGGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTT‑TCGTTT‑CTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCt                                                                                    <  5:233489/135‑1 (MQ=255)
                                  ccATTGCCGTGATCAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTT‑TCGTTT‑CTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGa                                                                                 <  6:318679/135‑1 (MQ=255)
                                              tCAACGCATGCGTGGTCGGGGGCGGAAACGCCGGACCCGCAACCCGCCTTCCTTGAGCGCCCCTTCCATCGGGGCGTTT‑TCGTTT‑CTGCGGGGCCGCCCCGCGCGGCCCTGGTCTTATCTGGAGCCCACATGACt                                                                     <  6:149536/135‑1 (MQ=255)
                                              tCAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTGT‑TGGTTT‑CTGCGGGGCCGCCCCGCGCGGCCCGTGTCTTTTCTGGGGCCCACATGAcg                                                                     >  8:406512/1‑134 (MQ=255)
                                                 aCGCAAGCGGGGTGGGTGTGGGAAACGGGGTACCCGCAACCCGCCTTCCTTCCGCGCCCCGTCCTGCGGGG‑GTTT‑TTGTTTGCTGCGGGGCCGCCCGGCGCGGCCGTTGTCTTATCTGGAGCCCACATGACTgg                                                                   <  8:170168/134‑1 (MQ=255)
                                                     aaGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTT‑TCGTTT‑CTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGa                                                                                 <  8:63171/116‑1 (MQ=255)
                                                     aaGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTT‑TCGTTT‑CTGCGGGGCCGCCCCGCGCGGCCCGTGTCTTATCTGGa                                                                                 >  7:63171/1‑116 (MQ=255)
                                                     aaGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTGT‑TGGTTT‑CTGCGGGGCCGCCCCGCGCCGCCCGGGTCGTATGTGGAGGCCACAATACGGGGTTCt                                                              >  7:202145/1‑135 (MQ=255)
                                                      cgCGTGGTGGGTGTCGGCAAGGCCTGACCCGAAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTT‑TCGTTT‑CTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCtt                                                             <  8:48915/134‑1 (MQ=255)
                                                      aGCGTGGGCGGTGGCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTT‑TCGTTT‑CTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCtt                                                             <  7:441855/135‑1 (MQ=255)
                                                           gggcGGGGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGCGCGCCCCGGCCATCGGGGCGTTT‑TCGTTG‑CTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACg                                                        <  1:293745/132‑1 (MQ=255)
                                                               ggTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTT‑TCGTTT‑CTGCGGGGCCGCCCCGCGCGGCCCTGGGCTTATCTGGGGCCCCCATGTCTGGATTCTTTCACGGCGt                                                    >  4:314924/1‑135 (MQ=255)
                                                                    cGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTT‑TGGGTT‑CTGCGGGGCCGCCCCGCGCGGCCCGTGTCGGATCGGGGGGCCCCATGGCGGGGTTGTTTGACGGCGGTGCCg                                               >  5:175491/1‑135 (MQ=255)
                                                                         aCGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTT‑TCGTTT‑GTGCGGGGCCGCCCCGCGCGGGCCTTGTGTTTTATGGCGCCCACATGACTGGATTGTTTAACGGCGTGACCGGaac                                           >  4:226030/1‑134 (MQ=255)
                                                                           gCCTGACCCGAAACCCGCCTTCCTTGAGCGCCCGGTCCTTGGGGGCGTTT‑GCTTTT‑CTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTaacgaa                                        <  6:464180/135‑1 (MQ=255)
                                                                                        cccGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTT‑TGGTTT‑CTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGa                                                                       >  7:243418/1‑91 (MQ=255)
                                                                                        cccGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTT‑TCGTTT‑CTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGa                                                                       <  8:243418/91‑1 (MQ=255)
                                                                                         ccGCCTTCCTTGCGCGCCCCGGCCATCGGGGCGTTT‑TCGTTT‑CGGCGGGGCCGCCCCGCGCGGCCCTGGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCgggg                          <  2:414960/135‑1 (MQ=255)
                                                                                          cGCCTTCCTTGAGCGCCCCGTCCATCGGGGGGTTT‑TCGGTT‑CTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGgc                         <  8:452324/135‑1 (MQ=255)
                                                                                          cGCCTTCCTTGAGCGCCCCGGCCATCGGGGAGTTT‑TCGTTG‑CTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGgc                         <  4:380028/135‑1 (MQ=255)
                                                                                                 cttGCGCGCCCGGTCCCGCGGGGCGTTT‑TCGTTG‑CGGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCgca                    <  2:265590/133‑1 (MQ=255)
                                                                                                 cttGAGCGCCCCGTCCCGCGGGGCGTTT‑TCGTTT‑CTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCgca                    <  1:30835/133‑1 (MQ=255)
                                                                                                 cttGAGCGCCCCGTCCATCGGGGCGTTT‑TCGTTT‑CTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCgca                    <  8:71527/133‑1 (MQ=255)
                                                                                                  ttccgcgcCCCGCCCAGGGGGGCGTTT‑GCGTTG‑CGGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCgcagca                 <  1:447690/131‑1 (MQ=255)
                                                                                                       cgcCCCGTCCATCGGGGCGTTT‑TCGTTT‑CTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCAGCATTgcg            <  8:76058/135‑1 (MQ=255)
                                                                                                          cccGTCCATCGGGGCGTTT‑TCGTTT‑CTGCGGGGCCGCCCCGCGCGGCCCGTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCAGCATTGCGCTg         <  6:220925/135‑1 (MQ=255)
                                                                                                                 ctCGGGGCGTTT‑TCGTGT‑CTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCAGCATTGCGCTGCCGtcct  <  2:366379/134‑1 (MQ=255)
                                                                                                                 aTCGGGGCGTTT‑TCGTTT‑CTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCAGCATTGCGCTGCCGtcct  <  8:363709/135‑1 (MQ=255)
                                                                                                                 aTCGCGGGGTTT‑TCGATT‑GTGCGGGGCCGCCCCGCTCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCAGCATTGCGCTGCCGtcct  <  8:71823/135‑1 (MQ=255)
                                                                                                                                     |                                                                                                                    
ACTCATGACGAGGCCGAAGCTATTGCCGACGAGGCCATTGCCGTGATCAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTT‑TCGTTT‑CTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCAGCATTGCGCTGCCGTCCT  >  ACB122/3431961‑3432208

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