Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,624,5170GT41.4% 75.5 / 26.6 58E492* (GAA→TAA) ppkpolyphosphate kinase, component of RNA degradosome
Reads supporting (aligned to +/- strand):  ref base G (7/27);  new base T (24/0);  total (31/27)
Fisher's exact test for biased strand distribution p-value = 2.48e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAAAACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTGAAGCTAAATAACCTTGTC  >  NC_000913/2624375‑2624666
                                                                                                                                                |                                                                                                                                                       
aGAAAACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCa                                                                                                                                                 <  1:55950/151‑1 (MQ=255)
aGAAAACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCa                                                                                                                                                 <  1:662706/151‑1 (MQ=255)
aGAAAACGGTAAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCa                                                                                                                                                 <  1:547517/151‑1 (MQ=255)
        gTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGTcc                                                                                                                                       <  1:339463/151‑1 (MQ=255)
        gTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGTcc                                                                                                                                       <  2:370316/151‑1 (MQ=255)
          gAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGTCCg                                                                                                                                      <  1:86193/150‑1 (MQ=255)
                  gCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGTCCGGTGACAttt                                                                                                                             <  1:342450/151‑1 (MQ=255)
                    gTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGa                                                                                                                           <  1:305043/151‑1 (MQ=255)
                      tATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTTAATTTTTTTAAAACCCATACCGccccggggcgttttttt                                                                                                                         >  2:377007/1‑137 (MQ=255)
                      tATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTTAATTTTTTTAAAAACCATAACGTCCCGTGGCAttt                                                                                                                             >  2:925342/1‑147 (MQ=255)
                       aTGCACACATCGGGGCGGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACTAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGatta                                                                                                                        <  2:898286/151‑1 (MQ=255)
                       aTGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCTCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGatt                                                                                                                         <  1:823548/150‑1 (MQ=255)
                         gCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTAACTTTTTTTAAAAACCATACCGTCCGGTGGCATTTGattatt                                                                                                                      >  2:972903/1‑151 (MQ=255)
                          cacaccTGGTGACCGGGACCTTTAACGAAAAATCCGCGCGTCTTTATACTGACTATTCGTTGCTGCCCGCCGATGCGCGCATCCCCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGattatt                                                                                                                      <  1:713750/150‑1 (MQ=255)
                                  gggACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTTAATTTATTTGAAAACCAAACCGTCCGGGGGCCTTTGATTTTTTTATTGTAt                                                                                                             >  1:835697/1‑151 (MQ=255)
                                  gggACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTAACTTTATTTACAACCCATACCGTCCGGGGACATTTGATTATTTTATTGt                                                                                                               >  1:172520/1‑149 (MQ=255)
                                   ggACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATc                                                                                                            <  1:940156/151‑1 (MQ=255)
                                   ggACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGTTTTTAATTTTTTTTAAACCCCTACCCTCCGGGGACATTTGTTTATTTTATgg                                                                                                                >  2:919331/1‑147 (MQ=255)
                                    gACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTTACTTTTTTTAAAAACCACACCGTCCGGGGGCATTTGATTATTTTATGGTATcg                                                                                                           >  2:107799/1‑151 (MQ=255)
                                      ccGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATcgcc                                                                                                         <  1:433594/151‑1 (MQ=255)
                                         ggAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCa                                                                                                      <  1:546023/151‑1 (MQ=255)
                                           aaCTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCaaa                                                                                                    <  2:654617/151‑1 (MQ=255)
                                           aaCTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTTACTTTTTTTAAAACCCATACCGTCCCGGGGCATTTGATTATTTTATTGTATCGCCGCaa                                                                                                     >  2:493947/1‑150 (MQ=255)
                                             cTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTTAATTTTTTTTAAAACCCTACCGGCCGGGGGGAATTTTTTTTTTAATTGTATTGCCCCcaaa                                                                                                   >  2:386182/1‑150 (MQ=255)
                                              tttAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTTAATTTTTTTAAAAACCCCACCCTCCGGGGGCATTTTATTTTTTTATGGGATCGCCGCcaaa                                                                                                   >  2:894849/1‑149 (MQ=255)
                                                tAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTc                                                                                                <  1:813317/150‑1 (MQ=255)
                                                 aaCGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTccc                                                                                              <  2:304142/151‑1 (MQ=255)
                                                   cGAAAAAACCGCGCGTCTTTAGTCTGACTATTCGTTGCTGACCGCCGATGCGCGCATCTCCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCccgc                                                                                            <  1:792911/151‑1 (MQ=255)
                                                     aaaaacCCCCGCGTCTTTATACTGACTATCCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCccgcc                                                                                           <  1:220282/150‑1 (MQ=255)
                                                            cgcgcgTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTAACTTTTTTTAAAACCCCTACCCTCCGGGGGCACTTTTTTTTTTTATTGTATCTCCCCAAAACcccccccccccctt                                                                                    >  1:322217/1‑146 (MQ=255)
                                                             gcgcgTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTTTTGAAAACCCCCACCCGCCGGGGGCATTTTTTTTTTTTATGGTATCGCCGCAAAACTCCCGCCGCCTCtt                                                                                    >  1:415017/1‑149 (MQ=255)
                                                              cgcgTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTTACTTTTTTTAAAAACCATAACCGCCCGGGGCATTTGATTATTTTATGGTATCCCCGCCAAAATCcccccccccttttt                                                                                  >  2:374856/1‑144 (MQ=255)
                                                               gcgTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTAt                                                                                <  1:95861/151‑1 (MQ=255)
                                                                 gTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTTATTCAATTTATTGAAAAACCATACCCCCCCGGCGAAATTGCTTTTTTTATGGTTTAACCCCCCAACTCCCGCCCCCTACTGttt                                                                                >  2:805673/1‑147 (MQ=255)
                                                                      tataCTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGTGGTTTTTATTTTTTGAAAAACCACACCCTCCCGGGACATTTTTTTTTTTTATTGGATTTCCCCAAAAAACCCCCCCCCTAtttttttgaa                                                                            >  1:563421/1‑140 (MQ=255)
                                                                          cTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGTTTTTACCTTTTTTTAAAACCCCACCCCCCCGGGGCAATTTGTTTTTTTATGGGTTCGCCCCAAAACTCCCCCCCCCTAttttttttaatgggggg                                                                     >  1:16318/1‑136 (MQ=255)
                                                                               tATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTTTTTAAAACCCAAACCGCCCGGGGGCATTTGATTATTTTATGGTATCGCCGCCAAAATCCCCCCCCCCATTGTTTGAAATGGTGGGcc                                                                   >  1:969369/1‑148 (MQ=255)
                                                                                     ttGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGGTTTTACTTTTTTTAAAACCCCTACCCCCCGGGGGCATTTTATTTTTTTATTGTATCGCCGCAAAACTCCCCCCGCCCCTTTTTTGAAATGGGGGGCCGCGAGGTc                                                           >  1:533587/1‑150 (MQ=255)
                                                                                       gCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTTAATTTTTTTAAAAACCAAACCCTCCGGGGGCATTTGTTTATTTTATGGTATCGCCCCCAAACTCCCCCCCCCCATTGTATGAAAAGGGGGACCGCGGGAACGcc                                                        >  2:624779/1‑151 (MQ=255)
                                                                                         tGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTAACTTTATTTAAAAACCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTTTATGAAATGGTGGACCGCGAGATCGcc                                                        >  2:897135/1‑149 (MQ=255)
                                                                                          gACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTTACTTTTTTTAAAACCCATACCGCCCGGTGGCATTTGATTTTTTTATTGTATTGCCGCCAAACTCCCCCCCCCTATTGTATGAAAATGTGGACCGCGGGGTCGCCCAc                                                     >  1:895130/1‑151 (MQ=255)
                                                                                          gACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCaa                                                      <  1:7312/150‑1 (MQ=255)
                                                                                           aCCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTTACTTTATTTAAAACCCAAACCGCCCGGTGGCATTTTATTTTTTATTGGAATCGCCGCCAAAATCCCGCCCCCCATTTTATAAAAATGTGGGACCCGAGACCGccca                                                      >  2:104738/1‑147 (MQ=255)
                                                                                             cgccgATGCGCGCATCACCAACGAAGTACGGCGGGGATTTTACTTTTTTTAAAAACCATACCCTCCGGGGGCATTTTATTTTTTTATGGGATCGCCGCAAAACTCCCCCCGCCCCTTGTTTGAAATGGGGGACCGCGGGATCGCCCAccc                                                   >  2:655918/1‑148 (MQ=255)
                                                                                                   tGCGCGCATCACCAACGAAGTACGGCGGGGTGTTAACTTTTTTTAAAAACCATACACTGCGGGGGGAATTTATTATTTAATTGGATCTGCCCAAAAAACCCGGCCCCCCTTTTATGAAAaggggggacccggggttg                                                          >  2:886663/1‑125 (MQ=255)
                                                                                                     cgcgcATCACCAACGAAGTACGGCGGGTATTTTAATTTTTTTAAAAACCATAACGTCCGGGGGCATTTGTTTTTTTTATAGGATCGCCGCAAAAATCCCGCCCCCTCTTGTATGAAATGGGGGACCGCGAGAACGCCAACGCGCAGCaa                                            >  2:401489/1‑149 (MQ=255)
                                                                                                        gcATCACCAACGAAGTACGGCGGGTGTTTTATTTTTTTTTAAACACAAACCCTCCCGGGGCATTTGGTTATTTAATTGTATTGCCGCAAAACTCCCCCCCCCCATTTTTTTAAAATGGGGAGCGCGGGGGCGCCAACCCCCCGCAAggg                                         >  1:624026/1‑149 (MQ=255)
                                                                                                         cATCACCAACGAAGTACGGCGGGTATTTTACTTTTTTTAAAAACCATAACGTCCCGGGGCATTTTATTTTTTTATTGTATCGCCCCAAAAAACCCCCCCCCCCTTTTATGAAAAGGTGGACCCCGAGGTCGCCCACCCCCCGCAAggg                                         >  2:603726/1‑148 (MQ=255)
                                                                                                              ccAACGAAGTACGGCGGGGATTTAACTTTATTGTAAAACCATACCGCCCGGGGGCATTTGATTATTTTATTGTATTGCCGCAAAAATCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCCGCAAGGGCTGccc                                   >  1:74359/1‑149 (MQ=255)
                                                                                                               cAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGt                                <  1:932659/151‑1 (MQ=255)
                                                                                                                aaCGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTg                               <  2:698068/151‑1 (MQ=255)
                                                                                                                   gAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTa                            <  2:649177/151‑1 (MQ=255)
                                                                                                                     aGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATc                          <  2:733779/151‑1 (MQ=255)
                                                                                                                              gggTATTTAACTTTATTGAAAACCCATACCGTCCGGTGACAGTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTGaa                  <  2:585908/151‑1 (MQ=255)
                                                                                                                                       aaCTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTGAAGCTAAAt           >  2:312962/1‑150 (MQ=255)
                                                                                                                                        aCTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTGAAGCTAAATa          >  2:71426/1‑150 (MQ=255)
                                                                                                                                          tttATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTGAAGCTAAATAAcc       >  2:180087/1‑151 (MQ=255)
                                                                                                                                               tGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTGAAGCTAAATAACCTTGt   <  2:449433/150‑1 (MQ=255)
                                                                                                                                                gAAACCCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTGAAGCTAAATAACCTTGTc  <  1:226399/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                       
AGAAAACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTGAAGCTAAATAACCTTGTC  >  NC_000913/2624375‑2624666

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