Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913303,6310AG22.9% 63.1 / 15.0 35intergenic (+26/+222)yagU/ykgJDUF1440 family inner membrane acid resistance protein/UPF0153 cysteine cluster protein
Reads supporting (aligned to +/- strand):  ref base A (13/14);  new base G (4/4);  total (17/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 = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGGACATTTAGTCTGGTTCTGGTCTATTGAAATTATTCGCAGAGATTTACGAAACAGAATTACTCATGAGCCAGACCCTGAGATCCCTTTAGGCTCAAACAGATAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCAGGATGATGTAGATCCGAAGAACGCTACAGAAGAGAGGCATTGTTGCTGGCAAATAGAAGAA  >  NC_000913/303500‑303774
                                                                                                                                   |                                                                                                                                               
tttgccagcaacaatgcctctcttctgtagcgttcttcggatctacatcatcctgaagtagcgcgatatactcactgactatctgcatcagaatataaaaagcaatgtttttaacctATAAAAATGGCGCAAATACAGCGCCAttt                                                                                                                                   >  4:106926‑M1/118‑146 (MQ=255)
           tctggttctggtCTATTGAAATTATTCGCAGAGATTTACGAAACAGAATTACTCATGAGCCAGACCCTGAGATCCCTTTAGGCTCAAACAGATAATGCATTGAATGATAAAAATGGCGCTGTATTTgcgccatttttatcattcaatg                                                                                                                      <  6:32381‑M2/148‑23 (MQ=255)
           tctggttctggtCTATTGAAATTATTCGCAGAGATTTACGAAACAGAATTACTCATGAGCCAGACCCTGAGATCCCTTTAGGCTCAAACAGATAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTaaaaa                                                                                                                     <  7:128092/149‑1 (MQ=255)
           gctggttctcgtCTATTGAAATTATTCGCAGAGATTTACGAAACAGAATTACTCAAGAGCCAGACCCTGAGATCCCTTTAGGCTCAAACAGATAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTaaaaa                                                                                                                     <  7:126269/148‑1 (MQ=255)
                        tATTGAAATTATTCGCAGAGATTTACGAAACAGAATTACTCATGAGCCAGACCCTGAGATCCCTTTAGGCTCAAACAGATAATGCATTGAATGATAAAAATGGCGCAAATACAgcgccattttta                                                                                                                                >  1:141655‑M2/1‑113 (MQ=255)
                        tATTGAAATTATTCGCAGAGATTTACGAAACAGAATTACTCATGAGCCAGACCCTGAGATCCCTTTAGGCTCAAACAGATAATGCATTGAATGATAAAAATGGCGCAAATACAgcgccattttta                                                                                                                                <  2:141655‑M2/125‑13 (MQ=255)
                             aaaTTATTCGCAGAGATTTACGAAACAGAATTACTCATGAGCCAGACCCTGAGATCCCTTTAGGCTCAAACAGATAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCt                                                                                                   >  5:708/1‑149 (MQ=255)
                               attattCGCAGAGATTTACGAAACAGAATTACTCATGAGCCAGACCCTGAGATCCCTTTAGGCTCAAACAGATAATGCATTGAATGATAAAAATGGCGCTGTATTTgcgccatttttatcattcaatgcattatctgtttgagcctaaa                                                                                                 >  5:32381‑M2/1‑106 (MQ=255)
                               attattCGCAGAGATTTACGAAACAGAATTACTCATGAGCCAGACCCTGAGATCCCTTTAGGCTCAAACAGATAATGCATTGAATGATAAAAATGGCGCTGTATTTgcgccatttttatcattcaatgcattatctgtttgagccta                                                                                                   <  7:46879‑M2/147‑42 (MQ=255)
                                                 cGAAACAGAATTACTCATGAGCCAGACCCTGAGATCCCTTTAGGCTCAAACAGATAATGCATTGAATGATAAAAATGGCGCTGTATTTgcgccatttttatcattcaatgcattatctgtttgagcctaaagggatctcagggtctggc                                                                               >  8:46879‑M2/1‑88 (MQ=255)
                                                  gAAACAGAATTACTCATGAGCCAGACCCTGAGATCCCTTTAGGCTCAAACAGATAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGta                                                                              >  3:60160/1‑149 (MQ=255)
                                                         aaTTACTCATGAGCCAGACCCTGAGATCCCTTTAGGCTCAAACAGATAATGCATTGAATGATAAAAATGGCGCAAATACAgcgccattttt                                                                                                                                 <  7:34396‑M2/91‑12 (MQ=255)
                                                         aaTTACTCATGAGCCAGACCCTGAGATCCCTTTAGGCTCAAACAGATAATGCATTGAATGATAAAAATGGCGCAAATACAgcgccattttt                                                                                                                                 >  8:34396‑M2/1‑80 (MQ=255)
                                                              cTCATGAGCCAGACCCTGAGATCCCTTTAGGCTCAAACAGATAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTAc                                                                   <  6:708/148‑1 (MQ=255)
                                                              cTCATGAGCCAGACCCTGAGATCCCTTTAGGCTCAAACAGATAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACt                                                                  >  3:28946/1‑149 (MQ=255)
                                                                 aTGAGCCAGACCCTGAGATCCCTTTAGGCTCAAACAGATAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCa                                                               <  4:60160/149‑1 (MQ=255)
                                                                        aGACCCTGAGATCCCTTTAGGCTCAAACAGATAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCAGgatgat                                                        <  2:168275/149‑1 (MQ=255)
                                                                        aGACCCTGAGATCCCTTTAGGCTCAAACAGATAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCAGgatgat                                                        >  1:168275/1‑149 (MQ=255)
                                                                                      cTTTAGGCTCAAACAGATAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCAGGATGATGTAGATCCgaagaa                                          <  7:23144/149‑1 (MQ=255)
                                                                                              tCAAACAGATAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGt                                                                               <  2:37096/104‑1 (MQ=255)
                                                                                              tCAAACAGATAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGt                                                                               >  1:37096/1‑104 (MQ=255)
                                                                                              tCAAACAGATAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCAg                                                              <  8:32665/121‑1 (MQ=255)
                                                                                              tCAAACAGATAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCAg                                                              >  7:32665/1‑121 (MQ=255)
                                                                                                aaaaagcaatgtttttaacctATAAAAATGGCGCTGTATTTGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCAGgatga                                                         <  3:84202‑M1/103‑1 (MQ=255)
                                                                                                aaaaagcaatgtttttaacctATAAAAATGGCGCTGTATTTGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCAGgatga                                                         >  4:84202‑M1/22‑124 (MQ=255)
                                                                                                     gATAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCAGGATGATGTAGATCCGAAGAACGc                                       >  3:8749/1‑137 (MQ=255)
                                                                                                     gATAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCAGGATGATGTAGATCCGAAGAACGc                                       <  4:8749/137‑1 (MQ=255)
                                                                                                       tAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCAGGATGATGt                                                      <  1:125760/120‑1 (MQ=255)
                                                                                                       tAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCAGGATGATGt                                                      >  2:125760/1‑120 (MQ=255)
                                                                                                        atgtttttaacctATAAAAATGGCGCTGTATTTGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCa                                                               <  7:81588‑M1/97‑1 (MQ=255)
                                                                                                        atgtttttaacctATAAAAATGGCGCTGTATTTGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCa                                                               >  8:81588‑M1/14‑110 (MQ=255)
                                                                                                              ttGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCAGGATGATGTAGATCCGAAGAACGCTACAGAAGAGAGGCAttgttg                  >  8:128979/1‑149 (MQ=255)
                                                                                                                           aTGGCGCAAATACAgcgcaatttttatcattaaatgcattatctgtttgagcctaaagggatctcagggtctatctcatgagtaattctgtttagtaaatatctgagaataatatcaatagaacagaac                         <  1:174576‑M2/129‑116 (MQ=255)
                                                                                                                            TGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCAGGATGATGTAGATCCGAAGAACGCTACAGAAGAGAGGCATTGTTGCTGGCAAATag       <  2:123008‑M1/146‑1 (MQ=255)
                                                                                                                              GCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCAGGATGATGTAGATCCGAAGAACGCTACAGAAGAGAGGCATTGTTGCTGGCAAATAgaagaa  >  1:123008‑M1/1‑149 (MQ=255)
                                                                                                                                   |                                                                                                                                               
TGGACATTTAGTCTGGTTCTGGTCTATTGAAATTATTCGCAGAGATTTACGAAACAGAATTACTCATGAGCCAGACCCTGAGATCCCTTTAGGCTCAAACAGATAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCAGGATGATGTAGATCCGAAGAACGCTACAGAAGAGAGGCATTGTTGCTGGCAAATAGAAGAA  >  NC_000913/303500‑303774

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