Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,389,5300TC65.0% 30.6 / 44.9 40K134E (AAA→GAA) nuoINADH‑quinone oxidoreductase subunit NuoI
Reads supporting (aligned to +/- strand):  ref base T (9/5);  new base C (13/13);  total (22/18)
Fisher's exact test for biased strand distribution p-value = 5.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TTGACGTCGATAGGCTTGGCTTCGTTCTCTGCTTCGCCCTTATCTTTGCCGTCGATTGCCATACCTGCCATCCGGTAGAAGTTATATTCCGGGTATTTGCCCGGACCGGAGATCAGCAGATCCTCTTTCTCGTAAACCAGATCCTGGCGCTTGTATTCCCCCATTTCGAAATCCGGGGTTAACTGAATCGCCGTGGTCGGACAGGCTTCTTCGCACAGACCACAGAAAATGCAGCGTGAGAAGTTGATGCGGAAAAAT  >  NZ_CP009273/2389403‑2389660
                                                                                                                               |                                                                                                                                  
ttGACGTCGATAGGCTTGGCTTCGTTCTCTGCTTCGCCCTTATCTTTGCCGTCGATTGCCATACCTGCCATCCGGTAGAAGTTATATTCCGGGTATTTGCCCGGACCGGAGATCAGCAGATCCTCTTTCTCGTAAAc                                                                                                                           >  2:629592/1‑137 (MQ=255)
ttGACGTCGATAGGCTTGGCTTCGTTCTCTGCTTCGCCCTTATCTTTGCCGTCGATTGCCATACCTGCCATCCGGTAGAAGTTATATTCCGGGTATTTGCCCGGACCGGAGATCAGCAGATCCTCTTTCTCGTAAAc                                                                                                                           >  2:238263/1‑137 (MQ=255)
ttGACGTCGATAGGCTTGGCTTCGTTCTCTGCTTCGCCCTTATCTTTGCCGTCGATTGCCATACCTGCCATCCGGTAGAAGTTATATTCCGGGTATTTGCCCGGACCGGAGATCAGCAGATCCTCTTCCTCGTAAAc                                                                                                                           >  1:703272/1‑137 (MQ=255)
        gATAGGCTTGGCTTCGTTCTCTGCTTCGCCCTTATCTTTGCCGTCGATTGCCATACCTGCCATCCGGTAGAAGTTATATTCCGGGTATTTGCCCGGACCGGAGATCAGCAGATCCTCTTCCTCGTAAACCAGATCCTgg                                                                                                                 >  2:605949/1‑139 (MQ=255)
            ggcttggcttCGTTCTCTGCTTCGCCCTTATCTTTGCCGTCGATTGCCATACCTGCCATCCGGTAGAAGTTATATTCCGGGTATTTGCCCGGACCGGAGATCAGCAGATCCTCTTTCTCGTAAACCAGATCCTGGCGCt                                                                                                             <  1:238263/139‑1 (MQ=255)
            ggcttggcttCGTTCTCTGCTTCGCCCTTATCTTTGCCGTCGATTGCCATACCTGCCATCCGGTAGAAGTTATATTCCGGGTATTTGCCCGGACCGGAGATCAGCAGATCCTCTTTCTCGTAAACCAGATCCTGGCGCt                                                                                                             <  1:629592/139‑1 (MQ=255)
            ggcttggcttCGTTCTCTGCTTCGCCCTTATCTTTGCCGTCGATTGCCATACCTGCCATCCGGTAGAAGTTATATTCCGGGTATTTGCCCGGACCGGAGATCAGCAGATCCTCTTCCTCGTAAACCAGATCCTGgcgc                                                                                                              <  2:703272/138‑1 (MQ=255)
            ggcttggcttCGTTCTCTGCTTCGCCCTTATCTTTGCCGTCGATTGCCATACCTGCCATCCGGTAGAAGTTATATTCCGGGTATTTGCCCGGACCGGAGATCAGCAGATCCTCTTCCTCGTAAACCAGATCCTGgcgc                                                                                                              <  1:535709/138‑1 (MQ=255)
            ggcttggcttCGTTCTCTGCTTCGCCCTTATCTTTGCCGTCGATTGCCATACCTGCCATCCGGTAGAAGTTATATTCCGGGTATTTGCCCGGACCGGAGATCAGCAGATCCTCTTCCTCGTAAACCAGATCCTGgcgc                                                                                                              <  2:288326/138‑1 (MQ=255)
                 ggcttCGTTCTCTGCTTCGCCCTTATCTTTGCCGTCGATTGCCATACCTGCCATCCGGTAGAAGTTATATTCCGGGTATTTGCCCGGACCGGAGATCAGCAGATCCTCTTTCTCGTAAACCAGATCCTGGCGCTTGTAt                                                                                                        >  2:844086/1‑139 (MQ=255)
                 ggcttCGTTCTCTGCTTCGCCCTTATCTTTGCCGTCGATTGCCATACCTGCCATCCGGTAGAAGTTATATTCCGGGTATTTGCCCGGACCGGAGATCAGCAGATCCTCTTTCTCGTAAACCAGATCCTGGCGCTTGTAt                                                                                                        >  1:466284/1‑139 (MQ=255)
                                   gCCCTTATCTTTGCCGTCGATTGCCATACCTGCCATCCGGTAGAAGTTATATTCCGGGTATTTGCCCGGACCGGAGATCAGCAGATCCTCTTTCTCGTAAACCAGATCCTGGCGCTTGTATTCCCCCATTTCGAAATcc                                                                                      <  2:466284/139‑1 (MQ=255)
                                       ttATCTTTGCCGTCGATTGCCATACCTGCCATCCGGTAGAAGTTATATTCCGGGTATTTGCCCGGACCGGAGATCAGCAGATCCTCTTCCTCGTAAACCAGATCCTGGCGCTTGTATTCCCCCATTTCGAAATCCgggg                                                                                  >  2:863509/1‑139 (MQ=255)
                                              tGCCGTCGATTGCCATACCTGCCATCCGGTAGAAGTTATATTCCGGGTATTTGCCCGGACCGGAGATCAGCAGATCCTCTTTCTCGTAAACCAGATCCTGGCGCTTGTATTCCCCCATTTCGAAATCCGGGGTTAACt                                                                            >  1:390335/1‑138 (MQ=255)
                                              tGCCGTCGATTGCCATACCTGCCATCCGGTAGAAGTTATATTCCGGGTATTTGCCCGGACCGGAGATCAGCAGATCCTCTTCCTCGTAAACCAGATCCTGGCGCTTGTATTccccc                                                                                                  <  1:433525/116‑1 (MQ=255)
                                              tGCCGTCGATTGCCATACCTGCCATCCGGTAGAAGTTATATTCCGGGTATTTGCCCGGACCGGAGATCAGCAGATCCTCTTCCTCGTAAACCAGATCCTGGCGCTTGTATTccccc                                                                                                  >  2:433525/1‑116 (MQ=255)
                                              tGCCGTCGATTGCCATACCTGCCATCCGGTAGAAGTTATATTCCGGGTATTTGCCCGGACCGGAGATCAGCAGATCCTCTTCCTCGTAAACCAGATCCTGGCGCTTGTATTCCCCCATTTCGAAATCCGGGGTTAACt                                                                            >  2:96185/1‑138 (MQ=255)
                                                ccGTCGATTGCCATACCTGCCATCCGGTAGAAGTTATATTCCGGGTATTTGCCCGGACCGGAGATCAGCAGATCCTCTTCCTCGTAAACCAGATCCTGGCGCTTGTATTCCCCCATTTCGAAATCCGGGGTTAACTGaa                                                                         <  1:658795/139‑1 (MQ=255)
                                                         gCCATACCTGCCATCCGGTAGAAGTTATATTCCGGGTATTTGCCCGGACCGGAGATCAGCAGATCCTCTTCCTCGTAAACCAGATCCTGGCGCTTGTATTCCCCCATTTCGAAATCCGGGGTTAACTGAATCGCCGTg                                                                 <  2:295698/138‑1 (MQ=255)
                                                                   ccATCCGGTAGAAGTTATATTCCGGGTATTTGCCCGGACCGGAGATCAGCAGATCCTCTTCCTCGTAAACCAGATCCTGGCGCTTGTATTCCCCCATTTCGAAATCCGGGGTTAACTGAATCGCCGTGGTCGGACAgg                                                       <  2:536943/138‑1 (MQ=255)
                                                                     aTCCGGTAGAAGTTATATTCCGGGTATTTGCCCGGACCGGAGATCAGCAGATCCTCTTCCTCGTAAAcc                                                                                                                          >  1:823634/1‑69 (MQ=255)
                                                                     aTCCGGTAGAAGTTATATTCCGGGTATTTGCCCGGACCGGAGATCAGCAGATCCTCTTCCTCGTAAAcc                                                                                                                          <  2:823634/69‑1 (MQ=255)
                                                                                   atatTCCGGGTATTTGCCCGGACCGGAGATCAGCAGATCCTCTTTCTCGTAAACCAGATCCTGGCGCTTGTATTCCCCCATTTCGAAATCCGGGGTTAACTGAATCGCCGTGGTCGGACAGGCTTCTTCGCACAGAcc                                       <  2:572920/138‑1 (MQ=255)
                                                                                   atatTCCGGGTATTTGCCCGGACCGGAGATCAGCAGATCCTCTTCCTCGTAAACCAGATCCTGGCGCTTGTATTCCCCCATTTCGAAATCCGGGGTTAACTGAATCGCCGTGGTCGGACAGGCTTCTTCGCACAGAcc                                       <  1:863509/138‑1 (MQ=255)
                                                                                   atatTCCGGGTATTTGCCCGGACCGGAGATCAGCAGATCCTCTTCCTCGTAAACCAGATCCTGGCGCTTGTATTCCCCCATTTCGAAATCCGGGGTTAACTGAATCGCCGTGGTCGGACAGGCTTCTTCGCACAGAcc                                       <  1:58562/138‑1 (MQ=255)
                                                                                          gggTATTTGCCCGGACCGGAGATCAGCAGATCCTCTTCCTCGTAAACCAGATCCTGGCGCTTGTATTCCCCCATTTCGAAATCCGGGGTTAACTGAATCGCCGTGGTCGGACAGGCTTCTTCGCACAGACCACAGaaaa                               >  1:619060/1‑139 (MQ=255)
                                                                                           ggTATTTGCCCGGACCGGAGATCAGCAGATCCTCTTTCTCGTAAACCAGATCCTGGCGCTTGTATTCCCCCATTTCGAAATCCGGGGTTAACTGAATCGCCGTGGTCGGACAGGCTTCTTCGCACAGACCACAGAAAAt                              >  1:843656/1‑139 (MQ=255)
                                                                                           ggTATTTGCCCGGACCGGAGATCAGCAGATCCTCTTTCTCGTAAACCAGATCCTGGCGCTTGTATTCCCCCATTTCGAAATCCGGGGTTAACTGAATCGCCGTGGTCGGACAGGCTTCTTCGCACAGACCACAGAAAAt                              >  1:484793/1‑139 (MQ=255)
                                                                                                       gaccggaGATCAGCAGATCCTCTTCCTCGTAAACCAGATCCt                                                                                                                   <  2:327315/42‑1 (MQ=255)
                                                                                                       gaccggaGATCAGCAGATCCTCTTCCTCGTAAACCAGATCCt                                                                                                                   >  1:327315/1‑42 (MQ=255)
                                                                                                            gagaTCAGCAGATCCTCTTTCTCGTAAACCAGATCCTGGCGCTTGTATTCCCCCATTTCGAAATCCGGGGTTAACTGAATCGCCGTGGTCGGACAGGCTTCTTCGCACAGACCACAGAAAATGCAGCGTGAGAAGTTGa             >  2:90492/1‑139 (MQ=255)
                                                                                                            gagaTCAGCAGATCCTCTTCCTCGTAAACCAGATCCTGGCGCTTGTATTCCCCCATTTCGAAATCCGGGGTTAACTGAATCGCCGTGGTCGGACAGGCTTCTTCGCACAGACCACAGAAAATGCAGCGTGAGAAGTTGa             >  2:196738/1‑139 (MQ=255)
                                                                                                                  agcagATCCTCTTCCTCGTATACCAGATCCTGGCGCTTGTATTCCCCCATTTCGAAATCCGGGGTTAACTGAATCGCCGTGGTCGGACAGGCTTCTTCGcaca                                           >  1:755250/1‑103 (MQ=255)
                                                                                                                  agcagATCCTCTTCCTCGTATACCAGATCCTGGCGCTTGTATTCCCCCATTTCGAAATCCGGGGTTAACTGAATCGCCGTGGTCGGACAGGCTTCTTCGcaca                                           <  2:755250/103‑1 (MQ=255)
                                                                                                                  agcagATCCTCTTCCTCGTAAACCAGATCCTGGCGCTTGTATTCCCCCATTTCGAAATCCGGGGTTAACTGAATCGCCGTGGTCGGACAGGCTTCTTCGCACAGAcc                                       >  2:155706/1‑107 (MQ=255)
                                                                                                                  agcagATCCTCTTCCTCGTAAACCAGATCCTGGCGCTTGTATTCCCCCATTTCGAAATCCGGGGTTAACTGAATCGCCGTGGTCGGACAGGCTTCTTCGCACAGAcc                                       <  1:155706/107‑1 (MQ=255)
                                                                                                                  agcagATCCTCTTCCTCGTAAACCAGATCCTGGCGCTTGTATTCCCCCATTTCGAAATCCGGGGTTAACTGAATCGCCGTGGTCGGACAGGCTTCTTCGCACAGACCACAGAAAATGCAGCGTGAGAAGTTGATGCgg        >  1:766161/1‑138 (MQ=255)
                                                                                                                       aTCCTCTTTCTCGTAAACCAGATCCTGGCGCTTGTATTCCCCCATTTCGAAATCCGGGGTTAACTGAATCGCCGTGGTCGGACAGGCTTCTTCGCACAGACCACAGAAAATGCa                           <  1:139548/114‑1 (MQ=255)
                                                                                                                       aTCCTCTTTCTCGTAAACCAGATCCTGGCGCTTGTATTCCCCCATTTCGAAATCCGGGGTTAACTGAATCGCCGTGGTCGGACAGGCTTCTTCGCACAGACCACAGAAAATGCa                           >  2:139548/1‑114 (MQ=255)
                                                                                                                       aTCCTCTTCCTCGTAAACCAGATCCTGGCGCTTGTATTCCCCCATTTCGAAATCCGGGGTTAACTGAATCGCCGTGGTCGGACAGGCTTCTTCGCACAGACCACAGAAAATGCAGCGTGAGAAGTTGATGCGGAAAAAt  >  2:381289/1‑139 (MQ=255)
                                                                                                                               |                                                                                                                                  
TTGACGTCGATAGGCTTGGCTTCGTTCTCTGCTTCGCCCTTATCTTTGCCGTCGATTGCCATACCTGCCATCCGGTAGAAGTTATATTCCGGGTATTTGCCCGGACCGGAGATCAGCAGATCCTCTTTCTCGTAAACCAGATCCTGGCGCTTGTATTCCCCCATTTCGAAATCCGGGGTTAACTGAATCGCCGTGGTCGGACAGGCTTCTTCGCACAGACCACAGAAAATGCAGCGTGAGAAGTTGATGCGGAAAAAT  >  NZ_CP009273/2389403‑2389660

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