Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,702,2450AC34.6% 49.5 / 10.3 26T5P (ACC→CCC) addadenosine deaminase
Reads supporting (aligned to +/- strand):  ref base A (16/1);  new base C (0/9);  total (16/10)
Fisher's exact test for biased strand distribution p-value = 3.20e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CTTTGTTTTTATAATGGTGCGCACTTTTATATCCAGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGA  >  NC_000913/1702182‑1702309
                                                               |                                                                
cTTTGTTTTTATAATGGTGCGCACTTTTATATCCAGAAAAAGAGTGCGACCATGATTGATaccac                                                                 <  1:3064279/65‑1 (MQ=255)
 tttGTTTTTATAATGGTGCGCACTTTTATATCCAGAAAAAGAGTGCGACCATGATTGATaccacc                                                                <  1:1843612/65‑1 (MQ=255)
    gTTTTTATAATGGTGCGCACTTTTATATCCAGAAAAAGAGTGCGACCATGATTGATACCACCCTg                                                             >  1:2114533/1‑65 (MQ=255)
     tttttATAATGGTGCGCATTTTTATATCCAGAAAAAGAGTGCGACCATGATTGATCCCCCCCTGc                                                            <  1:1181634/65‑1 (MQ=255)
     tttttATAATGGTGCGCACTTTTATATCCAAAAAAAGAGTGCGCCCATGATTGATACCCCCCTGc                                                            <  1:3409916/65‑1 (MQ=255)
      ttttATAATGGTGCGCACTTTTATATCCAGAAAAAGAGTGCGACCATGATTGATACTACCCTGcc                                                           >  1:712585/1‑65 (MQ=255)
             aTGGTGGGCTCTTTTATATCAAGAAAAAGAGTGCGACCATGATTGATCCCCCCCTGCCATTAACt                                                    <  1:2049107/65‑1 (MQ=255)
             aTGGTGCGCATTTTGATATCCAGAAAAAGAGTGCGACCATGATTGATCCCCCCCTGCCATTAACt                                                    <  1:1364078/65‑1 (MQ=255)
             aTGGTGCGCACTTTTATCTCCAGAAAAAGAGTGCGCCCATGATTGATACCCCCCTGCCATTAACt                                                    <  1:1271511/65‑1 (MQ=255)
             aTGGTGCGCACTTTTATATCAAAAAAAAGAGTGCGACCATGATTGATCCCCCCCTGCCATGAACt                                                    <  1:2094670/65‑1 (MQ=255)
               ggTGCGCACTTTTATATCCAGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGa                                                  >  1:3519289/1‑65 (MQ=255)
               ggTGCGCACTTTTATATCCAGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGa                                                  >  1:2998960/1‑65 (MQ=255)
               ggTGCGCACTTTTATATCCAGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGa                                                  >  1:1243082/1‑65 (MQ=255)
                gggcgcACTTTTATATCCAGAAAAGGGGTGGGGCCATGATTGATCCCCCCCTGCCATTAACTGat                                                 <  1:1403179/63‑1 (MQ=255)
                  gcgcATTTTTATATCCAGAAAAAGAGTGCGACCATAATTGATACCCCCCTGCCATTAACTGatat                                               <  1:2664703/65‑1 (MQ=255)
                       cTTTTATATCCAGAAAAAGAGTGCGCCCATGATTGATCCCCCCCTGCCATTAACTGATATCCATc                                          <  1:2162924/65‑1 (MQ=255)
                            atatCCAGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCAc                                     >  1:2175630/1‑65 (MQ=255)
                             tatCAAAAAAAAGGTTGGGACCGTGATTGATCCCCCCCTGCCATTAACTGATATCCATCGCCAcc                                    <  1:2616963/65‑1 (MQ=255)
                              atCCAGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCt                                   >  1:2258439/1‑65 (MQ=255)
                                ccAGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTg                                 >  1:2004779/1‑65 (MQ=255)
                                ccAGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTg                                 >  1:2231679/1‑65 (MQ=255)
                                    aaaaaGAGTGCGACCATGATTAACACCGCCCTGCCATTAACTGATATCCATCGCCACCTTGATgg                             <  1:1607650/65‑1 (MQ=255)
                                      aaaGAGTGCGCCCATGATTGACCCCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCa                           <  1:2654592/65‑1 (MQ=255)
                                        agagTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAAc                         >  1:732429/1‑65 (MQ=255)
                                        agagTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAAc                         <  1:266629/65‑1 (MQ=255)
                                         gagTGCGACCATGATTGATCCCCCCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACa                        <  1:1577982/65‑1 (MQ=255)
                                         gagTGCGACCATGATTGATACCCCCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACa                        <  1:707841/65‑1 (MQ=255)
                                          agTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACAt                       >  1:3909330/1‑65 (MQ=255)
                                          agTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACAt                       >  1:204881/1‑65 (MQ=255)
                                           ttGCGACCATTATTGATCCCCCCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACAtt                      <  1:4047293/64‑1 (MQ=255)
                                                aCCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTcc                 >  1:2201018/1‑65 (MQ=255)
                                                     gattgatCCCCCCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGa            <  1:975168/65‑1 (MQ=255)
                                                       ttgatACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGAcc          >  1:2021015/1‑65 (MQ=255)
                                                         gatACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCAt        >  1:581658/1‑65 (MQ=255)
                                                          atcccaccCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCAtt       <  1:3068100/62‑1 (MQ=255)
                                                           tACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTc      >  1:41338/1‑65 (MQ=255)
                                                            accaccCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCt     >  1:1149687/1‑65 (MQ=255)
                                                            accaccCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCt     >  1:3658837/1‑65 (MQ=255)
                                                            accaccCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCt     >  1:3529440/1‑65 (MQ=255)
                                                            accaccCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCt     >  1:3504092/1‑65 (MQ=255)
                                                            accaccCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCt     >  1:1185754/1‑65 (MQ=255)
                                                             ccaccCTGCCATTAACTGATATTCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCtc    >  1:611807/1‑64 (MQ=255)
                                                             ccaccCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCtt    >  1:2507766/1‑65 (MQ=255)
                                                             ccaccCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCtt    >  1:135131/1‑65 (MQ=255)
                                                               accCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGa  >  1:2415366/1‑65 (MQ=255)
                                                               accCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGa  >  1:1092218/1‑65 (MQ=255)
                                                               |                                                                
CTTTGTTTTTATAATGGTGCGCACTTTTATATCCAGAAAAAGAGTGCGACCATGATTGATACCACCCTGCCATTAACTGATATCCATCGCCACCTTGATGGCAACATTCGTCCCCAGACCATTCTTGA  >  NC_000913/1702182‑1702309

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: