Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,573,4960AG41.9% 53.8 / 22.3 43pseudogene (1083/1503 nt)yjiTpseudogene
Reads supporting (aligned to +/- strand):  ref base A (1/24);  new base G (18/0);  total (19/24)
Fisher's exact test for biased strand distribution p-value = 3.12e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGTGAAGATGGTGAGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGA‑TTGTTGGGAGCATTAAACTTGATGAC‑AG‑TG‑AAATTGCCA‑TTGGTGAGGTG‑CCGTTAACCTTCATCG‑TTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTCCC  >  NC_000913/4573349‑4573645
                                                                                                                                                   |                                                                                                                                                            
ttGTGAAGATGGTGAGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGGGCTgggggg                                                                                                                                                           >  1:128952/1‑151 (MQ=255)
   tgAAGATGGTGAGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGGGCTggggggg                                                                                                                                                          >  1:183957/1‑148 (MQ=255)
   tgAAGATGGTGAGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGGGCTGGGGGggtg                                                                                                                                                        >  1:333838/1‑151 (MQ=255)
        aTGGTGAGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGGTGGGGGGGTGG‑ttgt                                                                                                                                                  >  1:535474/1‑151 (MQ=255)
          ggTGAGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGA‑ttgtt                                                                                                                                                 <  2:671645/150‑1 (MQ=255)
          ggTGAGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGA‑ttgtt                                                                                                                                                 <  1:2813/150‑1 (MQ=255)
          ggTGAGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGGTTGTCGTTAGTGGCTCGTGGTGGGGGGGGGg‑ttgtt                                                                                                                                                 >  2:732100/1‑150 (MQ=255)
             gAGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGA‑TTGTTggg                                                                                                                                              <  1:430816/150‑1 (MQ=255)
               ggAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGA‑TTGTTGGGAGc                                                                                                                                           <  1:331746/151‑1 (MQ=255)
                 aaGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGA‑TTGTTGGGAGCAt                                                                                                                                         <  2:110430/151‑1 (MQ=255)
                           ccTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGA‑TTGTTGGGAGCATTAAACTtgat                                                                                                                               <  1:103069/151‑1 (MQ=255)
                                  ccAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGA‑TTGTTGGGAGCATTAAACTTGATGAC‑AG‑Tg                                                                                                                      <  1:130624/151‑1 (MQ=255)
                                    aGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGGGCTGGGGGGGTGA‑TTGTTGGGGGGATTTAAATTTATGAC‑a                                                                                                                          >  1:172773/1‑146 (MQ=255)
                                     gCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGA‑TTGTTGGGAGCATTAAACTTGATGAC‑AG‑TG‑aaa                                                                                                                  <  2:71036/151‑1 (MQ=255)
                                         aTGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGA‑TTGTTGGGAGCATTAAACTTGATGAC‑AG‑TG‑AAATTg                                                                                                               <  1:18840/150‑1 (MQ=255)
                                              tctcTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGA‑TTGTTGGGAGCATTAAACTTGATGAC‑AG‑TG‑AAATTGCCA‑tt                                                                                                         <  1:187539/150‑1 (MQ=255)
                                                 ctGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGA‑TTGTTGGGAGCATTAAACTTGATGAC‑AG‑TG‑AAATTGCCA‑TTg                                                                                                        <  1:100443/148‑1 (MQ=255)
                                                 ctGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGA‑TTGTTGGGAGCATTAAACTTGATGAC‑AG‑TG‑AAATTGCCA‑TTg                                                                                                        <  2:451806/148‑1 (MQ=255)
                                                     aaaaCAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGA‑TTGTTGGGAGCATTAAACTTGATGAC‑AG‑TG‑AAATTGCCA‑TTGGTGAGGt                                                                                                 <  1:638710/151‑1 (MQ=255)
                                                      aaaCAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGA‑TTGTTGGGAGCATTAAACTTGATGAC‑AG‑TG‑AAATTGCCA‑TTGGTGAGGTg                                                                                                <  1:436565/151‑1 (MQ=255)
                                                       aaCAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCCGGGTGGGGGGGTGG‑TTGTTGGGGGGATTTAACTTTATGAC‑AG‑GG‑AAATTTCCC‑tttggggggggg                                                                                               >  1:317814/1‑141 (MQ=255)
                                                               ggCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGGGCTGGGGGGGGGA‑TTTTTGGGGGCATTAAAATTGGTGAG‑AGAGG‑AAATTG‑CA‑TTTGGGGGGTG‑CCCTTAAcc                                                                                      >  1:3634/1‑151 (MQ=255)
                                                                    cAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGA‑TTGTTGGGAGCATTAAACTTGATGAC‑AG‑TG‑AAATTGCCA‑TTGGTGAGGTG‑CCGTTAACCTTCAt                                                                                 <  1:510979/151‑1 (MQ=255)
                                                                       ttCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGGTGGGGGGGTGG‑TTTTTTGGGGCATTTAAATTTATGAC‑AA‑TG‑GAAATGGCA‑TTGGTGGGGGG‑GCGTTTAcc                                                                                      >  2:70146/1‑143 (MQ=255)
                                                                        tCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGGGGTGGAGGGGTGA‑TTGTTTGGGGCATTTAACTTTATGAA‑AG‑GG‑GAAATGCCA‑TTGGGGGGGGG‑GCGTTTAACtt                                                                                    >  1:687645/1‑144 (MQ=255)
                                                                         cGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGA‑TTGTTGGGAGCATTAAACTTGATGAC‑AG‑TG‑AAATTGCCA‑TTGGTGAGGTG‑CCGTTAACCTTCATCG‑TTg                                                                           <  1:499542/151‑1 (MQ=255)
                                                                          gATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGGGGTGGGGGGGAGG‑TTTTTGGGGGGAATTAAATTGGTGACAAG‑GG‑AAATGGCAT‑TTGGGGGGGGG‑CGGTTAACC‑TCATCG‑tt                                                                            >  1:579697/1‑149 (MQ=255)
                                                                            ttACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGGGGTGGGGGGGTGG‑TTTTTTGGGGCATTTAAATTGATGAC‑AG‑GG‑GAAATGGCA‑TTTGTGAGGGG‑GCGTTTACCCTCATCG‑TTTAt                                                                         >  1:674672/1‑150 (MQ=255)
                                                                                  aaaaGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGGTGGGGGGGGGATTTTTTGGGGGCGTTAAA‑TTTATTAA‑AG‑TGAAAATTTCCT‑TTGGGGGGGGG‑CCGTTAA‑CTTCATCG‑t                                                                             >  2:729073/1‑140 (MQ=255)
                                                                                        gAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGGGGTGGGGGGGGGG‑TTGTTTGGGGGATTTAAATTTATGGC‑AG‑GG‑AAAATGTCA‑TTTGTGGGGGG‑GCGTTAAACCTTTTCG‑TTGATGGAGATCa                                                                 >  1:279367/1‑146 (MQ=255)
                                                                                        gAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGGGCTGGGGGGGTGG‑TTTTTGGGAGCATTAAAAATGATGAC‑AG‑GG‑GAATTGTCA‑TTTGTGAGGGGCCCGTTACCCTT‑ATCG‑TTTATTCAGATCaa                                                                >  1:423087/1‑147 (MQ=255)
                                                                                          aGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGGGGTGGGGGGGTGA‑TTGTTTGGGGGATTTAAATTGATGAC‑AG‑GG‑GAATTTCCC‑TTGGGGGGGGG‑CCGTTTAACTTTATCC‑TTGATGGAGATAAATg                                                              >  1:230083/1‑147 (MQ=255)
                                                                                            tGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGGGGGGTGA‑TTTTTGGGGGCATTTAAATTTGTGTC‑AG‑GG‑AAAATGCCATTTGGGGGGGGG‑CGGTTAA‑CTTTATCG‑GTGATGCAGATCAAtaggtg                                                          >  1:745520/1‑144 (MQ=255)
                                                                                            tGAGATTTGGCAGGCAAAATCCATCGGCAGGGTTGTCGTTAGTGGCTCCTGGTGGGGGGGGGA‑TTTTTGGGGGGATTTAAATTTTTTGC‑AG‑GG‑TAAAATTCA‑CTGGGGAGGTG‑GCGTTTACCTTTTTCT‑TTGGTGTCGATCAATTGCTGt                                                         >  2:315122/1‑150 (MQ=255)
                                                                                             gagaTTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGA‑TTGTTGGGAGCATTAAACTTGATGAC‑AG‑TG‑AAATTGCCA‑TTGGTGAGGTG‑CCGTTAACCTTCATCG‑TTGATGCAGATCAATGGCTGtt                                                        <  1:200224/150‑1 (MQ=255)
                                                                                             gagaTTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGA‑TTGTTGGGAGCATTAAACTTGATGAC‑AG‑TG‑AAATTGCCA‑TTGGTGAGGTG‑CCGTTAACCTTCATCG‑TTGATGCAGATCAATGGCTGTTa                                                       <  2:325400/151‑1 (MQ=255)
                                                                                                aTTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCTTGGTGGGGGGGGGG‑TTGTTTGGGGCATTAAAATTTATGAC‑AG‑GG‑AAAATTTCA‑TTGGTGGGGGG‑CCGGTAAACTTCCTTGTTTGAT‑CAGATCAAttggt                                                           >  2:716967/1‑140 (MQ=255)
                                                                                                  ttGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCCGGGTTGGGGGGTGA‑TTGTTTGGGGGACTTAAATTGTAGAG‑AG‑GG‑GAAATTCCA‑TTGGTGGGGGG‑GCGGTTAACTTCCTCG‑GTTGTGGAGATCAATTGCTGTTTTAggg                                                  >  1:488124/1‑151 (MQ=255)
                                                                                                       aggcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGGTGGGGGGGTGA‑TTTTTGGGAGGATTTAAATTGATGAC‑AG‑TG‑GAATTGCCA‑TTGGGGAGGGG‑GCGTTTACCTTCCTCG‑TTGATGGAGATCAATAGGTGTTACAggggaagg                                             >  1:717138/1‑146 (MQ=255)
                                                                                                             aaTCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGA‑TTGTTGGGAGCATTAAACTTGATGAC‑AG‑TG‑AAATTGCCA‑TTGGTGAGGTG‑CCGTTAACCTTCATCG‑TTGATGCAGATCAATGGCTGTTACAGGGACAGGcttctt                                       <  2:187046/151‑1 (MQ=255)
                                                                                                                      gCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGA‑TTGTTGGGAGCATTAAACTTGATGAC‑AG‑TG‑AAATTGCCA‑TTGGTGAGGTG‑CCGTTAACCTTCATCG‑TTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGc                              <  2:183951/151‑1 (MQ=255)
                                                                                                                        aGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGA‑TTGTTGGGAGCATTAAACTTGATGAC‑AG‑TG‑AAATTGCCA‑TTGGTGAGGTG‑CCGTTAACCTTCATCG‑TTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGa                            <  2:186487/151‑1 (MQ=255)
                                                                                                                                    aGTGGCTCGTGCTGGAGGGATGA‑TTGTTGGGAGCATTAAACTTGATGAC‑AG‑TG‑AAATTGCCA‑TTGGTGAGGTG‑CCGTTAACCTTCATCG‑TTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGtt                <  1:604359/151‑1 (MQ=255)
                                                                                                                                              gCTGGAGGGATGA‑TTGTTGGGAGCATTAAACTTGATGAC‑AG‑TG‑AAATTGCCA‑TTGGTGAGGTG‑CCGTTAACCTTCATCG‑TTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATtgtg       <  2:128944/150‑1 (MQ=255)
                                                                                                                                                  gAGGGATGA‑TTGTTGGGAGCATTAAACTTGATGAC‑AG‑TG‑AAATTGCCA‑TTGGTGAGGTG‑CCGTTAACCTTCATCG‑TTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTccc  <  1:235540/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                            
TTGTGAAGATGGTGAGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGA‑TTGTTGGGAGCATTAAACTTGATGAC‑AG‑TG‑AAATTGCCA‑TTGGTGAGGTG‑CCGTTAACCTTCATCG‑TTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTCCC  >  NC_000913/4573349‑4573645

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