Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077791,417,5240TC51.9% 34.3 / 79.3 52D526G (GAT→GGT) recEexonuclease VIII, 5' ‑> 3' specific dsDNA exonuclease
Reads supporting (aligned to +/- strand):  ref base T (20/5);  new base C (0/27);  total (20/32)
Fisher's exact test for biased strand distribution p-value = 4.22e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.11e-08
Rejected as consensus: Frequency below/above cutoff threshold.

CCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTCTG  >  NC_007779/1417429‑1417623
                                                                                               |                                                                                                   
ccGGATCATTTTCCTTAGCTTCAACCTGAATCTCTTTTGCGAATGGTTGCTAGTATGTTGCGTCACGCATTACCGCACCACAGGCAGGGCAGTTATCCCCg                                                                                                <  1:251310/101‑1 (MQ=255)
 cGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGc                                                                                               <  1:1400691/101‑1 (MQ=255)
   gATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCa                                                                                             >  1:618605/1‑101 (MQ=255)
       aTTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCCCAGTCAGGGCAGTTATCCCCGCCAGTCt                                                                                         >  3:290671/1‑101 (MQ=255)
         tttCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTgg                                                                                       >  3:1060077/1‑101 (MQ=255)
         tttCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTgg                                                                                       >  3:400593/1‑101 (MQ=255)
         tttCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTgg                                                                                       >  1:620303/1‑101 (MQ=255)
          ttCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGc                                                                                      >  3:1184831/1‑101 (MQ=255)
          ttCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGc                                                                                      >  3:24436/1‑101 (MQ=255)
           tCCTTAGCTTCACCCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCCCCGCCCCACAGTCGGGGCGGTTACCCCCGCCAGTCTGGcc                                                                                     <  4:29568/101‑1 (MQ=255)
           tCCTTAGCTTCACCCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCCCCCATCCCCGCCCCACAGTAAGGGAATTTACCCCCGCCAGTCTGGcc                                                                                     <  4:87319/101‑1 (MQ=255)
            ccTTAGCTTCACCCTGCCTTTTTTCATGGAATGTTTCCTGGTTTTTTGCGTCGCCCATCACCGCCCCACAGTCAGGCCGTTTACCCCCGCCAGTCTGGCCg                                                                                    <  1:138460/101‑1 (MQ=255)
             cTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGc                                                                                   >  3:1499101/1‑101 (MQ=255)
                aGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAgg                                                                                >  3:1591819/1‑101 (MQ=255)
                 gCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGc                                                                               >  1:778484/1‑101 (MQ=255)
                  cTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCa                                                                              >  1:913655/1‑101 (MQ=255)
                    tCACCCTGACTCTCTTCATCGAATTTTTCCTGGTATGTTGCGTCCCCCTTCCCCGCCCCAAAGTCAGGGCAGTTACCCCCGCCAGTCTGGCCGCAGGCAtt                                                                            <  1:1467963/101‑1 (MQ=255)
                    tCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCAtt                                                                            >  3:1603084/1‑101 (MQ=255)
                      acccTGCCTCTCTTCATCGAATGTTTCCTGGTATGTTGCGCCCCCCATCCCCGCCCCACAGTCGGGGCAGTTACCCCCGCCAGTCTGGCCGCAGGCATTGc                                                                          <  1:75146/99‑1 (MQ=255)
                       cccTGCCTCTCTTCATCGATTTTTTCCTGTTATGTTGCGTCCCCCATCCCCCCCCCACAGTAGGGGCATTTACCCCCGCCAGTCTGGCCGCAGGCATTGCa                                                                         <  1:1198696/100‑1 (MQ=255)
                        ccTGCCTCTCTTCATCGATTGTTTCCTGGTATGTTGCGCCCCCCACCCCCGCCCCACAGTCAGGGCAGTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAg                                                                        <  3:822090/101‑1 (MQ=255)
                                    catctattttttCTTGGTTTGTTGCGCCCCCCACCCCCGCCCCACAGTCGGGGCGTTTCCCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGt                                                            <  1:1061127/92‑1 (MQ=255)
                                    cATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGt                                                            >  3:1135303/1‑101 (MQ=255)
                                      tCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTc                                                          >  1:31447/1‑101 (MQ=255)
                                        gattGTTTCCTGGTATGTTGCGTCCCCCATCCCCGCACCACAGTCGGGGCATTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCt                                                        <  3:506253/98‑1 (MQ=255)
                                           tttttCTTGGTTTGTTGCGTCCCCCACCCCCGCCCCACAGTCAGGGCAGTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGtt                                                     <  1:633040/99‑1 (MQ=255)
                                           tGTTGCCTGGTTTGTTGCGCCCCCCACCCCCGCCCCAGAGTAAGGGCAGTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGtt                                                     <  3:276399/101‑1 (MQ=255)
                                            gTTTCTTGGTATGTTGCGTCGCCCATCCCCGCACCACAGTCAGGGCAGTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTg                                                    <  2:43439/101‑1 (MQ=255)
                                            gTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTg                                                    <  2:44366/101‑1 (MQ=255)
                                               tCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCac                                                 >  1:377400/1‑101 (MQ=255)
                                                ccTGGTTGGTTGCGCCGCCCACCCCCGCCCCACGGTCGGGGCAGTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCact                                                <  4:43263/101‑1 (MQ=255)
                                                ccTGGTATGTTGCGCCGCCCACCCCCGCCCCACAGTCGGGGCAGTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGGGCCTGGGGCact                                                <  4:16064/101‑1 (MQ=255)
                                                   ggTTTGTTGCGCCCCCCATCCCCGCCCCACAGTCGGGGCCGTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCactact                                             <  1:1050188/101‑1 (MQ=255)
                                                       tgtgGCGCCCCCCACCCCCGCCCCACAGTCGGGGCAGTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCt                                         <  1:55136/101‑1 (MQ=255)
                                                       tGTTGCGCCGCCCACCCCCGCCCCACAGTCGGGGCGGTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCt                                         <  4:39535/101‑1 (MQ=255)
                                                        gttggctCGCCCACCACCCCCCCACAGCCGGGGCAGTTTCCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTc                                        <  1:1483866/95‑1 (MQ=255)
                                                         ttGCGTCGCCCATCCCCCCACCACAGTCGGGGCAGTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCa                                       <  3:822248/101‑1 (MQ=255)
                                                             gccccccATCCCCGCACCACAGTCAGGGCAGTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGtt                                   <  3:786540/97‑1 (MQ=255)
                                                               cgccaacccccGCCCCACAGTAAGTTCTTTTACCCCCGCCAGTCTGGCCGCAGGCATNGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGa                                 <  3:1276589/92‑1 (MQ=255)
                                                                gCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGAt                                >  3:1091962/1‑101 (MQ=255)
                                                                 cccATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGAtt                               >  3:130992/1‑101 (MQ=255)
                                                                   cacccccGCACCACAGTCAGGGCAGTTCCCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCa                             <  3:592827/96‑1 (MQ=255)
                                                                             ccccaGTTAGGGCCTTTTTCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCt                   <  1:684613/98‑1 (MQ=255)
                                                                              cacaGTCGGGGCAGTTCCCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTg                  <  1:45723/101‑1 (MQ=255)
                                                                              cacaGTCAGGGCATTTCCCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTg                  <  3:633201/101‑1 (MQ=255)
                                                                              cacaGCCGGGGCAGTTCCCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTg                  <  3:908114/101‑1 (MQ=255)
                                                                               acaGTCAGGCCAGTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTgg                 <  3:1585445/101‑1 (MQ=255)
                                                                                   tgaGGCCAGTTCCCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCNGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGtttt             <  3:1453213/99‑1 (MQ=255)
                                                                                       ggCATTTACCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTcc         <  1:1225450/101‑1 (MQ=255)
                                                                                          aGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGtt      >  1:1038632/1‑101 (MQ=255)
                                                                                             tATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTct   >  1:962607/1‑101 (MQ=255)
                                                                                              aTCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTctg  >  1:1064740/1‑101 (MQ=255)
                                                                                               |                                                                                                   
CCGGATCATTTTCCTTAGCTTCAACCTGACTCTCTTCATCGAATGTTTCCTGGTATGTTGCGTCGCCCATCACCGCACCACAGTCAGGGCAGTTATCCCCGCCAGTCTGGCCGCAGGCATTGCAGGCTATTTCCGGTTCCTGTTGCACTACTGGCTCAGGTTGATTCATATCTGGGCTGGTTTTTTCCGTTTCTG  >  NC_007779/1417429‑1417623

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