Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,943,9110TC38.1% 13.8 / 15.7 23V192A (GTG→GCG) csdAcysteine sulfinate desulfinase
Reads supporting (aligned to +/- strand):  ref base T (2/11);  new base C (1/7);  total (5/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.93e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCTGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTT‑GCCCGGATCTGGCGCGAGCGATTACCTTTGCTCATTCAGCCGGGATGGTGGTGATGGTTGATGGTGCTC‑AGGGGGCAGTGCATTT‑CCCCGCGGATGTTCAGCAACTGGATATTGATTTCTATGCTTTTTCAGGTCACAAACTGTAT‑GGCCCGACAGGTATCGGCGTGCTGTATGGTAAATCAGAACTGCTGGAGGCGATGTCGCCCTGGCTGGGCGGCGGCAAAATGGTTCACGAAGTGAGTTTTGACGGCTTCACGACTCAATCTGCGCCGTGGAAACTGGAAGCTG  >  NC_000913/2943665‑2944150
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
aTGAGATTATTGTCAGCGTGGCAGAACACCACGCCCACCTCGTCCCCTGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTT‑GCCCGGATCTGGCGCGAGCGATTCCCTTTGCTCATTCAGCCGGGATGgtggtg                                                                                                                                                                                                                                                <  2:110393/251‑1 (MQ=255)
       tattGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCTGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTT‑GCCCGGATCTGGCGCGAGCGATTACCTTTGCTCATTCAGCCGGGATGGTGGtgatggtt                                                                                                                                                                                                                                          <  2:193972/250‑1 (MQ=255)
       tattGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCTGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTT‑GCCCGGATCTGGCGCGAGCGATTACCTTTGCTCATTCAGCCGGGATGGCGGtgatggtt                                                                                                                                                                                                                                          <  2:277803/250‑1 (MQ=255)
           gTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCTGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTT‑GCCCGGATCTGGCGCGAGCGATTACCTTTGCTCATTCAGCCGGGATGGCGGtgatggttgatgg                                                                                                                                                                                                                                     <  2:124742/251‑1 (MQ=255)
                       gAACACCACGCCAACCTCGTCCCCTGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTT‑GCCCGGATCTGGCGCGAGCGATTACCTTTGCTCATTCAGCCGGGATGGCGGTGATGGTTGATGGTGCTC‑AGGGGGc                                                                                                                                                                                                                        <  2:269710/251‑1 (MQ=255)
                                           ccccTGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTT‑GCCCGGATCTGGCGCGAGCGATTACCTTTGCTCATTCAGCCGGGATGGCGGTGATGGTTGATGGTGCTC‑AGGGGGCAGTGCATTT‑CCCCGCGGATg                                                                                                                                                                                                   <  2:164508/251‑1 (MQ=255)
                                                 gCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTT‑GCCCGGATCTGGCGCGAGCGATTACCTTTGCTCATTCAGCCGGGATGGTGGTGATGGTTGATGGTGCTC‑AGGGGGCAGTGCATTT‑CCCCGCGGATGTTCAGc                                                                                                                                                                                             <  2:108738/251‑1 (MQ=255)
                                                    gATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTTGATTTGTTGCCAGAACTGATTACTCCCTGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTT‑GCCCGGATCTGGCGCGAGCGATTACCTTTGCTCATTCAGCCGGGATGGTGGTGATGGTTGATGGTGCTC‑AGGGGGCAGTGCATTT‑CCCCGCGGATGTTCAGCa                                                                                                                                                                                            <  1:214658/249‑1 (MQ=255)
                                                        gTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGGAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTT‑GCCCGGATCTGGCGCGAGCGATTACCTTTGCTCATTCAGCCGGGATGGCGGTGATGGTTGATGGTGCTC‑AGGGGGCAGTGCATTT‑CCCCGCTGATGTTCAGCAACTGGa                                                                                                                                                                                      <  2:107556/251‑1 (MQ=255)
                                                                           ccAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTGCTCCCCGTAGTCGGATGCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTT‑GCACGGATCTGGCGCGAGCGATTACCTTTGCTCATTCAGCCGGGATGGTGGTGATGGTTGATGGTGCTC‑AGGGGGCAGTGCATTT‑CCCCGCGGATGTTCAGCAACTGGATATTGATTTCTATGCtttt                                                                                                                                                                   <  2:210149/251‑1 (MQ=255)
                                                                             aaaGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATGCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTT‑GCGCGGATCTGGGGCGAGCGAGTAGCTTTGCTCATTCAGCCGGGGTGGTGGTGATGGTGGATGGTGCTC‑AGGGGGCAGTGCATTT‑CCCCGCGTATGTTCAGCAACTGGATATTGATTTCTATGCTTTTTc                                                                                                                                                                 <  1:214998/251‑1 (MQ=255)
                                                                              aaGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTTTGCCCG‑ATTTGGGGCGAGGGATTACCTTTGTTCATTCTGCCGGGGGGGTGGGGGGGGGTGATGGTGGTC‑CGGGGGGCGTGCATTT‑TCCCGCGGGTGGTTAGCAACCGGATATTGATTTTTTTGCttttttt                                                                                                                                                                >  1:102617/1‑249 (MQ=255)
                                                                                                                     aTGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTT‑GCCCCGGTCTGGCGCGGGCGGTTACCTTTTCTCATTTAGCCGGGGGGGGGGGGGTGGGTGGTGGGGCTCCAGGGGGC‑GGGCATTTTCCCCCCGGTGGTTC‑GCAACTGGGTATTGATTTTTTTTCTTTTTTTTGTCACAAACTGTATTGGCCCgccgggatcggcggggtg                                                                                                                         >  1:104598/1‑234 (MQ=255)
                                                                                                                                 tgCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTT‑GCCCGGATCTGGCGCGAGCGATTACCTTTGCTCATTCAGCCGGGATGGTGGtgatgg                                                                                                                                                                                                                                            <  2:112632/126‑1 (MQ=255)
                                                                                                                                                  ccccGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTT‑GCCCGGATCTGGCGCGAGCGATTACCTTTGCTCATTCAGCCGGGATGGTGGTGATGGTTGATGGTGCTC‑AGGGGGCAGTGCATTT‑CCCCGCGGATGTTCAGCAACTGGATATTGATTTCTATGCTTTTTCAGGTCACAAACTGTAT‑GGCCCGACAGGTATCGGCGTGCTGTATGGTAAATCAGAACTGCTGGAGGCg                                                                                             <  1:365187/249‑1 (MQ=255)
                                                                                                                                                                                  gTCGATCGTTACTGGCGGTT‑GCCCGGATCTGGCGCGAGCGATTACCTTTGCTCATTCAGCCGGGATGGTGGTGATGGTTGATGGTGCTC‑AGGGGGCAGTGCATTT‑CCCCGCGGATGTTCAGCAACTGGATATTGATTTCTATGCTTTTTCAGGTCACAAACTGTAT‑GGCCCGACAGGTATCGACGTGCTGTATGGTAAATCAGAACTGCTGGAGGCGATGTCGCCCTGGCTGGGCGGCGGCAAAATGGtt                                                            <  1:274672/250‑1 (MQ=255)
                                                                                                                                                                                           tACTGGCGGTT‑GCCCGGATCTGGCGCGAGCGATTACCTTTGCTCATTCAGCCGGGATGGTGGTGATGGTTGATGGTGCTC‑AGGGGGCAGTGCATTT‑CCCCGCGGATGTTCAGCAACTGGATATTGATTTCTATGCTTTTTCAGGTCACAAACTGTAT‑GGCCCGACAGGTATCGGCGTGCTGTATGGTAAATCAGAACTGCTGGAGGCGATGTCGCCCTGGCTGGGCGGCGGCAAAATGGTTCACGAAGTg                                                   <  1:85401/250‑1 (MQ=255)
                                                                                                                                                                                                                 ggCGCGAGCGATTACCTTTGCTCATTCAGCCGGGATGGCGGTGATGGTTGATGGTGCTC‑AGGGGGCAGTGCATTT‑CCCCGCGGATGTTCAGCAACTGGATATTGATTTCTATGCTTTTTCAGGTCACAAACTGTAT‑GGCCCGACAGGTATCGGCGTGCTGTATGGTAAATCAGAACTGCTGGAGGCGATGTCGCCCTGGCTGGGCGGCGGCAAAATGGTTCACGAAGTGAGTTTTGACGGCTTCAc                                  <  1:69512/246‑1 (MQ=255)
                                                                                                                                                                                                                      gatcgtttccctgtGCTCATTCAGCCGGGATGACGGTGATCGTTCAGGGTGCGA‑AATGGGCAGTGCATTC‑CCCCGCGGATGTTCAGCAACTGGATATTGACGTCTATGCGTTTTCAGGTCACAACCTGTAT‑GGCCCGACCGGTATCGGCGTGCTGTATGGTAAATCAGAACTGCTGGAGGCGATGTCGCCCTGGCTGGGCGGCGGCAAAAAGGTTCACGAAGTGAGGTTTGACGGCTTCCCTACTCAATCt                        <  2:143542/238‑1 (MQ=255)
                                                                                                                                                                                                                         cgaTTACCTTTGCTCATTCAGCCGGGATGGCGCTGATGGTTGATGGTGCTC‑AGGGGGCCGTGCATTT‑CCCCGCGGATGTTCAGCAACTGGATATTGATTTCTATGCTTTTTCAGGTCACAAACTGTAT‑GGCCCGACAGGTATCGGCGTGCTGTATGGTAAATCAGCACTGCTGGAGGCGATGTCGCCCTGGCTGGGCGGCGGCAAAATGGTTCACGAAGTGAGTTTTGACGGCTTCACGACTCAATCTgcg                     >  2:358096/1‑251 (MQ=255)
                                                                                                                                                                                                                          gaTTACCTTTGCTCATTCAGCCGGGATGGTGGTGATGGTTGATGGTGCTC‑AGGGGGCAGTGCATTT‑CCCCGCGGATGTTCAGCAACTGGATATTGATTTCTATGCTTTTTCAGGTCACAAACTGTAT‑GGCCCGACAGGTATCGGCGTGCTGTATGGTAAATCAGAACTGCTGGAGGCGATGTCGCCCTGGCTGGGCGGCGGCAAAATGGTTCACGAAGTGAGTTTTGACGGCTTCACGACTCAATCTgcgc                    >  1:115665/1‑251 (MQ=255)
                                                                                                                                                                                                                                cTTTGCTCATTCCGCGGGGCTGGAGGTGATGGATGATGGTGATC‑AGGGGGCCGTGCATTT‑CCCCGCGGCTGTTCAGCAACTGGATATTGATTTCTATGCTTTTTCAGGCCACAAACTGTAT‑GGCCCGACCGGTCTCGGCGTGCTGAATGGTAAATCAGAACTGATGGAGGCGATGTCGCCCTGGCTGGGCGGCGGCAAAATGGTTCACGAAGTGAGTTTTGACGGCTTAACGACCCACTCTGCGCCGTgg               >  2:215784/1‑250 (MQ=255)
                                                                                                                                                                                                                                    gCTCATTCAGCCGGGATGGTGGTGATGGTTGATGGTGCTC‑AGGGGGCAGTGCATTT‑CCCCGCGGATGTTCAGCAACTGGATATTGATTTCTATGCTTTTTCAGGTCACAAACTGTAT‑GGCCCGACAGGTATCGGCGTGCTGTATGGTAAATCAGAACTGCTGGAGGCGATGTCGCCCTGGCTGGGCGGCGGCAAAATGGTTCACGAAGTGAGTGTTGACGGCTTCACGACTCAATCTGCGCCGTGGAAAc           <  2:172163/250‑1 (MQ=255)
                                                                                                                                                                                                                                            aGCCGGGATGGTGGTGATGGTTGATGGTGCTC‑AGGGGGCAGTGCATTT‑CCCCGCGGATGTTCAGCAACTGGATATTGATTTCTATGCTTTTTCAGGTCACAAACTGTAT‑GGCCCGACAGGTATCGGCGTGCTGTATGGTAAATCAGAACTGCTGGAGGCGATGTCGCCCTGGCTGGGCGGCGGCAAAATGGTTCACGAAGTGAGTTTTGACGGCTTCACGACTCAATCTGCGCCGTGGAAACTGGAAGCTg  <  2:283329/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
ATGAGATTATTGTCAGCGTGGCAGAACACCACGCCAACCTCGTCCCCTGGCTGATGGTCGCCCAACAAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTT‑GCCCGGATCTGGCGCGAGCGATTACCTTTGCTCATTCAGCCGGGATGGTGGTGATGGTTGATGGTGCTC‑AGGGGGCAGTGCATTT‑CCCCGCGGATGTTCAGCAACTGGATATTGATTTCTATGCTTTTTCAGGTCACAAACTGTAT‑GGCCCGACAGGTATCGGCGTGCTGTATGGTAAATCAGAACTGCTGGAGGCGATGTCGCCCTGGCTGGGCGGCGGCAAAATGGTTCACGAAGTGAGTTTTGACGGCTTCACGACTCAATCTGCGCCGTGGAAACTGGAAGCTG  >  NC_000913/2943665‑2944150

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