Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,269,9940GT36.4% 42.7 / 13.1 56intergenic (+392/+222)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base G (12/23);  new base T (0/20);  total (12/44)
Fisher's exact test for biased strand distribution p-value = 2.17e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACC  >  NC_000913/3269845‑3270116
                                                                                                                                                     |                                                                                                                           
gaatcgcgcctGAACGGTTTATTACGTACATCATGTAAAACAGACCGTTAAGACGCTTTCTTTTGGGTATAGTGGCGTGGACAGTCATTCATCTTTCTGCCGCTCCTAAAGTAAAAACCCGCAGAATCGGGTTTTTACGTAAAACAGGTg                                                                                                                             <  2:883195/144‑1 (MQ=17)
    cccGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAAc                                                                                                                         <  1:640744/150‑1 (MQ=255)
         cgtaaCGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAcc                                                                                                                    <  1:1312794/147‑1 (MQ=255)
         cgGAATTGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAcc                                                                                                                    <  2:1793477/150‑1 (MQ=255)
         cgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAcc                                                                                                                    <  1:1016836/150‑1 (MQ=255)
         cgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAcc                                                                                                                    <  1:1972342/150‑1 (MQ=255)
          ggAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCATTCTGGCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGGTTTTAAGTAAAACAGGTGTAAcggcacg                                                                                                                   >  1:903871/1‑144 (MQ=14)
           gAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGa                                                                                                                  <  1:1002057/150‑1 (MQ=255)
           gAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGAATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGa                                                                                                                  <  1:1002291/150‑1 (MQ=255)
           gAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTATTGGGTATAGAGTCGTGGACAATCATTAATAATTATTCCCCTCCAAAACAAAAAACCCACCGAAGCGGGGTTTTACGGAAATCAGGTGAAACTGCCCgc                                                                                                                  >  2:1920248/1‑149 (MQ=14)
             tcGGTTTATTACGTACATCAGGTAAAACTTACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATa                                                                                                                <  2:1444904/149‑1 (MQ=255)
                gTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCACATACTTTTGGGTATAGAGTCTTGTAAAGTCAATCAAAATTATGACCCTCAAAAAGAAAAAAACAACAGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCAAttcgg                                                                                                             >  2:1105723/1‑146 (MQ=9)
                         gTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACAGACCGATAAGCCGctttctt                                                                                                    <  2:642891/150‑1 (MQ=255)
                          taaaTCAGGTAAAAGTGACAGATAAGCAGCATTCATTTTTGTATTGTGTCGTGGACAGTGTTTCTTATTACTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGctttcttt                                                                                                   <  2:972613/147‑1 (MQ=11)
                             atcatttAATACTTAACGATAAGCCGCTTTCTTTTGGGTATAGTGTCTTCGACAGTCATTCATCTTTCTGCCCCTCCAATAGTATAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                                                <  2:583871/144‑1 (MQ=11)
                              tgaGGTAATACTGAGCGATAAGCTGCTTGCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCAGCCCCTCCAAAAGTAAAAACCCGCCGTAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTggg                                                                                               <  1:540342/148‑1 (MQ=11)
                              tCATGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCTTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTggg                                                                                               <  1:596520/150‑1 (MQ=255)
                              acatgTAATACTGACCGATAAGCCGCTTTCTTTTTTCTATAGTGTCGTGGAAAGTCTTTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTggg                                                                                               <  1:495866/146‑1 (MQ=11)
                                atgTAAAACTGACCGATAAGCCGCTTTCTTTTGCGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGta                                                                                             <  1:1955890/148‑1 (MQ=37)
                                aGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCAACTTTCAGCCCCTCCAAAAGAAAAAACCCACCGCATCGGGGTTGTTAGTATAAAAGGTGTTACTGAACGTAATGGCGCTTTCTTTTGTGta                                                                                             >  2:1527784/1‑150 (MQ=2)
                                       aCTGACCGATAATCCGCATTCTTTTGGGGATAGTGTCGTCGAAAGTCATTCATCTTTCTGCCCCTCCAAAAGTACAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTc                                                                                      <  1:435900/150‑1 (MQ=11)
                                         tGACCGATAAGTTGCTTTCATTTGTTTATTTTGTCGTGGACAGTGTTTCATATTTCTGCCACTCCAAAAGTAAAAACCCGCCGAATCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGt                                                                                    <  1:588727/150‑1 (MQ=11)
                                         tGACCGATAAGCCGCTTTCTTTTGGGTATAGTCTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGt                                                                                    <  1:1659970/150‑1 (MQ=37)
                                             cGATAAGCCGCAGGCTTTAGGATTATGTGTCGTGGAAAGTCTTTCTTCTTTCTGCCCCTCAAAAAGTAAAAACCCGCCGTAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGAc                                                                                <  1:1238642/150‑1 (MQ=11)
                                               ttatgCCGCTTTCTTTTTGGTTTAGTGTCGTGGACAGTCATTCATCTTTATGCCACTCAAAAATTATAAACCCGCCGAATCGGGTTTTTACGTAAAACAGGTTAAACTGACCGGTAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAg                                                                              <  1:985809/146‑1 (MQ=11)
                                                  agacGCTTTCTTTTGTTTAATGTGTCGTCGAAAGTCATTCATCTTTATGCCCCTCCAAAATTATAAACCCGCCGAATCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtca                                                                           <  1:1610817/147‑1 (MQ=11)
                                                   gCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAAc                                                                                                                         >  1:949192/1‑103 (MQ=17)
                                                      gCTTTCTTTTGTGGATCGTGTCGTGGACTGGCATTCTTCTTTCTGCCCCTCCAAAAGGAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGtcattca                                                                       <  2:1646601/150‑1 (MQ=11)
                                                            ttttGTGTTTAGTGTCGTGGACAGTCATTAATCTATATGCCCCTCAACAAGTAAAAACCCGACGTATCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTc                                                                 <  2:1812784/150‑1 (MQ=11)
                                                               tGGTTATAGTGTCGTGGACTGTGTTTCATCTTTCTGGCACTCCAAAATTAAAAACCCGCCGTATCGGGTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGc                                                              <  2:1709416/150‑1 (MQ=11)
                                                                  tgctaGTGTCGTCGACAGTGTATCTTATATATGCCACTCCAAAAGAAAAAACCCGCCGTAGCGGGTTTTTACGTACAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGcccc                                                           <  1:1754892/147‑1 (MQ=11)
                                                                  ggtttgtgtCGTCGACAGTCAATCTTATTTCTGCCCCTCCAAAAGTATAAACCCGCCGCAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGcccc                                                           <  1:1447532/145‑1 (MQ=11)
                                                                            gTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCgg                                                                                                          >  2:1790962/1‑92 (MQ=17)
                                                                                            atatatGCCACTCAAAAATCACAAACCCGCCGAATCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCACTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCg                                 <  1:1278735/145‑1 (MQ=11)
                                                                                              ttatcgcccTCAAAAAGTACAAACCCGCCGAATCGGTTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCggg                              <  1:1328099/144‑1 (MQ=11)
                                                                                               tctcgcccTCCAAAATAACAAACCCGCCGTAGCGGGTTATTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGt                             <  1:1990256/145‑1 (MQ=11)
                                                                                                   gcacTCAAAAAGTATAAACCCGCCGTATCGGGTTTTTACGTACAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttttt                         <  2:940398/147‑1 (MQ=11)
                                                                                                      cTCCAAAATTAAAAACCCGCCGTAGCGGGTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCTTCGGCGGGtttttgct                      <  2:1924392/150‑4 (MQ=11)
                                                                                                      cTCCAAAAGTATATAACCGCCGATGCGGGTTTTTACGTAAATCAGGTGAAACTGTCGGATAATCCTCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACg                      <  1:1837330/150‑1 (MQ=11)
                                                                                                      cTCAAAAATTATAAACCCGCCGTTTCGGTTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACg                      <  2:6916/150‑1 (MQ=11)
                                                                                                       tcaaaaaGTAAAAACCCGCCGTAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCTTCGGCGGGtttttgctt                     <  2:1665109/147‑5 (MQ=11)
                                                                                                        ccAAAATCACAAACCCGCCGAACCGGTTTATTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTa                    <  1:602045/150‑1 (MQ=11)
                                                                                                        ccAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGc                                                                                                          >  2:1504483/1‑65 (MQ=255)
                                                                                                         cAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTAATCTTTCTTCCCCTCCAAAACCAAAAACCCGCCGGAACGGGAATTTACTTaa                   >  2:1912503/1‑150 (MQ=11)
                                                                                                             atctcaaaCCCGCCGTATCGGTTTTTTACGTACAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCa               <  2:1160716/145‑1 (MQ=11)
                                                                                                              gTAAAAACCCGCTTCGGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAg                                                   >  1:1361857/1‑114 (MQ=11)
                                                                                                              gTAAAAACCCGCTTCGGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAg                                                   <  2:1361857/114‑1 (MQ=11)
                                                                                                                acaaaCCCGCAGCAACGGTTTTTTACGTAAAACAGGTCAAACTGACCTATAAGCCGCTTTCATTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCCCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGt            <  1:122822/148‑1 (MQ=11)
                                                                                                                aaaaaCCCGCCGTATCGGGTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGt            <  2:1598867/150‑1 (MQ=255)
                                                                                                                 aaaaCCCGCCGTATCGGGTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTg           <  2:1821410/150‑1 (MQ=255)
                                                                                                                      ccgccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTAGAGTCTCGTGGAGAGTCATTCATCGTTCATCCCCAACAAAAACAAAAACACGCAGAAGGGCGGTTTTAAGGTAAACACGTGaaaaa      >  1:1074683/1‑148 (MQ=2)
                                                                                                                          cgtaGCGGTTTATTACGTATAACTGTTTAAACTGTGCGATATGCCGCTTTCTTTTGTTTTTTGTGTCGTGGACTGTCTTTCTTCTTTCTGCCCCTCCAAAAGCATAAACCCGCCGTAGCGGGTTTTTACGTAAATCAGGTGAAACTGAcc  <  2:1623032/147‑1 (MQ=2)
                                                                                                                          cgtaGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGAcc  <  1:1700741/147‑1 (MQ=255)
                                                                                                                          cgtaGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGAcc  <  1:1511851/147‑1 (MQ=255)
                                                                                                                          cggatcGGTTTTTTACGTATATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGAcc  <  1:192673/145‑1 (MQ=11)
                                                                                                                                gggTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGAGTCGTGGACAGtcattcat                                                                      >  2:870568/1‑77 (MQ=21)
                                                                                                                                         cGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccgcc                                       >  1:344203/1‑99 (MQ=37)
                                                                                                                                                     |                                                                                                                           
GAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACC  >  NC_000913/3269845‑3270116

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