Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606476,7460TG36.5% 166.6 / 18.7 85Q6H (CAA→CACfsrpredicted fosmidomycin efflux system
Reads supporting (aligned to +/- strand):  ref base T (17/37);  new base G (31/0);  total (48/37)
Fisher's exact test for biased strand distribution p-value = 4.34e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 138 total aligned reads displayed.

TGATAGCACCTAAAATACCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTAC  >  REL606/476673‑476821
                                                                          |                                                                           
tGATAGCACCTAAAATACCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGggtg                                                                           >  2:4698315/1‑76 (MQ=255)
tGATAGCACCTAAAATACCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGggtg                                                                           >  2:5466864/1‑76 (MQ=255)
tGATAGCACCTAAAATACCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGggtg                                                                           >  1:4414734/1‑76 (MQ=255)
tGATAGCACCTAAAATACCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGggtg                                                                           >  1:6384917/1‑76 (MQ=255)
 gATAGCACCTAAAATACCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGggtgt                                                                          >  2:7074833/1‑76 (MQ=255)
 gATAGCACCTAAAATACCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGggtgt                                                                          >  2:779547/1‑76 (MQ=255)
 gATAGCACCTAAAATACCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGggtgt                                                                          >  1:5139708/1‑76 (MQ=255)
 gATAGCACCTAAAATACCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGggtgt                                                                          >  1:6498427/1‑76 (MQ=255)
   aaGCCCCCAAAATTCCAAAACGATGTCGGGGCGTGGGTTGTCGACGCCCCCCCCCCCGCCACAGGCTGGGgttgtt                                                                         <  1:4552148/75‑1 (MQ=255)
       aCCTAAAATACCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCAcac                                                                    >  2:6418369/1‑74 (MQ=255)
       aCCTAAAATACCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTc                                                                    >  2:467761/1‑76 (MQ=255)
         cTAAAATACCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGGGTTCACTCAt                                                                  >  2:1129946/1‑76 (MQ=255)
          tAAAATACCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCCCTCAtt                                                                 >  1:4974463/1‑76 (MQ=255)
           aaaaTACCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGGTTACTCATTg                                                                >  2:7233023/1‑76 (MQ=255)
            aaaTACCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCATTGc                                                               >  2:5993780/1‑76 (MQ=255)
                aCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCATTGGCATa                                                           >  1:5491104/1‑76 (MQ=255)
                aCCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCATTGCCCTa                                                           >  1:6682457/1‑76 (MQ=255)
                  cAAGCGTTGTTGGGGCCTGGGTCGTCCCCCCCCCCCCCCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAgg                                                         <  1:1128308/76‑1 (MQ=255)
                  cAAGCGTTGTGCGGCCTTTGGTCGTTGCCCCCCCCCCCCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAgg                                                         <  2:7376893/76‑1 (MQ=255)
                    aaaCGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCTCTCATTGCCCTTGGa                                                        >  1:4828020/1‑76 (MQ=255)
                    aaaCGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTTATTGCCATAGGa                                                        >  1:4925403/1‑76 (MQ=255)
                    aaaCGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCCTTGCCATAggg                                                        >  1:34371/1‑75 (MQ=255)
                    aaaCGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCATTGCCATAGGa                                                        >  2:5574568/1‑76 (MQ=255)
                    aaaCGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGGTTACTCATTGCCCTAGGa                                                        >  1:6317681/1‑76 (MQ=255)
                      aCGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTTACTCCTTGCCATAGGaaa                                                      >  1:5256804/1‑76 (MQ=255)
                      aCGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCATTGCCATAGGgaa                                                      >  1:6269899/1‑76 (MQ=255)
                        gATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCATTGCCATAGGGAAGt                                                    >  1:1486753/1‑76 (MQ=255)
                         aTGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGGTTACTCATTGCCATAGGAAAggc                                                   >  1:3372215/1‑74 (MQ=255)
                           gTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCATTGCCATAGGAAAGTCAc                                                 >  2:2053093/1‑76 (MQ=255)
                             tCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCCTTGGCATAGGAAAAGCACtt                                               >  2:2166840/1‑76 (MQ=255)
                               gggCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCCTTGGAAAGATACtttt                                             >  2:7282269/1‑76 (MQ=255)
                               gggCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCATTGCCATAGGAAAGTCACtttt                                             >  1:2925348/1‑76 (MQ=255)
                                ggCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCATTGCCATAGGAAAGTCACttttt                                            >  2:1002622/1‑76 (MQ=255)
                                 gCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCATTGCCATAGGAAAATCACtttttt                                           >  2:2069135/1‑75 (MQ=255)
                                  ccTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCCTTGCCATAGGGAAGTCACTTTTTCa                                          >  1:2544018/1‑76 (MQ=255)
                                   cTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTTACTCATTGGCCCTGGTAAGTTACttttttcg                                         >  2:6234156/1‑73 (MQ=255)
                                   cTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAg                                         >  1:3240606/1‑76 (MQ=255)
                                    ttGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGGCACTTTTTCAgg                                        >  1:5601119/1‑76 (MQ=255)
                                    ttGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAgg                                        >  1:6634879/1‑76 (MQ=255)
                                      ggTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCATTGCCCAAGGAAAGTCACTTTTTCAGGGt                                      >  1:3380193/1‑76 (MQ=255)
                                      ggTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCATTGCCATAGGAAAGTCACTTTTTTAGGGt                                      >  2:4802603/1‑76 (MQ=255)
                                      ggTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCCGGGt                                      >  1:1565680/1‑76 (MQ=255)
                                      ggTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCATTGCCATAAGAAAGTCACTTTTTCAgggg                                      >  2:2409501/1‑75 (MQ=255)
                                       gTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTTACTCATTTCCCTAGGGAAGTCACTTTTTTAggggt                                     >  1:4301952/1‑74 (MQ=255)
                                         cGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGc                                   >  1:3017558/1‑76 (MQ=255)
                                         cGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTTACTCATTGCCATAGGAAAGTCACTTTTTTAGGGGTGc                                   >  2:4062280/1‑76 (MQ=255)
                                             gACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCATTGCCATAGGAAAGTCACCTTTTTAGGGGTGCGATg                               >  1:2664297/1‑76 (MQ=255)
                                             gACGCAGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCATTGCCATAGGAAAGGCACTTTTTCAGGGGTGCGATg                               >  2:1419464/1‑76 (MQ=255)
                                             gACGCAGCCGCGCCCGCCACAGGCTGGGGATGTTCACACATTGCCCTAGGAAAGGCACTTTTTCAGGGTTGCGATg                               >  1:7106460/1‑76 (MQ=255)
                                               cGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTa                             >  1:6178271/1‑76 (MQ=255)
                                                gCAGCCGCGCCCGCCACAGGCTGGGGGTGTTTACTCATTGCCATAGGAAAGTCACTTTTTTAGGGTTGCGATGTaa                            >  2:4608764/1‑76 (MQ=255)
                                                 cAGCCGCGCCCGCCACAGGCTGGGGGTGTTCCCCCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTaaa                           >  1:4534483/1‑76 (MQ=255)
                                                 cAGCCGCGCACGCCACAGGCTGGGGGTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTaaa                           >  2:4184298/1‑76 (MQ=255)
                                                  cgCCCCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAg                          <  2:4937905/75‑1 (MQ=255)
                                                  aGCCGCGCCCGCCACAGGCTGGGGGTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAg                          >  1:1439643/1‑76 (MQ=255)
                                                      ccccccgccACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATg                      <  2:5129202/73‑1 (MQ=255)
                                                       cgcccgccACAGGCTGGGGGTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGa                     >  2:5027097/1‑76 (MQ=255)
                                                       cccccgccACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGa                     <  2:1006697/74‑1 (MQ=255)
                                                       cccccgccACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGa                     <  1:4896182/74‑1 (MQ=255)
                                                        ccccgccACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGAt                    <  2:26837/75‑1 (MQ=255)
                                                        ccccgccACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGAt                    <  2:3019207/75‑1 (MQ=255)
                                                        ccccgccACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGAt                    <  2:3413841/75‑1 (MQ=255)
                                                           cgccACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCtt                 <  1:2797414/76‑1 (MQ=255)
                                                           cgccACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCtt                 <  1:4491962/76‑1 (MQ=255)
                                                           cgccACAGGCTGGGGGTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCtt                 >  2:3193338/1‑76 (MQ=255)
                                                            gccACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTa                <  1:2897613/76‑1 (MQ=255)
                                                             ccACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTAt               <  1:2647261/76‑1 (MQ=255)
                                                             ccACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTAt               <  2:6423772/76‑1 (MQ=255)
                                                             ccACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTAt               <  2:7050560/76‑1 (MQ=255)
                                                             ccACAGGCTGGGGGTGTTCACTCATTGCCATAGGAAAGTCCCTTTTTCAGGGGTGCGATGTAAAGAATGATCTTAt               >  2:321579/1‑76 (MQ=255)
                                                              cacaGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTAtt              <  2:17703/76‑1 (MQ=255)
                                                               acaGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTAttt             <  2:4064088/76‑1 (MQ=255)
                                                               acaGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTAttt             <  2:6497704/76‑1 (MQ=255)
                                                               acaGGCTGGGGGTGTTCACTCATTGCCATAGGAAAGTCACTTTTTTAGGGTTGCGGTGTAAAGAATGATCTTAttt             >  2:4207522/1‑76 (MQ=255)
                                                                  ggCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTtgtg          <  1:6235118/76‑1 (MQ=255)
                                                                   gCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGa         <  2:912617/76‑1 (MQ=255)
                                                                   gCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGa         >  1:3229722/1‑76 (MQ=255)
                                                                   gCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGa         <  1:3951949/76‑1 (MQ=255)
                                                                    cTGGGGTTGTTCACTCATTGCCCTAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGat        >  2:4712041/1‑76 (MQ=255)
                                                                    cTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGat        <  1:1681266/76‑1 (MQ=255)
                                                                    cTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGat        >  1:4097348/1‑76 (MQ=255)
                                                                     tGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGatt       >  1:1411302/1‑76 (MQ=255)
                                                                     tGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGatt       <  1:2578669/76‑1 (MQ=255)
                                                                     tGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGatt       <  2:276443/76‑1 (MQ=255)
                                                                      ggggTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGatta      >  2:3676144/1‑76 (MQ=255)
                                                                       gtgttgttCTCTCATTGCCATATGAAAGTCGCTATTTCAGGGTTGCGATGAAAAGAATGATGTTATTTGTGtttat     <  2:3145916/74‑1 (MQ=255)
                                                                       gggTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGattat     <  2:694694/76‑1 (MQ=255)
                                                                       gggTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGattat     <  2:3145915/76‑1 (MQ=255)
                                                                        ggTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGattatt    <  1:5441428/76‑1 (MQ=255)
                                                                        ggTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGattatt    >  2:6997132/1‑76 (MQ=255)
                                                                        ggTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGattatt    <  2:2793863/76‑1 (MQ=255)
                                                                         gttgttCACTCCTTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGAttatta   >  2:609600/1‑76 (MQ=255)
                                                                         gttgttCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGAttatta   >  2:2207995/1‑76 (MQ=255)
                                                                         gttgttCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGAttatta   <  2:2174883/76‑1 (MQ=255)
                                                                         gttgttCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGAttatta   <  1:551585/76‑1 (MQ=255)
                                                                         gttgttCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGAttatta   <  1:4559928/76‑1 (MQ=255)
                                                                         gttgttCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGAttatta   <  2:3428300/76‑1 (MQ=255)
                                                                         gttgtgCACTCATTGCTATACGAAAGTGACGTTGTCAGGGTTGCGATTTAAATAATGATCTTATTTGTGAttatta   <  2:2174886/76‑1 (MQ=255)
                                                                          ttgttCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTAc  >  2:4057207/1‑76 (MQ=255)
                                                                          ttgttCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTAc  <  2:5741420/76‑1 (MQ=255)
                                                                          |                                                                           
TGATAGCACCTAAAATACCAAACGATGTTCGGGCCTTGGTCGTTGACGCAGCCGCGCCCGCCACAGGCTGGGGTTGTTCACTCATTGCCATAGGAAAGTCACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTAC  >  REL606/476673‑476821

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