Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077931,154,5830TC73.5% 23.3 / 11.3 34intergenic (+386/‑168)SAUSA300_RS05695/SAUSA300_RS05700fibrinogen‑binding protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (2/7);  new base C (1/24);  total (3/31)
Fisher's exact test for biased strand distribution p-value = 1.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.90e-01

CTTATAACAAGTTTTTTCAATACTCGTGGGGCCACAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGGCCCCAATATAGAAGCTGGCGAAAAGTCAGCTTACAAAAATGTGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGAACCCCAATTTCTACAGACAATGCA  >  NC_007793/1154489‑1154677
                                                                                              |                                                                                              
cTTATAACAAGTTTTTTCAATACTCGTGGGGCCACAACAGAGAGAAATAGGATAACCAATTCCAACAGACAATGCAAGTTGGC‑GGGGCCCCAACAgagag                                                                                          <  1:4625926/100‑2 (MQ=21)
cTTATAACAAGTTTTTTCAATACTCGTGGGGCCACAACAAAGAGAAATAGGATCACCAATTCCAACAGACAATGCAAGTTGGC‑GGGGCCCCAACAgagag                                                                                          <  1:4165763/100‑2 (MQ=21)
  tataACAAGTTTTTTCAATACTCGTGGGGCCACAACAGAGAGAAATAGGATCACCAATTCCAACAGACAATGCAAGTTGGC‑GGGGCCCCAACAgagagaa                                                                                        <  1:2025911/100‑4 (MQ=21)
    taACAAGTTTTTTCAATACTCGTGGGGCCACAACAGAGAGAAATAGGATCACCAATTCCAACAGACAATGCAAGTTGGC‑GGGGCCCCAACAgagagaaat                                                                                      <  1:1155660/100‑6 (MQ=21)
    taACAAGTTTTTTCAATACTCGTGGGGCCACAACAGAGAGAAATAGGATCACCAATTCCAACAGACAATGCAAGTTGGC‑GGGGCCCCAACAgagagaaat                                                                                      <  1:4834307/100‑6 (MQ=21)
     aaCAAGTTTTTTCAATACTCGTGGGGCCACAACAGAGAGAAATAGGTTCACCAATTCCAACAGACAATGCAAGTTGGC‑GGGGCCCCAACAgagagaaata                                                                                     <  1:1285760/100‑7 (MQ=21)
       cAAGTTTTTTCAATACTCGTGGGGCCACAACAGAGGGAAATAGGATCACCAATTCCAACAGACAATGCAAGTTGGC‑GGGGCCCCAACAgagagaaatagg                                                                                   <  1:1256871/100‑9 (MQ=21)
       cAAGTTTTTTCAATACTCGTGGGGCCACAACAGAGAGAAATAGGATCACCAATTCCAACAGACAATGCAAGTTGGC‑GGGGCCCCAACAgagagaaatagg                                                                                   <  1:3771437/100‑9 (MQ=21)
       cAAGTTTTTTCAATACTCGTGGGGACACAACAGAGAGAAATAGGATCACCAATTCCAACAGACAATGCAAGTTGGC‑GGGGCCCCAACAgagagaaatagg                                                                                   <  1:3060780/100‑9 (MQ=21)
         aGTTTTTTCAATACTCGTGGGGCCACAACAGAGAGAAATAGGATCACCAATTCCAACAGACAATGCAAGTTGGC‑GGGGCCCCAACAgagagaaataggat                                                                                 <  1:1303225/100‑11 (MQ=18)
         aGTTTTTTCAATACTCGTGGGGCCACAACAGAGAGAAATAGGATCACCAATTCCAACAGACAATGCAAGTTGGC‑GGGGCCCCAACAgagagaaataggag                                                                                 <  1:4618377/100‑11 (MQ=18)
                 gcaagttggcggggCCCCAACAGAGAGAAATAGGATCACCAATTTCAACAGGCAATGCAAGTTGGGGTGGGGCCCAACACAGAAGCTGGCGAAAAGTCAg                                                                          >  1:555075/11‑100 (MQ=11)
                 gcaagttggcggggCCCCAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGGCCCCAATATAGAAGCTGGCGAAAAGTCAg                                                                          >  1:3381702/11‑100 (MQ=21)
                               ccccaacaGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGGCCCCAATATAGAAGCTGGCGAAAAGTCAGCTTACAAAAAtgtg                                                            <  1:3931335/97‑1 (MQ=255)
                                cccaacaGAGAGAAATAGGATCACCAATTTCAACAGACAATGCAAGTTGGGGTGGGGCCCAACACAGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGc                                                           <  1:675480/98‑1 (MQ=14)
                                cccaacaGAGAGAAATAGGATCACCAATTTCAACAGACAATGCAAGTTGGGGTGGGGCCCAACACAGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGc                                                           <  1:2848739/98‑1 (MQ=14)
                                 ccaacaGAGAGAAATAGGATCACCAATTTCAACAGACAATGCAAGTTGGGGTGGGGCCCAACACAGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCa                                                          <  1:5627034/99‑1 (MQ=14)
                                  caacaGAGAGAAATAGGATCACCAATTTCAACAGACAATGCAAGTTGGGGTGGGGCCCAACACAGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCaa                                                         <  1:4997485/100‑1 (MQ=14)
                                     caGAGAGAAATAGGATCACCAATTTCAACAGACAATGCAAGTTGGGGTGGGGCCCAACACAGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGtt                                                      <  1:3223045/100‑1 (MQ=14)
                                            aaaTAGGATCACCAATTTCAACAGACAATGCAAGTTGGGGTGGGGCCCAACACAGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCgggg                                               <  1:957319/100‑1 (MQ=14)
                                                       ccAATTCCAACAAACAATGCAAGTTGGCGGGGGCCCCAATATAGAAGCTGGCGAAAAGTCAGCTTACAAAAATGTGCAAGTTGGCGGGGCCCCAACATag                                    <  1:3781607/100‑1 (MQ=255)
                                                       ccAATTCCAACAAACAATGCAAGTTGGCGGGGGCCCCAATATAGAAGCTGGCGAAAAGTCAGCTTACAAAAATGTGCAAGTTGGCGGGGCCCCAACATag                                    <  1:2001943/100‑1 (MQ=255)
                                                       ccAAGTCCAACAAACAATGCAAGTTGGCGGGGGCCCCAATATAGAAGCTGGCGAAAAGTCAGCTTACAAAAATGTGCAAGTTGGCGGGGCCCCAACATag                                    <  1:5152859/100‑1 (MQ=255)
                                                                acagacaaTGCAAGTTGGGGTGGGGCCCAACACAGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAgg                           <  1:282221/100‑1 (MQ=11)
                                                                acagacaaTGCAAGTTGGGGTGGGGCCCAACACAGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAgg                           <  1:2430767/100‑1 (MQ=11)
                                                                acaaacaaTGCAAGTTGGCGGGGGCCCCAATATAGAAGCTGGCGAAAAGTCAGCTTACAAAAATGTGCAAGTTGGCGGGGCCCCAACATAGAGAAATAgg                           <  1:2142470/100‑1 (MQ=255)
                                                                 cagacaaTGCAAGTTGGGGTGGGGCCCAACACAGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAggt                          <  1:3787344/97‑2 (MQ=11)
                                                                 cagacaaTGCAAGTTGGGGTGGGGCCCAACACAGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAggg                          <  1:2690014/97‑2 (MQ=11)
                                                                       aTGCAAGTTGGCGGGGGCCCCAATATAGAAGCTGGCGAAAAGTCAGCTTACAAAAATGTGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGAAcccct                    <  1:2498991/100‑2 (MQ=255)
                                                                        tGCAAGTTGGGGTGGGGCCCAGCACAGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCaa                   <  1:932478/100‑1 (MQ=14)
                                                                        tGCAAGTTGGGGTGGGGCCCAACACAGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCaa                   <  1:2306233/100‑1 (MQ=14)
                                                                         gCAAGTTGGCGGGGGCCCCAATATAGAAGCTGGCGAAAAGTCAGCTTACAAAAATGTGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGAACCCCAAt                  >  1:1937644/1‑100 (MQ=255)
                                                                                         cccAATATAGAAGCTGGCGAAAAGTCAGCTTACAAAAATGTGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGAACCCCAATTTCTACAGACAATGCa  <  1:3401722/100‑1 (MQ=255)
                                                                                         cccAACACAGAAGCTGGCGAAAAGTCAGCATACAAAAATGTGCAAGTTGGCGGGGCCCCAACAGAGAGAAATAGGATCACCAATTTCAACAGACAATGCa  <  1:824813/100‑1 (MQ=14)
                                                                                              |                                                                                              
CTTATAACAAGTTTTTTCAATACTCGTGGGGCCACAACAGAGAGAAATAGGATCACCAATTCCAACAAACAATGCAAGTTGGCGGGGGCCCCAATATAGAAGCTGGCGAAAAGTCAGCTTACAAAAATGTGCAAGTTGGCGGGGCCCCAACATAGAGAAATTGGAACCCCAATTTCTACAGACAATGCA  >  NC_007793/1154489‑1154677

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