Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,222,2330GA44.9% 51.1 / 27.2 49pseudogene (60/75 nt)ycgHpseudogene; putative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand):  ref base G (22/5);  new base A (0/22);  total (22/27)
Fisher's exact test for biased strand distribution p-value = 1.64e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTG  >  NC_000913/1222087‑1222382
                                                                                                                                                   |                                                                                                                                                      
gAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGaaa                                                                                                                                                    <  2:1021066/150‑1 (MQ=255)
gAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAat                                                                                                                                                   >  1:639131/1‑151 (MQ=255)
   aGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAatat                                                                                                                                                 >  2:325638/1‑150 (MQ=255)
          cAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTg                                                                                                                                         <  1:332079/151‑1 (MQ=255)
          cAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTg                                                                                                                                         <  1:728917/151‑1 (MQ=255)
            aGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAt                                                                                                                                       >  2:717673/1‑151 (MQ=255)
             gATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATc                                                                                                                                      >  2:780897/1‑151 (MQ=255)
                      aGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGtt                                                                                                                              >  1:528719/1‑150 (MQ=255)
                          aTTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCt                                                                                                                          >  2:622553/1‑150 (MQ=255)
                                   aaaTTAAACTGGGCGAAAAGGGCACCCTGAACCCGGCGGCCAGCGTCGGGGGAAGGGGGGCGTGCGGGGGGGTGATAATGGCTAAAATGACCCCGCAGTGAGGGTGGGCCTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAAcc                                                                                                                  <  1:346306/149‑1 (MQ=255)
                                        aatctgggggaaaaTGGCACCCTGAACCCGGGGGCCGGCGGGGGGGGAATGTGGGCGTGCGGCGGGTTAATAATGGCTACAATGACCCCGAGTTAATGGGGGGCCTAAAATATAGGTTCTGATCCCGCCGTTAGCTAAAAAACCGcgtcgt                                                                                                           <  2:273578/140‑1 (MQ=255)
                                               gCGAAAAGGGCACCCTGACGCCAGCGCCCAGCGTCGGGGGAAGGTGGGCGTGCAGGTGGGTGATAAGGGCTACAATCACCCCGAGGTGATGGGGGGCTTAAAATAAAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTAtt                                                                                                        <  2:43670/147‑1 (MQ=255)
                                                 gTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGACCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCg                                                                                                   >  2:424630/1‑150 (MQ=255)
                                                   aaaTGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGACGTTAGCTAAAAAACCGCGTCGTATTCATCGa                                                                                                  >  2:266/1‑149 (MQ=255)
                                                   aaaGGGCACCCTGAACCCAGGGGCCAGCGTGGGGGGAGTGGGGGCGTGCGGCGGGGTAATAATGGCTACAGTGCCCCCGAAGTGATGGGGGGCGTAAAATATATTTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGa                                                                                                  <  2:135314/149‑1 (MQ=255)
                                                       ggaacccTGAACCCAGCGGCCAGCGTGTGGGGATTGTGGGCGTGCAGTTGGGTGATAATGGCTACAATGACCCCGCAGTGATGGGGGGCCTAAAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCgg                                                                                             <  1:263169/145‑1 (MQ=255)
                                                          aaCCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTAc                                                                                          >  1:918353/1‑150 (MQ=255)
                                                          aaCCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTAc                                                                                          >  1:595634/1‑150 (MQ=255)
                                                                  tCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAat                                                                                 >  2:993635/1‑151 (MQ=255)
                                                                   ccAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAat                                                                                 >  2:847918/1‑150 (MQ=255)
                                                                        ggCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGTCTAAAATGACCCCGCAGTGAGGGTGGGCATAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtattat                                                                            <  1:406866/150‑1 (MQ=255)
                                                                        gcccAGCGTCGGGGGAATGTGGGGGTGCAGCTGGGTGATAAGGGCTACAATGCCCCCGCGGTGAGGGTGGGCCTAAAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATattatt                                                                           <  1:215386/149‑1 (MQ=255)
                                                                           cAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAttttt                                                                        >  1:527805/1‑151 (MQ=255)
                                                                            aGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAtttttt                                                                       >  2:382647/1‑151 (MQ=255)
                                                                             gcttCTGGGGATTGTGGGCGTCCGGGTGGGTGATAATGGCTACAAACCCCCCGCAGTGATGGTGGCCTGAAAATAAAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAtttttt                                                                       <  1:196452/147‑1 (MQ=255)
                                                                                   ggggAGGGTGGGGGGGCGGCGGGGTAATAATGGCTAAAAGGACACCGAAGGGGTGGTGGGCAAAAAAAATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATa                                                                <  2:639097/151‑1 (MQ=255)
                                                                                          gtgGGCGTGCAGGTGGGTGATAATGGCTAAAATGCCCCCGCAGTGAGGGGGGGCTTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACg                                                           <  1:982633/149‑1 (MQ=255)
                                                                                            gggCGTGCAGCGGGGTGAAAGGGCTTAAAATGACCCCGAAGTGGTGGTGGGCCTGAATTATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCAt                                                        <  1:984511/150‑1 (MQ=255)
                                                                                                gTGCAGCTGGTAGATAATGGCAAAAAGGACCCCGCAGGGATGGGGGGCTTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGgtg                                                    <  1:865484/150‑1 (MQ=255)
                                                                                                   cAGCGGGGTAATAATGCCTACAATGCCCCCGAAGTGAGGGTGGGCCTAAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGgt                                                     <  1:669935/146‑1 (MQ=255)
                                                                                                       tGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAAt                                            >  1:416991/1‑151 (MQ=255)
                                                                                                             aTAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGt                                      >  2:663337/1‑151 (MQ=255)
                                                                                                             aTAATGACAACAAGGACACGGCAGTGATGGGGGGCAAAAAATAAAGGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGc                                         <  2:275575/148‑1 (MQ=255)
                                                                                                               aagggCTAAAGTGCCACCGCAGTGAGGGGGGGCTTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGcc                                        <  1:630438/144‑1 (MQ=255)
                                                                                                                attgaataaaaTGACCCCGCAGGGAGGGTGGGTCTAAAAAATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAAt                                            <  2:335213/134‑1 (MQ=255)
                                                                                                                  gggttaaaaTGACCCCGCAGTGAGGGGGGGCCTAAAATATAATTTCTGCTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGAt                                   <  1:666078/143‑1 (MQ=255)
                                                                                                                  gggttaaaaTGACCCCGAAGAGATGGGGGGCCTAAAATATAGTTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGAt                                   <  2:923337/143‑1 (MQ=255)
                                                                                                                         aacaaccccGGAGTGATGGTGGGACTAAAATATAATTTGTGATTCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAg                             <  2:821420/142‑1 (MQ=255)
                                                                                                                          aTGACCCCGCAGTGATGGGGGGCGTAAAATATAATTTTTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCAt                          <  1:569327/151‑1 (MQ=255)
                                                                                                                             accccgaagttttggtggGCTTAAAATATATTTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCAttt                        <  1:780933/138‑1 (MQ=255)
                                                                                                                              cacCGCAGTGTTGGGGGCCCAAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCAttt                        <  1:294761/149‑1 (MQ=255)
                                                                                                                              cacCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGc                      >  2:364525/1‑151 (MQ=255)
                                                                                                                                 cGCAGTGATGGTGGGCCTAAAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCg                     <  2:947636/149‑1 (MQ=255)
                                                                                                                                    aGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCt                >  2:141350/1‑151 (MQ=255)
                                                                                                                                    aGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCt                >  1:786055/1‑151 (MQ=255)
                                                                                                                                        aTGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCttt              >  1:951278/1‑149 (MQ=255)
                                                                                                                                        aTGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCttt              >  1:461081/1‑149 (MQ=255)
                                                                                                                                        aTGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCtt               >  2:572404/1‑148 (MQ=255)
                                                                                                                                             gggCCTAAAATAAAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCa        <  1:770565/150‑1 (MQ=255)
                                                                                                                                                   gAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTg  >  2:201994/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
GAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTG  >  NC_000913/1222087‑1222382

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