Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 579,846 C→A 11.3% intergenic (‑252/‑34) ybcV ← / → ybcW DLP12 prophage; uncharacterized protein/DLP12 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913579,8460CA11.3% 97.8 / 2.3 35intergenic (‑252/‑34)ybcV/ybcWDLP12 prophage; uncharacterized protein/DLP12 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (21/10);  new base A (2/2);  total (23/12)
Fisher's exact test for biased strand distribution p-value = 5.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTTTTCAATTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAACCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAATGAATAAAGAACAATCTGCTGATGATCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAAC  >  NC_000913/579709‑579994
                                                                                                                                         |                                                                                                                                                    
tGTTTTCAATTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAACCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAag                                                                                                                                           >  2:195187/1‑149 (MQ=255)
tGTTTTCAATTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAACCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAag                                                                                                                                           >  1:584176/1‑149 (MQ=255)
 gTTTTCAATTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAACCTTCATACAGAAAGATTGGAAGTAATATTTTAACCGCTAGATGAaga                                                                                                                                          >  1:42613/1‑149 (MQ=255)
  ttttCAATTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAACCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAagag                                                                                                                                         >  1:553636/1‑149 (MQ=255)
                           tAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAACCTTCATACAGAAAGATTTGCAGGAATATTTTAACCGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAATGaa                                                                                                                >  1:358708/1‑149 (MQ=255)
                            aaCAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAACCTTCATACAGAAAGATTTGAAGTAATATTTTAACAGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAATGAAt                                                                                                               >  1:495087/1‑149 (MQ=255)
                            aaCAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAACCTTCATACAGAAAGATTTGAAGTAATATTTTAACAGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAATGAAt                                                                                                               >  1:495807/1‑149 (MQ=255)
                             aCAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAACCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAATGAATa                                                                                                              <  1:749431/149‑1 (MQ=255)
                               aaaGTGCGGTCCTGCTGGCATTCTGGACGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAACCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAATGAATaaa                                                                                                            <  2:246834/149‑1 (MQ=255)
                                              tgGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAACCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAATGAATAAAGAACAATCTGCtgat                                                                                             >  2:659371/1‑149 (MQ=255)
                                                  aTTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAACCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAATGAATAAAGAACAATCTGCTGATGATc                                                                                         >  1:2810/1‑149 (MQ=255)
                                                  aTTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAACCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAATGAATAAAGAACAATCTGCTGATGATc                                                                                         >  1:583696/1‑149 (MQ=255)
                                                       ggAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAACCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAATGAATAAAGCACAATCTGCTGATGATACCTcc                                                                                    >  2:408992/1‑149 (MQ=255)
                                                          ggtaaaTACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAACCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAATGAATAAAGAACAATCTGCTGATGATCCCTCCGTg                                                                                 <  2:114453/146‑1 (MQ=255)
                                                             aaaTACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAACCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAATGAATGAAGAACAATCTGCTGATGATCCCTCCGTGGAt                                                                              >  1:115682/1‑149 (MQ=255)
                                                             aaaTACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAACCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAATGAATAAAGAACAATCTGCTGATGAGCCCTCCGTGGAt                                                                              >  2:542671/1‑149 (MQ=255)
                                                               aTACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAACCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGCAAGCGAATGACCTGTGTTAAGTAATAAAGAACAATCTGCTGATGATCCCTCCGTGGATCt                                                                            >  1:619761/1‑149 (MQ=255)
                                                                    aCCGACAGATGTATGTAAGGCCAACGTGCTCAAACCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAATGAATAAAGAACAATCTGCTGATGATCCCTCCGTGGATCTGATTc                                                                       <  1:694339/149‑1 (MQ=255)
                                                                    aCCGACAGATGTATGTAAGGCCAACGTGCTCAAACCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAATGAATAAAGAACAATCTGCTGATGATCCCTCCGTGGAGCTGAGTc                                                                       >  1:636300/1‑149 (MQ=255)
                                                                    aCCGACAGATGTATGTAAGGCCAACGTGCTCAAACCTTCATACAGAAAGATTTGAAGTAATATTTTAACAGCTAGATGAAGAGCAAGCGCATGGTGCTACAAAATGAATAAAGAACAATCTGCTGATGATCCCTCCGTGGATCTGATTc                                                                       <  2:495807/149‑1 (MQ=255)
                                                                    aCCGACAGATGTATGTAAGGCCAACGTGCTCAAACCTTCATACAGAAAGATTTGAAGTAATATTTTAACAGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAATGAATAAAGAACAATCTGCTGATGATCCCTCCGTGGATCTGATTc                                                                       <  2:495087/149‑1 (MQ=255)
                                                                    aCCGACAGATGTATGTAAGGCCAACGTGCTCAAACCGTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAATGAATAAAGAACAATCTGCTGATGATCCCTCCGTGGATCTGATTc                                                                       <  2:584176/149‑1 (MQ=255)
                                                                                     aGGCCAACGTGCTCAAACCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAATGAATAAAGAACAATCTGCTGATGATCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCtt                                                      <  1:528600/149‑1 (MQ=255)
                                                                                     aGGCCAACGTGCTCAAACCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCAAGATGAAGAGCAAGCGCATGGAGCGACAAAATGAATAAAGAACAATCTGCTGATGATCCCTCCGTGGATCTGATTCGTGTAAAAAATCTGCtt                                                      >  1:55361/1‑149 (MQ=255)
                                                                                            cGTGCTCAAACCTTCATACAGAAAGATTTGAAGTACTATTTTAACCGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAATGAATAAAGAACAATCTGCTGCTGATCCCTCCGTGGATCTGATTCCTGTAAAAAATATGCTTAATAGCa                                               <  1:420186/149‑1 (MQ=255)
                                                                                                cTCAAACCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAATGAATAAAGAACAATCTGCTGATGATCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATACCACCAt                                           >  2:484784/1‑149 (MQ=255)
                                                                                                cTCAAACCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAATGAATAAAGAACAATCTGCTGATGATCCCTACGTGGATCGGATTCGTGTAAAAAATATGCTTAATAGCACCAt                                           >  2:576430/1‑149 (MQ=255)
                                                                                                 tCAAACCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAATGAATAAAGAACAATCTGCTGATGATCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCAtt                                          >  1:456146/1‑149 (MQ=255)
                                                                                                  cAAACCTTCATACAGAAAGATTTGAAGTAATATGTTAACCGCTAGATGCAGGGCAAGCGCATGGAGCGACAAAATGAATAAAGAACAATCTGCTGATGATCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCAttt                                         <  2:42613/149‑1 (MQ=255)
                                                                                                   aaaCCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAATGAATAAAGAACAATCTGCTGATGATCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATGTc                                        >  1:612900/1‑149 (MQ=255)
                                                                                                        ttCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGCAAGCGCATGGAGCAACAAAATGAATAAAGAACAATCTGCTGATGATCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGa                                   <  1:734137/149‑1 (MQ=255)
                                                                                                             aCAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAATGAATAAAGAACAATCTGCTGATGATCCCTCCGTGGATCTGATTCGAGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTAc                              >  2:678954/1‑149 (MQ=255)
                                                                                                               aGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAATGAATAAAGAACAATCTGCTGATGATCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCCCCATTTCTATGAGTTAccc                            >  2:30536/1‑149 (MQ=255)
                                                                                                                gAAAGATTTGAAGAAATATTTTAACCGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAATGAATAAAGAACAATCTGCTGATGATCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCt                           <  2:193378/149‑1 (MQ=255)
                                                                                                                                  ttttAACCGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAATGAATAAAGAACAATCTGCTGATGATCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGt         >  1:347862/1‑149 (MQ=255)
                                                                                                                                         cGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAATGAATAAAGAACAATCTGCTGATGATCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAAc  <  1:169664/149‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                    
TGTTTTCAATTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAACCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAATGAATAAAGAACAATCTGCTGATGATCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAGTTACCCTGATGTTGTAATTGCATGTATAGAAC  >  NC_000913/579709‑579994

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