Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,151,5890AT48.1% 50.7 / 32.6 56L135F (TTA→TTTyegJuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (3/25);  new base T (26/0);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 1.96e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

GGCGTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCA  >  NC_000913/2151452‑2151736
                                                                                                                                            |                                                                                                                                                     
ggCGTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATc                                                                                                                                          <  1:413125/151‑1 (MQ=255)
   gTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATTCCATTTATTCTTTTTTTGTTACCTATcc                                                                                                                                         >  1:289733/1‑149 (MQ=255)
     aaTCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTTATCATTTATCCTTTTTTGGTCACCt                                                                                                                                             >  1:720067/1‑143 (MQ=255)
     aaTCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATAACAATTATCCttttttttttgccaattcaaa                                                                                                                                      >  2:466927/1‑135 (MQ=255)
        cAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATAACATTTATTCTTTTTTTGTTAccccacccaaaaa                                                                                                                                   >  1:453667/1‑139 (MQ=255)
          gCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTTAACCTTTTTTTGTTACCTATTCCAaagaaa                                                                                                                                 >  1:360633/1‑149 (MQ=255)
            aaTGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAaagaagaag                                                                                                                              <  1:56113/151‑1 (MQ=255)
               gACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTTAACCTTTTTTCTTTTTTTGTTACCTATCCCAaaaaaaaa                                                                                                                               >  1:306172/1‑147 (MQ=255)
                 cccTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCcgac                                                                                                                         <  1:57911/151‑1 (MQ=255)
                 cccTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGATATTACAATTATTCTTTTTTTGTTGCCTAAACaaaaaaaaaa                                                                                                                               >  2:356638/1‑145 (MQ=255)
                  ccTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGATATAACACTGCTCCCTTTTTTGTCACCTCTTCaaaaaaaaaaaaccac                                                                                                                         >  2:294285/1‑144 (MQ=255)
                     ggTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAg                                                                                                                      <  1:760451/150‑1 (MQ=255)
                     ggTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAg                                                                                                                      <  2:1043906/150‑1 (MQ=255)
                       tATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATAAAATTGATCCTTTTTTTGTGGCCTATACAAAAAAAAAAAcc                                                                                                                            >  2:759117/1‑142 (MQ=255)
                         ttGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGGTATAACATTTATACTTTTTTTGTTACCTATCCCAAAAAAAAAAACGACGAGCTa                                                                                                                   >  1:745050/1‑149 (MQ=255)
                           gTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAaga                                                                                                               <  1:909341/150‑1 (MQ=255)
                              aCAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCaa                                                                                                           <  1:670284/151‑1 (MQ=255)
                                 aaCGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAAt                                                                                                        <  2:713535/151‑1 (MQ=255)
                                  aCGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAAtt                                                                                                       <  1:852877/151‑1 (MQ=255)
                                     tAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATAACAATGTTCCTTTATTTGTTACCCATCCAAAAAAAGAAGACGACGAGCTAAAAGAAAAAATa                                                                                                      >  1:796458/1‑149 (MQ=255)
                                       gAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATAACCTTTATTCTTTTTTTGTGACCCATCCCAAAAAAAAAACCGACGAGGGAAGAGCaaaa                                                                                                         >  1:244395/1‑144 (MQ=255)
                                       gAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGGTTTAACATTTATCCTTTTTTTGGGGCCTATTCAAAAAAAAAAGGCGGACAGCTCAGGGGAAAAtattttt                                                                                                  >  2:16827/1‑145 (MQ=255)
                                        aaTATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGGTTTACAAATGGTTCTTTTTTGGTGGCCCTTGCAAAAAAAAAAGCCCAGGAGCTTAGAGGAAAAATaat                                                                                                    >  1:934984/1‑146 (MQ=255)
                                         atatGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTTATCATTTATCCTTTTTTTGTCACCTATCCCAAAGAAAAAACCGACGAGCTAAGAGCAAAAATAGtt                                                                                                   >  1:382660/1‑148 (MQ=255)
                                           atGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATATCATTTATCCTTTTTTTGTTACCTATTCAAAAAAAAAAAACCACGAACTAAAAACaaaa                                                                                                         >  1:977147/1‑140 (MQ=255)
                                             ggCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATAACACTTGTCCTTTTTTTTGCACCTATCCCAAAGAAGGAACAGGAGAGGTAAAAACAAAATTTGtt                                                                                                   >  2:477040/1‑144 (MQ=255)
                                              gCCAGGAGTGGAAAATCAAGAGAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATAACCATGATTCTTTTTTTGTTGCCCATTCAAAAAAAGAAGCCGACGGACCAAGAAGAAAAATTAttttccaa                                                                                              >  2:284964/1‑144 (MQ=255)
                                              gCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTTACCCTTTATAATTTTTTTTTCACCCATATCAAAAAAGAAAACGGACAGAGAAGAACAAAATTAGTTCGCTa                                                                                              >  2:458336/1‑148 (MQ=255)
                                               ccAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCa                                                                                          <  1:247085/151‑1 (MQ=255)
                                                cAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTTAACATTTATTCTTTTTTTGTTACCTATCCCAAAAAAAAAAACGACGAGGCAAAAGGAAAAATAATTCGCTaa                                                                                             >  1:559524/1‑147 (MQ=255)
                                                cAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCa                                                                                          <  2:559024/150‑1 (MQ=255)
                                                 aagggTGGAAAATCAAAGAAGGAAATTATTTCAGCTGGATTGTTTAGCCCCGGGGTTAAATTTACCGGCGGTAATCCCTTTAACCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGTCGAGCTAAGAGCAAAATTAGTTCGCTAATCAGc                                                                                        <  1:431839/147‑1 (MQ=255)
                                                  ggAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTTACCCTTGATCCTTTATTTGTCACCTCTCCAAAAAAAAAAACCGACCAGCTAAAAACAAAAATAATTTGCTAATCAAcc                                                                                       >  1:125789/1‑151 (MQ=255)
                                                    aGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGATATACCCTTGATCCTTTCTTTGGCGCCCAACCAAAAAAAAAAGACGGCCATCGAAGCGAAAAAATAGTTTGCtcatcaac                                                                                        >  2:349473/1‑146 (MQ=255)
                                                       ggAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTTATCATTGAGCCTTTTTTTGTCACCTATCCCAAAAAAAAAACCGCCCAGCTAAGAACAAAAATAGTTTGCTCATTAGcc                                                                                       >  1:761530/1‑146 (MQ=255)
                                                         aaaaTCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCGTGTATTGGTCACCTATCCAAAAGAAGAAACCGACGAGCGAAGAGCAAAATTAGTTCGCTAATCAGCCGTTatat                                                                                <  2:697770/151‑1 (MQ=255)
                                                          aaaTCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAATCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTAtatt                                                                               <  2:699394/151‑1 (MQ=255)
                                                           aaTCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTAtattt                                                                              <  2:348439/151‑1 (MQ=255)
                                                              cAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTAtatttat                                                                            <  1:696269/150‑1 (MQ=255)
                                                                  gaaagaaGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGtt                                                                        <  2:1052674/151‑1 (MQ=255)
                                                                   aaagaaGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTTTACAATTGATTCTTTTTTTGTTACCCATCCAAAAAAAGAAAACCACGAGCCAAAAGCAAAAATAGTTTGCTAATCAAACGCTATATTTTTTGttt                                                                       >  1:439935/1‑151 (MQ=255)
                                                                   aaagaaGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATATCAATTATTCTTTTTTTGTGACCTATACAAAAAAAAAAGACGACGAGCTAAGAGCAAAAATAATTCGCTAATTAACCGCTTTAtttttttttt                                                                       >  1:1011808/1‑151 (MQ=255)
                                                                    aagaagATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAATCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTa                                                                      <  2:786486/151‑1 (MQ=255)
                                                                    aagaagATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTa                                                                      <  1:605233/151‑1 (MQ=255)
                                                                     agaagaTATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTaa                                                                     <  2:706401/151‑1 (MQ=255)
                                                                          atatTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGCTTTATCCTTTATCCCTTTTTTGTTACCCATTCCAAAAAAAAAAACGGCGGGCTAAAAGCAAAAATAGTTTGCTTATAAACCCGTATATTTTTTGtttttaatg                                                                 >  1:580791/1‑144 (MQ=255)
                                                                                   aCTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATAACAATTTTTCCTTTTTTGTTCCCCCTTCCAAAAAAAAAAACCGCGAGCCAAAAGGAAAAATAGGTTGCTAATAAACCGTTATATTTTTTTTTTTAAATGGCgggt                                                           >  1:298284/1‑145 (MQ=255)
                                                                                      ggATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAAttg                                                    <  1:203319/150‑1 (MQ=255)
                                                                                          tGTATACGCGCGGGGATAAAATTTACGGCGGGTTTACCATTTATCCTTTTATGGTCACCTATACAAAAGAAGAAGACGGCGAGATAAGAGCAAAAATAATTTGCTCATCACCCGCTATATTTTTTGTTTAAAATGGCCGGTGGATttttgt                                                >  2:246113/1‑151 (MQ=255)
                                                                                                cgcgcgGGGATAAAATTTACGGCGGTTTTAACATTTATCCTTTTTTTGTCACCTATCCCAAAAAAGAAACCGACGAGGTAAAAGGAAAAATAATTTGCTAATCAGCCGTTTTATTTTTTTTTTTAAATTGGGGGTGTAttgttttttt                                             >  1:774825/1‑143 (MQ=255)
                                                                                                cgcgcgGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCaa                                          <  1:369052/151‑1 (MQ=255)
                                                                                                 gcgcgGGGATAAAATTTACGGCGGATATAACATTTATTCTTTTTTTGGCACCTATCCCAAAAAAAAAACCGACGAGCCAAAAGGCAAAATAGTTTGCTAATCAACCGTTATATTTTTTTTTTTAAATGGCGGGTgtaattttttttt                                              >  1:830250/1‑135 (MQ=255)
                                                                                                             aaTTGAGGGGGGTTATACATTTATTCTTTTGTTGGTCCCCTATCCAAAAGAAGAACCCGAGAGGCTAAGAGAAAATTTGGTTCCATAACCACCCTTTATATTTATTGTTAAATCGGGCGTTTAATTTGTTGTCTTCAACCGCCATGAGaaa                             <  2:325227/151‑1 (MQ=255)
                                                                                                                               ctatttacccTTTATTTGTCCCCTACCCAAAAAAAGAACCGGACAGGCAAAAAGAAAATTTAGTTCCTTAATCAGCCTTTATTTTTATTTTTAAATCGGGTGTTGAAATTTTTGTCTTCACCCGCCATGAGAAACTCACTCTGCCcaca              <  2:976509/141‑1 (MQ=255)
                                                                                                                                  ttGACCTTTTATTGGTCCCCTACCAAAAAAAAGAAGCCGACGAGTAAAGAGAAAAATTAGTTCGCTAATCAGCCTTTATATTTATTTTTTAAACTGGCGGTTGATTTTTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAgcg           <  1:653128/149‑1 (MQ=255)
                                                                                                                                      tCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAAACGCCATGAGAAACTCACTCTGCCCACAGCCCCCGTg     >  1:884269/1‑151 (MQ=255)
                                                                                                                                          ttATTGGTCACCTATCAAAAAGAAGAAGCCGACGAGCAAAGAGAAAAATTAGTTCGCTAATCAGCCTTTATATTTTTTGTTTAAACGGGCGGTTGATTTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCa  <  1:135031/150‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                     
GGCGTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCA  >  NC_000913/2151452‑2151736

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