Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,097,5810AG42.6% 90.3 / 16.8 55noncoding (72/88 nt)serXtRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base A (4/27);  new base G (23/0);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 3.23e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGTCAACAAACTTTTT  >  NC_000913/1097444‑1097725
                                                                                                                                          |                                                                                                                                                     
tCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGGGGGGGGGGTTCGCAacccc                                                                                                                                            >  2:383067/1‑149 (MQ=255)
 cACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGa                                                                                                                                          <  2:160626/150‑1 (MQ=38)
  aCGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATa                                                                                                                                        <  2:324505/151‑1 (MQ=39)
   cGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGAGGGGGGGGtttaatcccctcgc                                                                                                                                         >  1:831507/1‑137 (MQ=1)
         gTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGGGGGGGGGGTTTGGACCCCCCATACCtt                                                                                                                                    >  1:803610/1‑148 (MQ=11)
                      tCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACtt                                                                                                                   <  2:438314/150‑1 (MQ=255)
                         ccAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCa                                                                                                               <  2:831692/151‑1 (MQ=255)
                               gTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGAAACCCCCATACGTTGCCGTATAAACACTTTCCAGGCg                                                                                                           >  2:573978/1‑149 (MQ=255)
                                   aaCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATCGAACACCCGACACGTTGTCGTATACACCCCTTCCCgggg                                                                                                           >  1:158321/1‑143 (MQ=14)
                                       ccGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCa                                                                                                <  2:325653/151‑1 (MQ=255)
                                         ggTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGGGGGGGGGGTTTGAACCCCCCCCACGGTTTCGGATAAACACTTTTCAGGCGGGGTCCCTCTGc                                                                                              >  1:326211/1‑151 (MQ=11)
                                              gggCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTc                                                                                         <  2:654119/151‑1 (MQ=255)
                                                   tctcATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGac                                                                                     <  2:866625/150‑1 (MQ=255)
                                                   tctcATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGac                                                                                     <  1:649778/150‑1 (MQ=255)
                                                    ctcATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGAGGGGGGGGGTCGCACCCCCCCTACGTTTCCGTTTAACCACCTTCCAGGGGGGCTCCCTctccccccc                                                                                         >  1:52803/1‑137 (MQ=14)
                                                       aTCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTTAGGGGGGGGTTTGAACCCCCCCTACGTTTCCGTATACACACTTTTCCAGGGGGCTCCTTCATCCACCCCGAGACCTc                                                                                >  2:964172/1‑151 (MQ=17)
                                                       aTCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGaca                                                                                    <  1:856937/147‑1 (MQ=255)
                                                         ccccccGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTCGGACCCCCCCTACGTTTCCGTATAAACACTTTCTCGGGGGGCTCCCTCCGCCACCCCGGCAcccc                                                                                >  1:652743/1‑147 (MQ=14)
                                                             ccGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCaaa                                                                         <  1:1024675/151‑1 (MQ=255)
                                                                 tgtgtgCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAAACCCCCCTACGTTGTCGGATACACACTTTCCAGGCGTGCTCCTTCCGCCACCCGGACACCCCCCCaaaattt                                                                     >  2:802015/1‑147 (MQ=14)
                                                                  gtgtgCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAACCCCCCATACGTTGCCGTATACACACTTTTCCGGCGTGCTCCTTCAACCCCCCGGACACCTCCCCaaaattt                                                                     >  2:656119/1‑146 (MQ=14)
                                                                   tgtgCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGAGGGGGGGGTTTGAAACCCCCATAACTTTCCGTATAAACCCCTTCCAGGCGGGGTCCTTCCACCACCCCGGCACCCCACCCAAtttt                                                                     >  1:713340/1‑147 (MQ=9)
                                                                        aaTATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAAACCCCCATACGTTGCCGTATACAAACCTTTCAGGCGGGCTCCCTCAGCCCCTCGGACACCCCACCAAAATTTCGTTc                                                                >  2:713503/1‑149 (MQ=14)
                                                                        aaTATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCt                                                              <  1:677214/151‑1 (MQ=255)
                                                                             aCGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTAtt                                                         <  1:666097/151‑1 (MQ=255)
                                                                             aCGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGGGTTGGACCCCCCATACGTTGCCGCGTAAAAACCTTTCCCGGGGGGTCCCTTCACCCCCCCGGCACCCCACCCAAATTTTGGTTCCttctt                                                          >  2:66847/1‑147 (MQ=2)
                                                                              cGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGCGAGGGGGGGGTTTGGACCCCCCCTACCTTTCCGGTTAACCAATTTCCCGGCGGGCCTCCTCCACCCCCCCGGCCCCCCCCCCAAATTTTCGt                                                                  >  1:655628/1‑142 (MQ=2)
                                                                                 aaaaaaaGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGAGGGGGGGAGTTGAACCCCCCCCAAGTTTCCGTATACACACATTCCCGGCGGGCCCCTTCCCCCACCCGGGCACCCCACCCAAATTTTGTTCCTGTCTTTCt                                                      >  1:378806/1‑150 (MQ=9)
                                                                                  aaaaaaGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTgg                                                   <  2:372460/151‑1 (MQ=255)
                                                                                     aaaGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAAc                                                <  2:58627/151‑1 (MQ=255)
                                                                                      aaGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACgg                                              <  2:6475/151‑1 (MQ=255)
                                                                                      aaGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACgg                                              <  1:665250/151‑1 (MQ=255)
                                                                                       aGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACgg                                              <  2:704826/151‑1 (MQ=255)
                                                                                         cccGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGGGGGGGGGGTTTGGAACCCCCGTACATTTTCGTATAAACACTTTTCAGGCGTGGGCCTTCCCCCCCCCCGGCACCTCCCCAAAATTTCGTTTCTGTCTTTCTGGGAACgg                                              >  1:922584/1‑150 (MQ=9)
                                                                                          ccGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGAACCCCCCACACGTTTGCGTATTCACACTTTTCAGGCGGGGTACTTCCGCCCCTCGGACACCTCCCCAAAATTTTGTTCCTGTCTTTCTGGGACGGgcg                                           >  1:230848/1‑151 (MQ=14)
                                                                                           cGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGGACCCCCCATACGTTTCCGTATTCACACCTTTCAGGCGGGCTCCCTCAGCCCCTCCGACACCTCCCCAAATTGTCGTTCCTGTCTTTCTGGGACggg                                             >  1:280609/1‑148 (MQ=14)
                                                                                                 ttCGTACGAGCTCTTCTTTAAATATGGCGGGGGGGGGGGGGGTTGAACCCCCCCTCAGTTTTCGTATTAACACTTTTCAGGCGGGGTCCCTCCGACCCCCCGGCACCCCCCCCAAATTTCGTTCCTTTTTTGTTGGAAAGGGGGCTCAttt                                    >  1:232677/1‑151 (MQ=2)
                                                                                                   cGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGCAACCCCCACACGTTGCCGTATAAACACCTTTCAGGCGGGGTCCTTTTACCCCTCCGGCACCCCCCCAAAATTTTGTTTCTCTCTTGCTGGGACGGGGGCCAAttt                                    >  2:960910/1‑149 (MQ=9)
                                                                                                      aCGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAggg                                <  2:204534/150‑1 (MQ=255)
                                                                                                        gAGCTCTTCTTTAAATATGGCGGTGGGGGGGGGGGTTGAACCCCCCCTAAGTTGCCCGATACACACTTTCCCGGCGGGCTCCTTCAGCCACCCGGACACCCCCCCCAAATTTCGGTCCTGTCTTGCTGGGAAGGGCGCCAATTTTgggga                              >  2:330603/1‑148 (MQ=9)
                                                                                                        gAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAACCCCCCATACGTTTCCCTATAAACACTTTTCCGGCGTGCTCCTTCAGCCACTCGGACACCCCACCCAAATTTTGTTTCTTTCTTTCTGGAAAGGGCGCCAATTTTggg                                >  2:16208/1‑148 (MQ=9)
                                                                                                            tcttctTTAAATATGGCGGTGAGGGGGGGGGTTGAAACCCCGACTTCTTTTCGTATAAAAACTTTCTCGGCGGGGGCCTTCAAGCACTCGGACACCTCCCCCAAATTTCGTTTCTGTCTTTCTGGGAAGGGCGCGCATTTTGGGGAAtta                          >  1:87094/1‑148 (MQ=9)
                                                                                                             cttcttTAAATATGGCGGTGAGGGGGGGGTTCGAACCCCCGATAAGTTGCCGTATACACACTTTCCAGGCGGGGTCCCTCAGCCACCCCGACACCCCACCAAAATGTCGTTTCTGTCTTTCTTGAACGGGGGCCAATTTTGGGgaa                             >  2:756214/1‑146 (MQ=14)
                                                                                                               tcttTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAGGGAAATCATg                        <  1:803870/149‑1 (MQ=255)
                                                                                                                    aaaTATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAGGGAAATCATGACCTGa                  <  2:185390/150‑1 (MQ=255)
                                                                                                                    aaaTATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAGGGAAATCACGACCTGa                  <  1:165659/150‑1 (MQ=255)
                                                                                                                      atatGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGt               <  1:388637/151‑1 (MQ=255)
                                                                                                                        atGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGTca             <  2:675668/151‑1 (MQ=255)
                                                                                                                          ggCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGTcaa            <  1:934175/150‑1 (MQ=255)
                                                                                                                            cGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGTcaaca          >  1:305933/1‑150 (MQ=255)
                                                                                                                                 gAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGTCAACAAACt      <  1:494281/150‑1 (MQ=255)
                                                                                                                                  aGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGTCAACAAACttt    <  2:427163/151‑1 (MQ=255)
                                                                                                                                  aGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGTCAACAAACttt    >  1:428222/1‑151 (MQ=255)
                                                                                                                                   gggggggATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGTCAACAAACtttt   <  2:491172/151‑1 (MQ=255)
                                                                                                                                    ggggggATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGTCAACAAACttttt  >  2:467180/1‑151 (MQ=255)
                                                                                                                                          |                                                                                                                                                     
TCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGTCAACAAACTTTTT  >  NC_000913/1097444‑1097725

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