Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913580,3020CT55.6% ‑3.8 / 42.2 37R48R (AGG→AGA) ylcIDUF3950 domain‑containing protein YlcI
Reads supporting (aligned to +/- strand):  ref base C (4/12);  new base T (7/13);  total (12/25)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.15e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCGTATTCAGTGTCGCTGATTTGTATTGTCTGAAGTTGTTTTTACGTTAAGTTGATGCAGATCAATTAATACG  >  NC_000913/580085‑580527
                                                                                                                                                                                                                         |                                                                                                                                                                                                                                 
gATCACCATAACTGCTAATCAGTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATGTCTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGGCCTGTTCGGAGGGAAGAACGCGGGATGTTCAATCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACgg                                                                                                                                                                                                                          <  1:129453/227‑1 (MQ=255)
     ncATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGTCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAAtt                                                                                                                                                                                             >  1:197613/2‑251 (MQ=255)
             gCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTGCAATATCGTCCGGTGCGAAGGAAAGAACGGGGGATTTTCATTCTTCATCACTTTTAATTTAGGAATACGCTCTCTTTTCGGACGTGTGCCCCCGAGGCCAGGCTTCAATGACCCAGGCTGAGAAATTTCCGGAcc                                                                                                                                                                                     <  1:316209/251‑1 (MQ=255)
             gCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGGGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCGGACGTTAGCCTCCTACGGCAGGCTTCAATGCCCCAGGCTGAGAAATTCCCGGAcc                                                                                                                                                                                     <  1:133700/251‑1 (MQ=255)
                         aaCTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACCTAATTACTGCAATGCCCTCGTAATTAAGTTAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGGGTGATGTTCATTCTTCATCACGTTTAATTGATGTATATGCTCTCTTTGCTGACGTGAGCCTCGGACTGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGAc                                                                                                                                                                                      <  1:166444/238‑1 (MQ=255)
                             actTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTTTTCGGAGGGAAGATCGCGGGATGTTCATTCTTCATCACTTTTAATGTATGTATATGCTCTCTTTTCTTACTTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAaga                                                                                                                                                                       <  1:302527/249‑1 (MQ=255)
                                tagtcTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGTCTCCGAGGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCgg                                                                                                                                                                                        <  1:98802/225‑1 (MQ=255)
                                  gtcTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGTCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGcn                                                                                                                                                                    <  1:83536/245‑2 (MQ=255)
                                         aCAGAGCCAACACTCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGAAATTAAGTGAAGTTACAATATCCTCCGGGTCGGACGGAACAACTCGGGATGCCCATCCTTCATCCCTTTTAATTGATGTATATGCTCCCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGTCTGAGAAATTCCCGGa                                                                                                                                                                                       <  1:330108/221‑1 (MQ=255)
                                                             gtcactgtcaGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGTCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAgn                                                                                                                                       <  1:96381/249‑2 (MQ=255)
                                                                 ctgtcaGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGACTTTACAATATCGTCCTGTTCGGAGGGCCGAACGCGGGATGTTCATTCTTCCTCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTCGTCTCCGACGGCAGGCTTCAATGACCCCGGCTGAGAAATTCCCGGACCCTTTTTGCTCAag                                                                                                                                                                        >  1:281856/1‑212 (MQ=38)
                                                                            gTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGTCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCGCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCggg                                                                                                                      <  1:303765/251‑1 (MQ=255)
                                                                              ggTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATACTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGgtt                                                                                                                    >  1:157316/1‑251 (MQ=255)
                                                                                           cTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGATGTTCTGCggg                                                                                                       >  1:15420/1‑251 (MQ=255)
                                                                                             cAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTGAATTGATGTATATGCTCTCTTGTCTGACGTTAGTCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTcc                                                                                                                                                                                           <  1:150162/165‑1 (MQ=255)
                                                                                                   cTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTcnc                                                                                                                                                                                          <  1:255588/160‑1 (MQ=255)
                                                                                                          gCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCGCTTGTAATGGATGTATATGCTCTCTTTTCTGACGTTAGTCTCCGACGGCAGGCTTCAATGACCCAGGCTGATAAATTCCCGGAACCTCTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCTGGTTCTGTTCTTCGTTg                                                                                        <  1:246933/251‑1 (MQ=21)
                                                                                                                     nnaGNNNNNTTANAATATCGTTCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTGTAATTGATGTATATGCTCTCTTTTCTGACGTTAGTCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGAc                                                                                                                                                                                      <  2:197613/144‑1 (MQ=38)
                                                                                                                      aaGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTg                                                                             <  1:151386/250‑1 (MQ=255)
                                                                                                                       aGTGAATTTACAATAGCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCAGTCTTCATCACTTTTAATTGATGTAGATGCTCTCTTTTCTGACGTTAGTCTCCGACGGCAGGCTTCAATGACCCAGGCTga                                                                                                                                                                                                   <  1:109194/131‑1 (MQ=255)
                                                                                                                        gTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGTCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTCGTTGACATGAGGTTGccc                                                                          <  1:26223/251‑1 (MQ=255)
                                                                                                                                        gTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGTCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTCGTTGACATGa                                                                                  >  1:89978/1‑227 (MQ=255)
                                                                                                                                        gTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGAAGTATATGCTCTCTTTTCTGACGTTAGTCTCCGACGGCAGGCTTCAATGACCCAGGCTAAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTCGTTGACATGAGGTTTCCACGTATTCAGTGTCGCt                                                          >  1:48059/1‑251 (MQ=255)
                                                                                                                                          ncTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGAccc                                                                                                                                                                                    >  1:205214/2‑127 (MQ=255)
                                                                                                                                          ccTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGAccc                                                                                                                                                                                    <  2:205214/127‑1 (MQ=255)
                                                                                                                                           cTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGACGCATAAGCTCTCTTTTCTGACGTTAGTCTCCGACGGCCGGCTTCAATGACCCAGGATGAGAAATTCCCGGACCCTTTTTGCTCCAGAGCGATGTTCATTTGTTCAATAATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCACGTATTCAGTGTCGCTGAt                                                       >  1:281029/1‑251 (MQ=18)
                                                                                                                                             gTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGTCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCttcg                                                                                           <  1:166844/213‑3 (MQ=255)
                                                                                                                                                        agaaCGCGGGATGTTCATTCTTCATCTCTTTTAATTGATGTATATGCTCTCTTTTCTGACGGTAGTCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTCGTTGACATGAGGTTGCGCCGTATGCAGTTTCGCTTATTTTTATTGTCTGa                                          <  1:87332/251‑1 (MQ=17)
                                                                                                                                                         gaaCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACNNAGGCTgnna                                                                                                                                                                                                 >  2:1912/1‑99 (MQ=255)
                                                                                                                                                         gaaCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGgttctgttc                                                                                               <  1:1912/197‑1 (MQ=255)
                                                                                                                                                              cgNGATGTNCATTNTTCATCACTTTTAATNGATGTATAGGCTCTCTTTTNNGNCGTTANCCTcc                                                                                                                                                                                                                               <  2:135951/64‑1 (MQ=255)
                                                                                                                                                              cgGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGACTcc                                                                                                                                                                                                                               >  1:135951/1‑64 (MQ=255)
                                                                                                                                                                     ttCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTGGCCCCGTATTCAGTGTCGc                                                           <  1:270548/221‑1 (MQ=255)
                                                                                                                                                                                                 atatGCTCTCTTTTCTGACGTTAGTCTCCGACGGCAGGCTTCAATGACCCAGTCTGAGAAATTCACGGTCCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGTGTTCTGTTCTTCGTTGACATGAGTTTGCCCCGTATTCAGTGTCGCTGATTTGTATTGTCTGAATTTGTTTTTACGTTAAGTTGATGCAGATCAATTAATAcn  <  1:63240/250‑2 (MQ=11)
                                                                                                                                                                                                    tGCTCTCTTTTCTGACGTTAGTCTCCGACGGCAGGCTTCAATGAcc                                                                                                                                                                                                           <  1:18822/46‑1 (MQ=255)
                                                                                                                                                                                                    tGCTCTCTTTTCTGACGTTAGTCTCCGACGGCAGGCTTCAATGAcc                                                                                                                                                                                                           >  2:18822/1‑46 (MQ=255)
                                                                                                                                                                                                      ctctctTTTCTGACGTTAGTCTCCGACGGCAGGCTTCACTGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCCATCCTTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTCGTTGACATGAGGTTGCCCCGTATTCAGTGTCGCTGATTTGTATTGTCTGAAGTTGttttt                                >  1:300275/1‑215 (MQ=255)
                                                                                                                                                                                                                         |                                                                                                                                                                                                                                 
GATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCAGGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCGTATTCAGTGTCGCTGATTTGTATTGTCTGAAGTTGTTTTTACGTTAAGTTGATGCAGATCAATTAATACG  >  NC_000913/580085‑580527

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