Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913300,7110GC54.3% 13.1 / 41.3 35P75R (CCG→CGG) paoCPaoABC aldehyde oxidoreductase, Moco‑containing subunit
Reads supporting (aligned to +/- strand):  ref base G (8/8);  new base C (19/0);  total (27/8)
Fisher's exact test for biased strand distribution p-value = 5.47e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.72e-04

GTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGCGG  >  NC_000913/300579‑300853
                                                                                                                                    |                                                                                                                                              
gTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCgcggcgg                                                                                                                              >  2:301317/1‑150 (MQ=255)
      gCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGACGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAACCAGTACGCCCCGCGCTTTTTTCGCGGCGTCCGTAt                                                                                                                        >  2:640695/1‑151 (MQ=255)
          cTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGGACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAg                                                                                                                    >  2:163302/1‑151 (MQ=255)
           tACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCCGCGCTTTTTGCGCGGCGGCCGTATCAAgg                                                                                                                   >  2:330154/1‑151 (MQ=255)
                       ggCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCCGGACGCCCCGCGCTTTTTTCGCGGCGTCCGTATCAAGGGGGGGGAGGCgg                                                                                                       >  2:334741/1‑150 (MQ=255)
                       ggCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGAAAGCCCGTACGCCCCGCGCCTTTTTTGCGGCGGCCCGATCAAgggggggggggggg                                                                                                       >  1:204155/1‑143 (MQ=255)
                       ggCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGACAGCCAGTACGCCCGGCGCTTTTTTCGCGGCGGCCCTATCAAGGGGGGGGGGGCgg                                                                                                       >  2:289576/1‑150 (MQ=255)
                        gCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCCGCGCCTTTTTCGCGGCGTCCGTATTCAGGGCGGGGAGGCGt                                                                                                       >  2:5981/1‑150 (MQ=255)
                         ccTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCTGAGTGCCCCGGCGTTACTGGCGGTTATGAAAGCCAGTACGCCCCGCGCTTTTTGCGCGGCGTGCGGATCAAGGGCGGGGAGGCggg                                                                                                      >  2:264079/1‑148 (MQ=255)
                         ccTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCCGTACGCCCCGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCggggggggggc                                                                                                      >  2:451643/1‑141 (MQ=255)
                            gatgatAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCCGTACGCCCCGCGCCTTTTTCGCGGCGTCCGGATCAAAGGCGGGGAGGCGTCCctt                                                                                                  >  1:308094/1‑151 (MQ=255)
                            gatgatAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTTCGCGGCGTCCGTATCAAGGGCGGTGAGGCGGCCttt                                                                                                  >  2:482597/1‑151 (MQ=255)
                                   gTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTG‑CCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCCGCGCCTTTTTCGCGGCGGCCGTATCAAGGGCGGTGAGGGGTCCTTTTGCAATgg                                                                                          >  2:132760/1‑151 (MQ=255)
                                     gCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAAAGACAGCCAGGACGCCCCGCGCCTTTTTCGCGGCGGCCGGCTCCAGGGCGGGGGGGCGTCCCTTGGCCATgg                                                                                          >  1:639677/1‑150 (MQ=255)
                                        cAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCCGCGCTTTTTTCGCGGGGGCCGGATCAAGGGGGGGGAGGGGTCCTTTTGGAATGGCggg                                                                                      >  2:32083/1‑150 (MQ=255)
                                        cAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCCGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGGGAGGCGTCCTTTGGCAATGGCggg                                                                                      >  2:67663/1‑150 (MQ=255)
                                          aTAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCCGCGCTTTTTTCGCGGCGTCCGTATCAAGGGGGGGGGGGCGGCCTTTGGCAATGGCgggaa                                                                                    >  2:277807/1‑148 (MQ=255)
                                          aTAGTGGGGCCGCCTAACAGCCTGGCGGTGGTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGACAGCCCGTACGCCCCGCGCTTTTTGCGCGGCGTCCGGATCAAGGGCGGGGGGGGGGCCCTTGGCAATGGGGGaa                                                                                     >  2:265402/1‑150 (MQ=255)
                                                  gccgccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATTATACCCCGTACGCCCGGCGCTTTTTTCGCGGCGGCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATTGCGGCAccgccggca                                                                            >  1:627758/1‑148 (MQ=255)
                                                   ccgccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCCGCGCTTTTTTCGCGGCGTCCGTATCAAGGGCGGGGGGGGGGCCTTTTGCAATGGCGGGAACGACGgtat                                                                           >  1:686292/1‑151 (MQ=255)
                                                        tAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCa                                                                      <  2:691532/151‑1 (MQ=255)
                                                                     gtgtTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCCGCGCTTTTTTCGCGGCGTCCGTATCAAGGGGGGGGAGGCGGCCTTTTGGAATGGCGGGAACGACGATATAGCCATAGGCGGCGGTgg                                                         >  2:96904/1‑151 (MQ=255)
                                                                     gtgtTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCCGCGCTTTTTGCGCGGCGTCCGGATCAAGGGCGGGGAGGGGTCCCTTGGGAATGGCGGGACCGGCGATATAAcccaaggggggggtgg                                                         >  1:527528/1‑137 (MQ=255)
                                                                         ttttGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGc                                                     <  1:272684/151‑1 (MQ=255)
                                                                         ttttGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAAGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGc                                                     <  2:239744/151‑1 (MQ=255)
                                                                                 ccTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGGGGGGAGGCGGCCTTTGGCAATGGCGGGACCGACGATATAGCCATAGGCGGCGGTGGGGGCTTCTTCAt                                             >  2:102114/1‑151 (MQ=255)
                                                                                 ccTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCAt                                             <  2:349689/151‑1 (MQ=255)
                                                                                 ccTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGGGAGGCGTCCTTTGGCAATGGCGGGAACGACGATATAGCCATAGGCGGCGTTGGGGGGTTCTTCAt                                             >  2:207058/1‑151 (MQ=255)
                                                                                             gCCCCGGCGTTACTGGCGGTAAAGACAGCCAGTACGCCCCGCGCTTTTTTCGCGGCGTCCGCATCAAGGGCGGGGAGGCGGCCTTTGGCAATGGGGGAACCGACGATATAGCCATAGGCGGGGTTGGGGGGTTCTTCATGCCCTTCGTAgg                                 >  2:532785/1‑151 (MQ=255)
                                                                                              ccccGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAgg                                 <  2:487300/150‑1 (MQ=255)
                                                                                                  ggCGTTACTGGCGGTAATGACAGCCAGTACGCCCCGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGGACCGACGATATAGCCCTAGGCGGCGGTGGGGGGTTCTTCATGGCATTCGTAGGCGGAg                            >  1:256523/1‑151 (MQ=255)
                                                                                                   gCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAg                            <  1:207056/150‑1 (MQ=255)
                                                                                                      ttACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGcgt                         <  2:714187/150‑1 (MQ=255)
                                                                                                             cGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGcgtgccgtgc                  <  1:225870/150‑1 (MQ=255)
                                                                                                                            gTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGCgg  >  2:520562/1‑151 (MQ=255)
                                                                                                                                    |                                                                                                                                              
GTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGCGG  >  NC_000913/300579‑300853

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