Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913272,3301.G20.0% 126.1 / 45.8 50intergenic (+141/‑250)insN/eyeApseudogene, IS911 transposase A;IS, phage, Tn; Transposon‑related functions; extrachromosomal; transposon related/novel sRNA, function unknown, CP4‑6 propahge
Reads supporting (aligned to +/- strand):  ref base . (19/21);  new base G (5/5);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGTTGAAAACGTCACCA  >  NC_000913/272184‑272473
                                                                                                                                                   |                                                                                                                                               
aCCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCg                                                                                                                                                <  2:38325/148‑1 (MQ=255)
  cTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCa                                                                                                                                             >  7:104108/1‑149 (MQ=255)
       gtgtGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGGCTGTCTGGAATAGTGAACCGGTAACTCACGCTGAGAACCCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGcc                                                                                                                                         <  3:91138/148‑1 (MQ=255)
               ggAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGGGGAGAACCTGAGCGGCAGTGACCTGCGGGCATGCCGCAggtca                                                                                                                                <  4:51095/149‑4 (MQ=255)
                gAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGGCATGCCGCAggtcac                                                                                                                               >  4:139363/1‑145 (MQ=255)
                 aGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGGCATGCCGCAggtcact                                                                                                                              <  3:139363/149‑6 (MQ=255)
                 aGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGAt                                                                                                                              <  8:91925/149‑1 (MQ=255)
                   ttCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGGCATGccgcag                                                                                                                                     <  6:18723/141‑5 (MQ=255)
                   ttCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGGCATGccgcag                                                                                                                                     >  5:18723/1‑137 (MQ=255)
                        ggcgggcATGGATGCTTAAATGAACCGCAAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAAcc                                                                                                                        <  6:29013/148‑1 (MQ=255)
                                   aTGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGGCATGccgcag                                                                                                                                     >  5:145629/1‑121 (MQ=255)
                                   aTGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGGCATGccgcag                                                                                                                                     <  6:145629/125‑5 (MQ=255)
                                     gCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGGCATGCCGCAggtca                                                                                                                                >  3:60456/1‑124 (MQ=255)
                                     gCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGGCATGCCGCAggtca                                                                                                                                <  4:60456/127‑4 (MQ=255)
                                      cTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGc                                                                                                                     >  6:24449/1‑137 (MQ=255)
                                      cTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGc                                                                                                                     <  5:24449/137‑1 (MQ=255)
                                       ttAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTCATGTAACCCGCTGACACAGGGGa                                                                                                         >  6:46196/1‑148 (MQ=255)
                                                 cgcgAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCgaga                                                                                              <  6:91572/149‑1 (MQ=255)
                                                  gcgAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGAt                                                                                             <  7:109846/149‑1 (MQ=255)
                                                           tctgGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGcc                                                                                    >  1:131224/1‑149 (MQ=255)
                                                           tctgGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGcc                                                                                    >  4:32001/1‑149 (MQ=255)
                                                             tgGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATc                                                                                            >  8:133603/1‑139 (MQ=255)
                                                             tgGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATc                                                                                            <  7:133603/139‑1 (MQ=255)
                                                               gAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCgaga                                                                                <  2:131224/149‑1 (MQ=255)
                                                                  tatTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCga                                                                                                <  6:12070/130‑1 (MQ=255)
                                                                          cGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAg                                                                     >  5:55182/1‑149 (MQ=255)
                                                                            gTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTaa                                                                                       >  8:47303/1‑129 (MQ=255)
                                                                            gTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTaa                                                                                       <  7:47303/129‑1 (MQ=255)
                                                                                cTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAg                                                               <  3:98638/149‑1 (MQ=255)
                                                                                cTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAg                                                               >  4:98638/1‑149 (MQ=255)
                                                                                 tCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAg                                                                                                    >  8:18849/1‑111 (MQ=255)
                                                                                 tCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAg                                                                                                    <  7:18849/111‑1 (MQ=255)
                                                                                   aCGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCAc                                                                                          >  5:155853/1‑119 (MQ=255)
                                                                                   aCGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCAc                                                                                          <  6:155853/119‑1 (MQ=255)
                                                                                         agaAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGaaa                                                       <  8:104108/148‑1 (MQ=255)
                                                                                          gaAACCCAACAATCCCAGCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGATATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAgg                                                     <  5:46196/149‑1 (MQ=255)
                                                                                                 aacaaTCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTgaag                                               >  4:72631/1‑148 (MQ=255)
                                                                                                    aaTCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACa                                           <  6:55182/149‑1 (MQ=255)
                                                                                                     aTCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACt                                                                                         <  1:104213/102‑1 (MQ=255)
                                                                                                     aTCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACt                                                                                         >  2:104213/1‑102 (MQ=255)
                                                                                                      tCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACAt                                          >  3:73939/1‑148 (MQ=255)
                                                                                                         cACCTGGTGTCAAGGAGAAGAACCTTAGCGGCAGTGACCTGCGGCATGCCCTCAGGATGATGTAACCCTCAGACAACGGGGATTTAGTCGAGATCACTAAGCCGATATGATc                                                                           >  4:106199/1‑112 (MQ=255)
                                                                                                         cACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAAc                                      >  2:61453/1‑149 (MQ=255)
                                                                                                                       gTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCct                        <  4:73939/149‑1 (MQ=255)
                                                                                                                         ggAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCctct                      <  3:32001/149‑1 (MQ=255)
                                                                                                                          gagaACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGt                                                                   >  3:124843/1‑103 (MQ=255)
                                                                                                                          gagaACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGt                                                                   <  4:124843/103‑1 (MQ=255)
                                                                                                                                cTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGtt               >  8:152004/1‑149 (MQ=255)
                                                                                                                                             cccTGCGGGCATGCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGTTGAAAACGTcacca  >  4:37668/2‑149 (MQ=255)
                                                                                                                                                tGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTa                                >  3:79632/1‑116 (MQ=255)
                                                                                                                                                tGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTa                                <  4:79632/116‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
ACCTATAGTGTGCCCGGAGTTCAGGGCGGGCATGGATGCTTAAATGAACCGCGAGTCTGTCTGGAATATTGAACCGGTAACTCACGATGAGAAACCCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACATGAACCCGTTAATCCGCCTCTGTGGGTTGAAAACGTCACCA  >  NC_000913/272184‑272473

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