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

TTGTGAAGATGGTGAGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTCCC  >  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)
        aTGGTGAGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGGTGGGGGGGTGGttgt                                                                                                                                                  >  1:535474/1‑151 (MQ=255)
          ggTGAGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGAttgtt                                                                                                                                                 <  2:671645/150‑1 (MQ=255)
          ggTGAGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGAttgtt                                                                                                                                                 <  1:2813/150‑1 (MQ=255)
          ggTGAGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGGTTGTCGTTAGTGGCTCGTGGTGGGGGGGGGgttgtt                                                                                                                                                 >  2:732100/1‑150 (MQ=255)
             gAGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTggg                                                                                                                                              <  1:430816/150‑1 (MQ=255)
               ggAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGc                                                                                                                                           <  1:331746/151‑1 (MQ=255)
                 aaGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCAt                                                                                                                                         <  2:110430/151‑1 (MQ=255)
                           ccTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTtgat                                                                                                                               <  1:103069/151‑1 (MQ=255)
                                  ccAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTg                                                                                                                      <  1:130624/151‑1 (MQ=255)
                                    aGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGGGCTGGGGGGGTGATTGTTGGGGGGATTTAAATTTATGACa                                                                                                                          >  1:172773/1‑146 (MQ=255)
                                     gCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGaaa                                                                                                                  <  2:71036/151‑1 (MQ=255)
                                         aTGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTg                                                                                                               <  1:18840/150‑1 (MQ=255)
                                              tctcTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCAtt                                                                                                         <  1:187539/150‑1 (MQ=255)
                                                 ctGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTg                                                                                                        <  1:100443/148‑1 (MQ=255)
                                                 ctGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTg                                                                                                        <  2:451806/148‑1 (MQ=255)
                                                     aaaaCAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGt                                                                                                 <  1:638710/151‑1 (MQ=255)
                                                      aaaCAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTg                                                                                                <  1:436565/151‑1 (MQ=255)
                                                       aaCAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCCGGGTGGGGGGGTGGTTGTTGGGGGGATTTAACTTTATGACAGGGAAATTTCCCtttggggggggg                                                                                               >  1:317814/1‑141 (MQ=255)
                                                               ggCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGGGCTGGGGGGGGGATTTTTGGGGGCATTAAAATTGGTGAGAGAGGAAATTGCATTTGGGGGGTGCCCTTAAcc                                                                                      >  1:3634/1‑151 (MQ=255)
                                                                    cAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCAt                                                                                 <  1:510979/151‑1 (MQ=255)
                                                                       ttCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGGTGGGGGGGTGGTTTTTTGGGGCATTTAAATTTATGACAATGGAAATGGCATTGGTGGGGGGGCGTTTAcc                                                                                      >  2:70146/1‑143 (MQ=255)
                                                                        tCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGGGGTGGAGGGGTGATTGTTTGGGGCATTTAACTTTATGAAAGGGGAAATGCCATTGGGGGGGGGGCGTTTAACtt                                                                                    >  1:687645/1‑144 (MQ=255)
                                                                         cGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTg                                                                           <  1:499542/151‑1 (MQ=255)
                                                                          gATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGGGGTGGGGGGGAGGTTTTTGGGGGGAATTAAATTGGTGACAAGGGAAATGGCATTTGGGGGGGGGCGGTTAACCTCATCGtt                                                                            >  1:579697/1‑149 (MQ=255)
                                                                            ttACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGGGGTGGGGGGGTGGTTTTTTGGGGCATTTAAATTGATGACAGGGGAAATGGCATTTGTGAGGGGGCGTTTACCCTCATCGTTTAt                                                                         >  1:674672/1‑150 (MQ=255)
                                                                                  aaaaGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGGTGGGGGGGGGATTTTTTGGGGGCGTTAAATTTATTAAAGTGAAAATTTCCTTTGGGGGGGGGCCGTTAACTTCATCGt                                                                             >  2:729073/1‑140 (MQ=255)
                                                                                        gAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGGGGTGGGGGGGGGGTTGTTTGGGGGATTTAAATTTATGGCAGGGAAAATGTCATTTGTGGGGGGGCGTTAAACCTTTTCGTTGATGGAGATCa                                                                 >  1:279367/1‑146 (MQ=255)
                                                                                        gAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGGGCTGGGGGGGTGGTTTTTGGGAGCATTAAAAATGATGACAGGGGAATTGTCATTTGTGAGGGGCCCGTTACCCTTATCGTTTATTCAGATCaa                                                                >  1:423087/1‑147 (MQ=255)
                                                                                          aGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGGGGTGGGGGGGTGATTGTTTGGGGGATTTAAATTGATGACAGGGGAATTTCCCTTGGGGGGGGGCCGTTTAACTTTATCCTTGATGGAGATAAATg                                                              >  1:230083/1‑147 (MQ=255)
                                                                                            tGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGGGGGGTGATTTTTGGGGGCATTTAAATTTGTGTCAGGGAAAATGCCATTTGGGGGGGGGCGGTTAACTTTATCGGTGATGCAGATCAAtaggtg                                                          >  1:745520/1‑144 (MQ=255)
                                                                                            tGAGATTTGGCAGGCAAAATCCATCGGCAGGGTTGTCGTTAGTGGCTCCTGGTGGGGGGGGGATTTTTGGGGGGATTTAAATTTTTTGCAGGGTAAAATTCACTGGGGAGGTGGCGTTTACCTTTTTCTTTGGTGTCGATCAATTGCTGt                                                         >  2:315122/1‑150 (MQ=255)
                                                                                             gagaTTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGtt                                                        <  1:200224/150‑1 (MQ=255)
                                                                                             gagaTTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTa                                                       <  2:325400/151‑1 (MQ=255)
                                                                                                aTTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCTTGGTGGGGGGGGGGTTGTTTGGGGCATTAAAATTTATGACAGGGAAAATTTCATTGGTGGGGGGCCGGTAAACTTCCTTGTTTGATCAGATCAAttggt                                                           >  2:716967/1‑140 (MQ=255)
                                                                                                  ttGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCCGGGTTGGGGGGTGATTGTTTGGGGGACTTAAATTGTAGAGAGGGGAAATTCCATTGGTGGGGGGGCGGTTAACTTCCTCGGTTGTGGAGATCAATTGCTGTTTTAggg                                                  >  1:488124/1‑151 (MQ=255)
                                                                                                       aggcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGGTGGGGGGGTGATTTTTGGGAGGATTTAAATTGATGACAGTGGAATTGCCATTGGGGAGGGGGCGTTTACCTTCCTCGTTGATGGAGATCAATAGGTGTTACAggggaagg                                             >  1:717138/1‑146 (MQ=255)
                                                                                                             aaTCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGcttctt                                       <  2:187046/151‑1 (MQ=255)
                                                                                                                      gCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGc                              <  2:183951/151‑1 (MQ=255)
                                                                                                                        aGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGa                            <  2:186487/151‑1 (MQ=255)
                                                                                                                                    aGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGtt                <  1:604359/151‑1 (MQ=255)
                                                                                                                                              gCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATtgtg       <  2:128944/150‑1 (MQ=255)
                                                                                                                                                  gAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTccc  <  1:235540/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                            
TTGTGAAGATGGTGAGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTCCC  >  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: