Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1225,525,8830AC36.7% 62.6 / 10.3 49intergenic (‑53/+102)livH_5/braC_2High‑affinity branched‑chain amino acidtransport system permease protein LivH/Leucine‑, isoleucine‑, valine‑, threonine‑, and alanine‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (8/23);  new base C (17/1);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 3.05e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGTAGCCGATGGCGATCAGGCCATAGACCGACCCGAGGGTCAGGCCGTTGACCAGTTGCTGCAGGAAAATACCATCCATAACGC‑AAT‑CTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCACCTGTGCGGGGCCCTCGCAGCGCAGAACTGTGCAGATCTACGAATAAAGACAGGTACCGCGGGGCTTGGGCCCGAGCAACCGCCCGGGCCGCGGGCCGTTGTTATTTTTGTTTTTCCAGCTGGTGGTATTTAC  >  ACB122/5525752‑5526016
                                                                                                                                     |                                                                                                                                     
tgtAGCCGATGGCGATCAGGCCATAGACCGACCCGAGGGTCAGGCCGTTGACCAGTTGCTGCAGGAAAATACCATCCATAACGC‑AAT‑CTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCccccc                                                                                                                                    >  3:116883/1‑134 (MQ=255)
tgtAGCCGATGGCGATCAGGCCATAGACCGACCCGAGGGTCAGGCCGTTGACCAGTTGCTGCAGGAAAATACCATCCATAACGC‑AAT‑CTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCCCCt                                                                                                                                    >  2:68603/1‑135 (MQ=255)
tgtAGCCGATGGCGATCAGGCCATAGACCGACCCGAGGGTCAGGCCGTTGACCAGTTGCTGCAGGAAAATACCATCCATAACGC‑AAT‑CTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCACCt                                                                                                                                    >  6:160341/1‑135 (MQ=255)
tgtAGCCGATGGCGATCAGGCCATAGACCGACCCGAGGGTCAGGCCGTTGACCAGTTGCTGCAGGAAAATACAATCCATAACGC‑AAT‑CTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCACCt                                                                                                                                    >  4:291955/1‑135 (MQ=255)
 gtAGCCGATGGCGATCAGGCCATAGACCGACCCGAGGGTCAGGCCGTTGACCAGTTGCTGCAGGAAAATACCATCCATACCGC‑AAT‑CTCACCTGATTGAGATTGCACGAGCGCCGCCCACGGCGCGGGCCACCtg                                                                                                                                   >  6:308161/1‑135 (MQ=255)
 gtAGCCGATGGCGATCAGGCCATAGACCGACCCGAGGGTCAGGCCGTTGACCAGTTGCTGCAGGAAAATACCATCCATAACGC‑AAT‑CTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCccccag                                                                                                                                   >  2:11885/1‑133 (MQ=255)
     ccGATGGCGATCTGGCCATAGACCGACCCGAGGGTCAGGCCGTTGACCAGTTGCTGCAGGAAAATACCATCCATAACGC‑AAT‑CTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCCCCTGTGCg                                                                                                                               >  5:184371/1‑135 (MQ=255)
        aTGGCGATCAGGCCATAGACCGACCCGAGGGTCAGGCCGTTGCCCAGTTGCTGCAGGAAAATACCATCCATAACGC‑AAT‑CTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCCCCTGTGCgggg                                                                                                                            >  5:135275/1‑135 (MQ=255)
           gCGATCAGGCCATAGACCGACCCGTGGGTCAGGCCGTTGACCAGTTGCTGCAGGAAAATACCATCAAAAACGC‑AAC‑CAACCCTGATTTAGATTGCACTAGCGCCGACCACGGCGCGGGCCACCTGTGCGGGGccc                                                                                                                         <  8:84891/135‑1 (MQ=255)
           gCGATCAGGCCATAGACCGACCCGAGGGTCAGGCCGTTGACCAGTTGCTGCAGGAAAATACCATCCATAACGC‑AAT‑CTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCACCTGTGCGGGAccc                                                                                                                         <  8:91101/135‑1 (MQ=255)
           gCGAGCCGGCCATACACCGCCCCGGGGGGCAGGCCGTTGCCCAGTGGCTCCAGCAAAATCCCATCCATACCGC‑AAT‑CCCACCTGATTGACATTGCACGAGCCCCGCCCCCGGCGCGGGCCACCTGTGCGGGGccc                                                                                                                         <  4:116883/135‑1 (MQ=255)
           agagacAGGCGATAGACCTCCCCGCGGGTCAGGGCGTTGACCAGTGGCTGCAGTAAAATACGATCCATAACTC‑AAT‑CTCACCTTATTGTGATTGCAGGAGCGCCGACCACGGCGCGGGCCACCTGTGCGGGGccc                                                                                                                         <  8:377741/130‑1 (MQ=255)
             gATCAGGCCATAGACCGACCCGAGGGTCAGGCCGTTGACCAGTTGCTGCAGGAAAATACCATCCATAACGC‑AAT‑CTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGGCCCCGGTGCGGGGcccgc                                                                                                                       >  5:51989/1‑133 (MQ=255)
             gATCAGGCCATAGACCGACCCGAGGGTCAGGCCGTTGACCAGTTGCTGCAGGAAAATACCATCCATAACGC‑AAT‑CTCACCTGATTGAGATTGCACGAGAGCCGACCACGGCGCGGGCCCCCGGTGCGGGGgccca                                                                                                                       >  5:176023/1‑133 (MQ=255)
                cAGGCCATAGACCGACCCGAGGGTCAGGCCGTTGACCAGTTGCTGCAGGAAAATACCATCCATAACGC‑AAT‑CTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCACCTGTGCGGGGCCCTcgca                                                                                                                    <  3:305509/135‑1 (MQ=255)
                    ccATAGCCCCACCGGTGGGCCTGGCCGTGGACCCGTTGCTGATCTAAAATCCCACCAATAAC‑C‑AATCCGCCCCTTATTCCCATTGCAGGCCCGCCCACCCCGCCGCGGGCCACCTGTGCGGGGCCCTcgcagcgca                                                                                                               <  8:32487/136‑1 (MQ=255)
                     aaTAGACCGACCCGAGGGCCCGGCCGTTGCCCAGTTGCGGCAGTAAAAACCCACCCATAACCC‑AAT‑CACACCTGATTGAGATTGCACGCGCGCCGACCACGGCGCGGGCCACCTGTGCGGGGCCCTcgcagcgca                                                                                                               <  6:169744/134‑1 (MQ=255)
                                ccGAGGGTCAGGCCGTTGACCAGTTGCTGCAGGAAAATACCATCCATAACGC‑AAT‑CTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCCCCTGTGCGGGGCCCTCGCCCCGCAGCACAGTGcc                                                                                                     >  6:92914/1‑133 (MQ=255)
                                ccGAGGGTCAGGCCGTTGACCAGTTGCTGCAGGAAAATACCATCCATAACGC‑AAT‑CTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCCCCTGTGCGGGGCCCTCGCAGCGCAGAACtgtg                                                                                                       >  8:137419/1‑132 (MQ=255)
                                ccGAGGGTCAGGCCGTTGACCAGTTGCTGCAGGAAAATACCATCCATAACGC‑AAT‑CTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCCCCGGTGGGGGGCCCTCGCAGCGCAGAACtgt                                                                                                        >  2:231270/1‑131 (MQ=255)
                                ccGAGGGTCAGGCCGTTGACCAGTTGCTGCAGGAAAATACCATCCATAACGC‑AAT‑CTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCCCCCGTGCGGGGGCCGCGCCGCGCCGAAcggggcc                                                                                                     >  7:225677/1‑128 (MQ=255)
                                 cGAGGGTCAGGCCGTTGACCAGTTGCTGCAGGAAAATACCATCCATAACGC‑AAT‑CTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCCCCAGTGCGGGGCCCGCCCCGCGCAGCACtggggagg                                                                                                   >  3:86201/1‑129 (MQ=255)
                                  gAGGGTCAGGCCGTTGACCAGTGGCGCCAGTAAAATACAATCCAAAACGC‑AAT‑CACACCTTATTGAGATTCCACGACCGCCGACCACGGCGCGGGCCACCTGTGCGGGGCCCTCGCAGCGCAGAACTGTGCAGAt                                                                                                  <  6:51989/135‑1 (MQ=255)
                                        taGGCCGTTGACCTGTTGCTGCTGTAAAAT‑CCCTACCTAACGCAAAT‑CTCACATTTTTGATATTCCACTAGCGCCCACCAGGGGGCGGGCCACCTGTGCGGGGCCCTCGCAGCGCAGAACTGTGCAGATCTACGa                                                                                            <  6:340220/134‑1 (MQ=255)
                                        aaGGCCTTTCACCTGTTGTTCCATAAAAATACCATCAATACCCC‑AAT‑CTCCCGTGTTTTCTATTGCAGCCGCGCGACCCACGGCGCGGGCCACCTGTGCGGGGCCCTCGCAGCGCAGAACTGTGCAGATCTACGa                                                                                            <  2:248704/134‑1 (MQ=255)
                                          ggCCGTGGACCCTTTGCTGCAGGAAAATACCATCCATACCGC‑AAG‑CTCACCTGATTGCGATTGCACGAGCGCCGACCCCGGCGCGGGCCACCTGTGCGGGGCCCTCGCAGCGCAGAACTGTGCAGATCTACGAAt                                                                                          <  6:77577/135‑1 (MQ=255)
                                            ccGTTGACCAGTTGCTGCAGGAAAATACCATCCATAACGC‑AAT‑CTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGGCACCAGTTCGGGGCCCTCGCAGAGCCCAACATTGTAGATCTACGAATaa                                                                                        >  8:7490/1‑135 (MQ=255)
                                             cGTTGACCAGTTGCTGCAGGAAAATACCATCCATAACGC‑AAT‑CTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCCCCGGTGCGGGGGCCGAGCACCGCCGAACGGTTGAGCTCTTCGAATa                                                                                         >  8:24814/1‑133 (MQ=255)
                                                   ccAGTTGCTGCAGGAAAATACCATCCATAACGC‑AAT‑CTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGGCCCCTGTGCGGGGCCCTCGCCCCGCAGAACTGTGCAGATGTACGAATAAATACAg                                                                                  >  4:79228/1‑134 (MQ=255)
                                                        tgctgcAGTAAAATACCAACCAAAACGC‑AAT‑CCCACCTGATTGAGATTGCACGCGCCCCGACCACGGCGCGGGCCACCTGTGCGGGGCCCTCGCAGCGCAGAACTGTGCAGATCTACGAATAAAGACAGGTAcc                                                                             <  4:86201/134‑1 (MQ=255)
                                                         gctgcAGGAAAATACCATCCATAACGC‑AAT‑CTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCACCTGTTCGGGGCCCTCGCAGCGCAGAACGGTGCAGATCTACGAATAAAGACAGGTACcgc                                                                           >  4:180071/1‑135 (MQ=255)
                                                              aGGAACATACCATCCATAACGC‑AAT‑CTCCCCTGATTGAGATTGCACGCGCGCCGCCCCCGGCGCGGGCCACCTGTGCGGGGCCCTCGCAGCGCAGAACTGTGCAGATCTACGAATAAAGACAGGTACCGCggg                                                                        <  6:335582/133‑1 (MQ=255)
                                                                             atAACGC‑AAT‑CTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGGCACCGGTGCGGGGGCCGCGCAGCGCCGAACTGGGCAGCTCTTCGAATAAAGACAGGTACCGCGGGGCTTGGGCCCGAGCAAcc                                                      >  7:96017/1‑136 (MQ=255)
                                                                             atAACGC‑AAT‑CTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCACCAGTGCGGGGCCCGCGCACCGCCGAACGGTGGAGATCTAACAATAAATAAAGGTAACGCGGGGGGTTGGCCCGAGCCAc                                                       >  6:424749/1‑135 (MQ=255)
                                                                                           cACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGGCCCCGGTGCGGGGGCCGCGCCGCGCCGCACGGTGGAGATCTACGAATAAAGACACGTACCGCGGGGGGTGGGGCCGAGCCACCCCCcgggcccc                                           >  6:334183/1‑133 (MQ=255)
                                                                                           cACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCCCCTGTGCGGGGCCCTCGCAGCGCAGAACTGTGCAGCTGTACGAATTAAGACAGGTACCGGGGGGCGTTGGCCCGAGCAACCGCCccggccgcg                                          >  7:325496/1‑136 (MQ=255)
                                                                                           cACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCACCAGTGGGGGGCCCTCGCCGCGCCGAACTGTGCAGATCTTCGTATTAAGACAGGTTCCGCCGGGCGTGGGCCCGCGCAACCGCccgcgccgcg                                          >  8:333782/1‑136 (MQ=255)
                                                                                              cTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCCCCCGTGCCGGGCCCGCGCCGCGCCGAACAGTGCAGAGCAACGCAAAAACACAGGTACCGCGGGGCGTGGGCCCCCGCCACCCCCccggccgcggg                                        >  3:384744/1‑135 (MQ=255)
                                                                                                 aTTGCGATGGCACGCGCGCCGCCCCCGGCGCGGGCCACCTGTGCGGGGCCCGCGCAGCGCAGAACTGTGCAGATCTACGAATAAAGACAGGTACCGCGGGGCTTGGGCCCGAGCAACCGCCcgggccgcgggccg                                     <  3:120833/135‑1 (MQ=255)
                                                                                                     agaTTGCACGAGCGCCGACCACGGCGCGGGCCCCCTGTGCGGGGCCCTCGCCGCGCAGCACTGTGCAGATCTACGAATAAAGACAGGTACCGCGGGGCTTGGGCCCGAGCACCCGCCccggccgcgggccgggggt                                >  7:355687/1‑131 (MQ=255)
                                                                                                              gAGCGCCGACCACGGCGCGGGGCCCCGGTGCGGGGCCCGCGCAGCGCAGGCCTGTGGAGATAGACTTATAAAGACAGGTTCCGCGGGGCGTGGGCCCGCGCCACACCCCGGGCCGCGGGGCCTGGTGTTTTTTGt                        >  6:295113/1‑135 (MQ=255)
                                                                                                                   ccGACCACGGCGCGGGCCACCTGTGCGGGGCCCTCGCAGCGCAGAACTGTGCAGATCTACGAATAAAGACAGGTACCGCGGGGCTTGGGCCCGAGCAACCGCCCGGGCCGCGGGCCGTTGTTATTTTTGTTTTTcc                  <  7:24814/136‑1 (MQ=255)
                                                                                                                       ccccGGCGCGGGCCACCTGGGCGGGGCCCTCGCAGCGCAGAACTGTGCAGATCTACGAATAAAGCCCGGTACCGCGGGGCTTGGGCCCGAGCAACCGCCCGGGCCGCGGGCCGTTGTTATTTTTGTTTTTCCAGc               <  2:351984/132‑1 (MQ=255)
                                                                                                                                gggCCACCTGTGCGGGGCCCTCGCAGCGCAGAACTGTGCAGATCTACGAATAAAGACAGGTACCGCGGGGCGTGGGCCCGAGCAACCGCCCGGGCCGCGGGCCGTTGTTATTTTTGTTTTTCCAGCTGGTGGTAtt     <  8:225677/136‑1 (MQ=255)
                                                                                                                                 ggCCCCCTGGGCGGGGGCCTCGCAGCGCAGAACTGTTCAGATCTACGAATAAAGACAGGTACCGCGGGGCTTGGGCCCGAGCAACCGCCCGGGCCGCGGGCCGTTGTTATTTTTGTTTTTCCAGCTGGTGGTAtt     <  3:32583/135‑1 (MQ=255)
                                                                                                                                 ggCCAGCTGTGCGGGGCCCTCGCAGCGCAGAACTGTGCAGATCTACGAATAAAGACAGGGTCGGCGGGGCTTGGGCCCGAGCAACCGCCCGGGCCGCGGGCCGTTGTTATTTTTGTTTTTCCAGCTGGTGGTAtt     <  6:131538/135‑1 (MQ=255)
                                                                                                                                 ggCCACCTGTGCGGGGCCCTCGCAGCGCAGAACTGTGCAGATCTACGAATAAATAAAGGTACCGCGGGGCTTGGGCCCGTGCAACCGCCCGGGCCGCGGGCCGTTGTTATTTTTGTTTTTCCAGCTGGTGGTAtt     <  5:5830/135‑1 (MQ=255)
                                                                                                                                 ggCCACCTGTGCGGGGCCCTCGCAGCGCAGAACTGTGCAGATCTACGAATAAAGACAGGTACCGCGGGGCTTGGGCCCGAGCAACCGCCCGGGCCGCGGGCCGTTGTTATTTTTGTTTTTCCAGCTGGTGGTAtt     <  1:231270/135‑1 (MQ=255)
                                                                                                                                 ggCCACCTGTGCGGGGCCCTCGCAGCGCAGAACTGTGCAGATCTACGAATAAAGACAGGTACCGCGGGGCTTGGGCCCGAGCAACCGCCCGGGCCGCGGGCCGTTGTTATTTTTGTTTTTCCAGCTGGTGGTAtt     <  6:176023/135‑1 (MQ=255)
                                                                                                                                 ggCCACCTGGCCGGGGCCCTCGCAGCGCAGAACTGTGCAGATCTACGAATAAAGCGAGGTACCGCGGGGCTTGGGCCCGAGCAACGGCCCGGGCCGCGGGCCGTGGTTATTTTTGTTTTTCCAGCTGGTGGTAtt     <  6:443401/135‑1 (MQ=255)
                                                                                                                                    cACCTGTGCGGGGCCCTCGCAGCGCAGAACTGTGCAGATCTACGAATAAAGACAGGTACCGCGGGGCTTGGGCCCGAGCAACCGCCCGGGCCGCGGGCCGTTGTTATTTTTGTTTTTCCAGCTGGTGGTATTTAc  <  4:384744/135‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                     
TGTAGCCGATGGCGATCAGGCCATAGACCGACCCGAGGGTCAGGCCGTTGACCAGTTGCTGCAGGAAAATACCATCCATAACGC‑AAT‑CTCACCTGATTGAGATTGCACGAGCGCCGACCACGGCGCGGGCCACCTGTGCGGGGCCCTCGCAGCGCAGAACTGTGCAGATCTACGAATAAAGACAGGTACCGCGGGGCTTGGGCCCGAGCAACCGCCCGGGCCGCGGGCCGTTGTTATTTTTGTTTTTCCAGCTGGTGGTATTTAC  >  ACB122/5525752‑5526016

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