Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913303,6320AT22.9% 63.9 / 11.4 35intergenic (+27/+221)yagU/ykgJDUF1440 family inner membrane acid resistance protein/UPF0153 cysteine cluster protein
Reads supporting (aligned to +/- strand):  ref base A (13/14);  new base T (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.

TGGACATTTAGTCTGGTTCTGGTCTATTGAAATTATTCGCAGAGATTTACGAAACAGAATTACTCATGAGCCAGACCCTGAGATCCCTTTAGGCTCAAACAGATAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCAGGATGATGTAGATCCGAAGAACGCTACAGAAGAGAGGCATTGTTGCTGGCAAATAGAAGAAGTATCG  >  NC_000913/303500‑303780
                                                                                                                                    |                                                                                                                                                    
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                                                                    >  7:32665/1‑121 (MQ=255)
                                                                                              tCAAACAGATAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCAg                                                                    <  8:32665/121‑1 (MQ=255)
                                                                                                aaaaagcaatgtttttaacctATAAAAATGGCGCTGTATTTGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCAGgatga                                                               >  4:84202‑M1/22‑124 (MQ=255)
                                                                                                aaaaagcaatgtttttaacctATAAAAATGGCGCTGTATTTGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCAGgatga                                                               <  3:84202‑M1/103‑1 (MQ=255)
                                                                                                     gATAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCAGGATGATGTAGATCCGAAGAACGc                                             <  4:8749/137‑1 (MQ=255)
                                                                                                     gATAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCAGGATGATGTAGATCCGAAGAACGc                                             >  3:8749/1‑137 (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)
                                                                                                                                    aTACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCAGGATGATGTAGATCCGAAGAACTCTACAGAAGAGAGGCATTGTTGCTGGCAAATAGAAGAAGTATCg  <  7:128979/149‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                                    
TGGACATTTAGTCTGGTTCTGGTCTATTGAAATTATTCGCAGAGATTTACGAAACAGAATTACTCATGAGCCAGACCCTGAGATCCCTTTAGGCTCAAACAGATAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCAGGATGATGTAGATCCGAAGAACGCTACAGAAGAGAGGCATTGTTGCTGGCAAATAGAAGAAGTATCG  >  NC_000913/303500‑303780

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: