Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,576,0770CA45.8% 22.3 / 14.9 24D515Y (GAT→TAT) yddBputative TonB‑dependent outer membrane receptor
Reads supporting (aligned to +/- strand):  ref base C (9/4);  new base A (0/11);  total (9/15)
Fisher's exact test for biased strand distribution p-value = 5.89e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

TCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATACCGCACGCCGGGC  >  NC_000913/1575942‑1576219
                                                                                                                                       |                                                                                                                                              
tCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGc                                                                                                                                  <  1:646605/150‑1 (MQ=255)
                ctgcAACCCCATTGCCATTTCATCTTTATAAGGTTTTTTAAAATCCTAAAAACGCGTCATTTTTTTATTCCCTAACCCCAATTCCCCCCAGCTATTGCGAATTTCACGTAACCCCATAACAAGAATATTCCCGCCATAGTAACGGtt                                                                                                                     <  1:149859/147‑1 (MQ=255)
                   cAATCCCATTGCCAGTTCATCTTAATAAGGCTTTTCCAAACCCTGATATCGCGTTATATTTTTATTCCCTGATCCCGATCCCCTCCATCTATGGCGGATATCCCATAACCCAAAATAAAAAATTTTCCCGCCATAGTAACGGt                                                                                                                      <  1:715872/143‑1 (MQ=255)
                           ttccaaGTTCTTCTTTAAAAGGTTTTTAAAAATCCTAATACCGCGTCAGATTTTTTTTCCCTAATACCTATTCCGTCCATTTTTTGCGTATATCATATAACCCAATAAAAAAAATATCCCCGCCATAGTAACGGTTATAACCTGCTGTaa                                                                                                       <  2:819152/145‑1 (MQ=255)
                           ttGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAAt                                                                                                      >  2:630043/1‑151 (MQ=255)
                                        gTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATGTCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTg                                                                                          >  1:789724/1‑150 (MQ=255)
                                             agggttttttCAAATCTTGAAACCGCGTAAGTTTTTTTTTCCCTAATACCGATTCCGTCCAGTTATTGGGGATATCACTTAACCCAATATAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGc                                                                                    <  2:666031/145‑1 (MQ=255)
                                             aaggttttttCAAATCCTAAAAACGCGTCAAATTTTTTTTCCCTAATCCCTATCCCTCCCAGCTATTGCGTATATCACATACCCCAATAAAAAAATTTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTa                                                                                      <  1:16820/143‑1 (MQ=255)
                                                   ttttAAAATCTTAATACCGCTTCAGATTTTTATTCCCGGATCCCGATCCCGCCCAGCTATTGGGAATACCACGTACTCCCAAAAAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCaaaa                                                                                <  2:292919/149‑1 (MQ=255)
                                                          aTCCTGATAATGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCa                                                                        >  1:478542/1‑150 (MQ=255)
                                                                aTACCGCGTCAGATTTTTTTTCCCTACTACCTATTCCGTCCAGTTATTGAGAATTTCACGTATTCCAATATAAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAata                                                                             <  1:257926/139‑1 (MQ=255)
                                                                  aaCGCTTCAGATTTTTATTCCCTGATCCCGATCCCCTCCAGCTTTTGCGGATATCACTTACTCCCATAAAAAGAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTc                                                                 <  1:172889/149‑1 (MQ=255)
                                                                  aaCGCGTAAGATTTTTTTTCCCTAATCCCTATTCCGTCCAGTTATGGCGGATATCACGTAACCCAATATAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCa                                                                <  1:36616/150‑1 (MQ=255)
                                                                          agagTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGcgcg                                                       >  2:463380/1‑151 (MQ=255)
                                                                              ttttATTCCCTGATACCGATCCCGTCCAGCTATTGGGAATATAACAAAACCCCATAAAAAAAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCgg                                                      <  1:457964/148‑1 (MQ=255)
                                                                               tttATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCgggg                                                    >  1:752784/1‑149 (MQ=255)
                                                                                        ttatccctatcccGCCCAGCTTTGGGGGATATCACAAAACCCAAAAACAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATAtt                                             <  1:405594/136‑1 (MQ=255)
                                                                                             ccGTTTCCGTCCAGCTATGGCGAATATCACGTAACCCAATATAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGttt                                      <  2:542165/149‑1 (MQ=255)
                                                                                                     gTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGt                            >  1:246593/1‑151 (MQ=255)
                                                                                                       ccAGTTATTGCGGAATTCACGAACACCAAAAAAAAAAATTTCCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAg                             <  1:62641/148‑1 (MQ=255)
                                                                                                        cAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTgt                         >  1:715691/1‑151 (MQ=255)
                                                                                                           cTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTgtcat                      >  2:821016/1‑151 (MQ=255)
                                                                                                                         aCGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATACCGCACg        >  1:774828/1‑151 (MQ=255)
                                                                                                                                cccATATCAAAAATTTCCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATACCGCACGCCGGGc  <  2:434219/150‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                              
TCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATACCGCACGCCGGGC  >  NC_000913/1575942‑1576219

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