Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,423,2970GC43.4% 83.3 / 29.5 53intergenic (+103/+126)yhdZ/rrfFputative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon
Reads supporting (aligned to +/- strand):  ref base G (10/20);  new base C (23/0);  total (33/20)
Fisher's exact test for biased strand distribution p-value = 2.22e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCG  >  NC_000913/3423150‑3423446
                                                                                                                                                    |                                                                                                                                                      
aaaTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCggg                                                                                                                                                      <  1:686315/150‑1 (MQ=255)
    cAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCCCCCGGCCg                                                                                                                                                   >  1:714769/1‑149 (MQ=255)
      gagCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGCTTCTCCCCCCCCCCCTCCCGGcc                                                                                                                                                    >  2:108976/1‑146 (MQ=255)
           tACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCACCCTCCCGGCCGGTGCtt                                                                                                                                             >  1:388594/1‑148 (MQ=255)
            aCGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCa                                                                                                                                         <  1:275120/151‑1 (MQ=255)
               aGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCCGCCGCTGCtt                                                                                                                                             >  2:286358/1‑144 (MQ=255)
                 ggCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTGCTCCCCCCCCCCCCCCCGGGCGGTGGTTctccacagt                                                                                                                                    >  2:53223/1‑143 (MQ=25)
                      ttttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCGCCCCCCTCCCGGCCGGTGCTTCGCCAccattaaa                                                                                                                                >  2:688465/1‑143 (MQ=255)
                       tttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAt                                                                                                                              <  1:9080/150‑1 (MQ=255)
                          aTCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCCCCCGGCCGTTGCTTCGCAACGTTTAAAAcccc                                                                                                                          >  1:366678/1‑149 (MQ=255)
                          aTCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCGGGCCGTGCTTTTTAACGTTCAAAtac                                                                                                                            >  2:177662/1‑147 (MQ=255)
                               aGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcg                                                                                                                     <  1:290970/151‑1 (MQ=255)
                                ggTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGTCCGCCCCCCCCCCCCCCGGCCCCTGCTTCCCCACCTTCAAAAtccctcccc                                                                                                                     >  2:260383/1‑149 (MQ=255)
                                 gTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCACCCCCCCGGCCGTTGCTTCGCAAAGTTCAAATCCCCAccc                                                                                                                      >  1:215759/1‑148 (MQ=21)
                                     tCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCTCCCCCCCCCCCCCCGGCCGGTGCTTTGCAACGGTCAAATCCGCcccc                                                                                                                      >  2:895254/1‑144 (MQ=21)
                                       cATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTCCCGCCCCCCCCCCCCCCGGCCCGTGCTTTGCAACAGTCAAAATCGCCCCccgccggtt                                                                                                              >  2:811707/1‑145 (MQ=255)
                                               gAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTAc                                                                                                      <  1:785590/150‑1 (MQ=255)
                                                aaTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTc                                                                                                    <  2:463048/151‑1 (MQ=255)
                                                 aTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCg                                                                                                   <  1:552472/151‑1 (MQ=255)
                                                   gTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCggg                                                                                                 <  1:385786/151‑1 (MQ=255)
                                                      aGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCGCCCCCCCCCCGGCCGTTGCTTTGCAACGCTTCAAACCGCCCCCCGCGGGTTTGTTCTACTcccggg                                                                                               >  1:690313/1‑145 (MQ=17)
                                                       gTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCGGCCGGTGCTTCGCAACGTTTAAATTCCCCCCCCGCGGGGTTTTCCTCCTCcgggg                                                                                               >  2:429747/1‑147 (MQ=255)
                                                         tCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAgt                                                                                            <  2:838772/150‑1 (MQ=255)
                                                           gAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTCCCGCCCCCCACCCCCCGGGCCGTGTCGTCGCCACGTTCAAAATCGCTCCCCGCGGGTTTTTTCCACTccgcggaggatt                                                                                         >  2:139822/1‑140 (MQ=255)
                                                             aaaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCCCCCCCCCCCCCCCCGGCCGCTGCTTCTCCAAGGTTAAATCCCCCCCCCGCGGAGTTGTTCCACTCCGGGGAgagt                                                                                          >  1:657964/1‑148 (MQ=12)
                                                              aaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTGCCGCCCCCCACCCACCCGGCCGGTGCTTCGCAACGGTCCAATCCGCTCCCCGCGGGGTTGTGCTACTCCGGAGAGGGttt                                                                                        >  1:548654/1‑148 (MQ=17)
                                                                  aaaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGac                                                                                  <  2:694080/151‑1 (MQ=255)
                                                                  aaaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCCCCCCCCCGGCCGGTGCTTCGCAACCGTCACATCCGCCCCCCGCGGGGTTGTTCTTCTCCGGGGAGGGTTTACCgcc                                                                                  >  2:52728/1‑149 (MQ=9)
                                                                      aGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaac                                                                              <  2:150841/151‑1 (MQ=255)
                                                                       ggCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaaca                                                                             <  1:821266/151‑1 (MQ=255)
                                                                        gCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaacaa                                                                            <  2:370690/151‑1 (MQ=255)
                                                                        gCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCCCCCCCCCGGCCCGTGCTTCGCAACGTTTAAAACCCCCCCCCGCGGGGTTTTTCTACTCCGGAGAGGGGTTACCGacaaaca                                                                             >  1:184802/1‑150 (MQ=9)
                                                                          cATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTTCCCCCCCCCCCCCCGGCCGTTTCTTCGCAACGTTCAAATCCCCCCCCCGCGGGTTTTTCCTTCTCCGGAGAGGGTTTACCCaccaaaaa                                                                            >  1:475418/1‑149 (MQ=14)
                                                                           aTCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCACCCCCCCGGCCGGTTCTTCGCAACGTTCAAATCCGCCCCCCGCGGGTTTTTCCTACTCCGGAGAGGGGTTACCCacaaaaaa                                                                            >  1:589438/1‑148 (MQ=14)
                                                                               tttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCCCCGGGCCGTTGCTTTCCAACGTTCAAATCCGCCCCCCGCGGGTTTTTCCTACCCGGGAGAGTGTTTACCCACAAaccaca                                                                          >  1:445181/1‑146 (MQ=14)
                                                                                   ggAGGGCCTCTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAAc                                                                  <  1:571585/150‑1 (MQ=255)
                                                                                                cTTGATTTGATGTCTGGCAGTTTATGGCGGGCGCCCCGCCCCCCCCCCCCCCGGCCGCTGTTTCGCAAAGTTTAAAACCGCCCCCCGCGGGTTTTTCCCTATCGGGGGGGAGGTCTCCGCCAAACAACAAAAAAAACAAAAAGCCAAggtt                                                    >  1:508793/1‑148 (MQ=255)
                                                                                                  tgatttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCt                                                    <  2:114086/149‑1 (MQ=40)
                                                                                                    atttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCg                                                <  2:205311/151‑1 (MQ=38)
                                                                                                        gatGTCTGGCAGTTCATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACt                                             <  1:829481/150‑1 (MQ=38)
                                                                                                         atGCCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCGGCCGGTGCTTTTCAAACTTTAAATCCCCCCCCCGCGGGGTTTTTCTTCTCCGGGGAGGGGTTACCCACAAACAAAAAATAAAAAAAAAAGGCCCGGCTTcccgc                                              >  2:338722/1‑145 (MQ=0)
                                                                                                            tCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGTTTCTTTGCAACAGTCAAAACCGCCCCCCGCGGGTTTTTTCTACTCCGGGGAGGGGTCTCCCACAAAAAACAGAAAAAAAAAAAAGCCCAGGCTTTCGGCTGGGcc                                        >  1:131438/1‑151 (MQ=2)
                                                                                                                  aGTTTATGGCGGGCGTCCTGCCCGCCCCCCCCCCGGCCGGTGCTTCGCAACGTTCAAATCCGCTCCCCGCGGGTTTGTCCTCCTCCGGGGAGGGGTCACCGACAAAAAAAAAATAAAAAAAAAAGGCCAGGCTTCCCACTTGGCCcttttt                                  >  2:999981/1‑149 (MQ=2)
                                                                                                                      tATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGtttt                               <  2:43270/150‑1 (MQ=38)
                                                                                                                         ggcgggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAttt                           >  2:429807/1‑151 (MQ=38)
                                                                                                                             ggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATg                       >  1:335800/1‑151 (MQ=38)
                                                                                                                                gTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCt                    >  2:679118/1‑151 (MQ=40)
                                                                                                                                     gCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGc                 <  1:453485/149‑1 (MQ=40)
                                                                                                                                     gCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGc                 >  1:425888/1‑149 (MQ=40)
                                                                                                                                     gCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCa                <  2:98371/150‑1 (MQ=40)
                                                                                                                                        ccGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGtt             <  2:421880/151‑1 (MQ=40)
                                                                                                                                            cACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCt         >  2:102610/1‑151 (MQ=40)
                                                                                                                                            cACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCt         >  2:730549/1‑151 (MQ=40)
                                                                                                                                                   cGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCg  >  2:668911/1‑151 (MQ=40)
                                                                                                                                                    |                                                                                                                                                      
AAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCG  >  NC_000913/3423150‑3423446

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: