Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,392,6170A.41.7% 55.4 / 52.6 38coding (425/1530 nt)nuoMNADH:quinone oxidoreductase subunit M
Reads supporting (aligned to +/- strand):  ref base A (7/14);  new base . (10/5);  total (19/19)
Fisher's exact test for biased strand distribution p-value = 8.95e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.30e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCCAGAAGAAGAAGAACAGGAACATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAACATAGTACCGCCAGCACACCGAGCAGA  >  NC_000913/2392476‑2392757
                                                                                                                                             |                                                                                                                                            
attgCAATCAACACCACCAGGCCGCTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGACGTGCTACTCCTCTAACCCACCCCGCGATAATCCCTCCACAGTGCGATAAGGAAGTACACCGGCACCAGCATCATTTCCCagaa                                                                                                                                             <  2:852631/140‑1 (MQ=255)
aTGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCCagaa                                                                                                                                             <  1:1013396/143‑1 (MQ=255)
   gCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCCagaagaa                                                                                                                                          <  2:184827/143‑1 (MQ=255)
     aaTCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCCagaagaaga                                                                                                                                        >  1:910687/1‑143 (MQ=255)
       tCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGCTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCCagaagaagaag                                                                                                                                      <  1:64653/143‑1 (MQ=255)
               aCCAGGCCACTCGCCGGGGTGTAAATGAAGAACGTGGTTGCCGCCGTGATACGCGTTTTAGCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCCAGAAGAAGAAGAACAGGaa                                                                                                                              <  2:932200/143‑1 (MQ=255)
                 cAGGCAACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCGGGAAGTACAGCGGCACCAGCATCATTTCCCCGCAGAAGCAGAACAGGAACa                                                                                                                            >  2:801083/1‑143 (MQ=255)
                       aCTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCCAGAAGAAGAAGAACAGGAACATGTCGa                                                                                                                      >  1:831169/1‑143 (MQ=255)
                        cccgcgTGGGTGTAAATGAAGAACTTGGGTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGGGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCCAGAAGAAGAAGAACAGGAACATGTCGAt                                                                                                                     <  1:1084035/141‑1 (MQ=255)
                            ccTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTGTGGGGCTATACGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCC‑‑‑AGAAGAAGAACAGGAACATGTCGATGGCAAgg                                                                                                              <  1:1193492/143‑1 (MQ=255)
                                   ttAAATGTTGAAATTGAGTGTCGCCGTGATACCCGTTTTACCGTCAGCGGCTTTATGCCCCCACAGTGCGATCAGAAAGTACATCGCCACCAGCATCATTTCCC‑‑‑AGAAGAAGAACAGGAACATGTCGATGGCAAGGACCACGc                                                                                                       <  2:281884/142‑1 (MQ=255)
                                         gaagaaCTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTGTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCC‑‑‑AGAAGAAGAACAGGAACATGTCGATGGCAAGGAACACGCCGATaa                                                                                                 <  2:1003255/143‑1 (MQ=255)
                                           agaaCTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCCAGAAGAAGAAGAACAGGAACATGTCGATGGCAAGGAACACGCCGATa                                                                                                  <  2:910687/143‑1 (MQ=255)
                                            gaaCTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCC‑‑‑AGAAGAAGAACAGGAACAGGTCGATGGCAAGCAACCCGCCGACAAcgc                                                                                              >  2:267668/1‑143 (MQ=255)
                                                  ggTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGCTCAGGAAGTACATCGGCACCAGCATCATTTCCC‑‑‑AGAAGAAGAACAGGAACATGTCGAGGGCAACGAACACGCCGATCACCCCGCCCa                                                                                        >  1:484729/1‑143 (MQ=255)
                                                    ttGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCCAGAAGAAGAAGAACAGGAACATGTCGATGGCAAGGAACACGCCGATAACGCCGccc                                                                                         <  1:295096/143‑1 (MQ=255)
                                                         gccgTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCCAGAAGAAGAAGAACAGGAACATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGAt                                                                                    >  1:74226/1‑143 (MQ=255)
                                                          ccgTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCC‑‑‑AGTAGAAGACCAGGACCAGGGCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCcac                                                                                >  1:1072853/1‑143 (MQ=255)
                                                                tgcgcgTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCC‑‑‑AGAAGAAGAACAGGAACATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAgg                                                                          <  1:738035/141‑1 (MQ=255)
                                                                      ttttACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCCGGAAGTACATCGGCACCAGCAGCATTTCCC‑‑‑AGACGAAGAACAGGCACATGTCGATGGCACGGAACACGCCGATAACGCCGCCCAGGCTCCGCATCAGGtagcgg                                                                    >  1:512961/1‑138 (MQ=255)
                                                                           ccGTCAGAGGCTTTATGCCCCCACAGTGCGACCAGGAAGTACATCGGCACCAGCATCCTTGCGCTGGAGAAGAAGCACAGGAACATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCGGCTCGAggcg                                                                  >  2:300950/1‑141 (MQ=255)
                                                                            cGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCC‑‑‑AGAAGAAGAACAGGAACATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGTGGTGgatg                                                              >  1:155714/1‑141 (MQ=255)
                                                                                      tttATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCC‑‑‑AGAAGAAGAACAGGAACATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGAACCCAATCCGGCTAAGGTGGAAGAAGGCCTGat                                                    >  2:879720/1‑143 (MQ=255)
                                                                                      tgtacgCCCGCACAGGTCGATCAGGAAGAACAGCGGCACCAGCATCATGCCCCAGAAGAAGAAGAACAGGAACAGGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCt                                                       <  2:1072155/138‑1 (MQ=255)
                                                                                      ctcccgCCCCCCCGTTGCGTGCAGGAAGTACATCGGCACTAGCATCATTGCCCAGAAGAAGAAGAACAGGAACATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTTAGGTGGAAGAAGCCCt                                                       <  2:873988/138‑1 (MQ=255)
                                                                                           gCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCCAGAAGAAGAAGAACAGGAACATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGatat                                                  <  1:1069192/143‑1 (MQ=255)
                                                                                            cccccACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCCAGAAGAAGAAGAACAGGAACATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGatatg                                                 >  1:645669/1‑142 (MQ=255)
                                                                                             ccccACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCCAGAAGAAGAAGAACAGGAACATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATAttt                                                >  1:909350/1‑143 (MQ=255)
                                                                                                   ggTGCGATCGGGAAGTACATCGGCACCAGCCTTTTTGCCCAGAAGAAGAAGAACAGGAACAGGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGAt                                          <  2:246822/142‑1 (MQ=255)
                                                                                                      gCGATCAGGAAGTACATCGGCACCAGCATCATTTCCC‑‑‑AGAAGAAGAACAGGAACATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCttt                                    >  1:95068/1‑143 (MQ=255)
                                                                                                        gATCAGGAAGTACATCGGCACCAGCATCATTTCCC‑‑‑AGAAGAAGAACAGGAACATGTCGCTGACAAAGAACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTcc                                  >  1:169668/1‑143 (MQ=255)
                                                                                                        gATCAGGAAGTACATCGGCACCAGCATCATTTCCC‑‑‑AGAAGAAGAACAGGAACATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTcc                                  >  1:844502/1‑143 (MQ=255)
                                                                                                            aGGAAGTACATCGGCACCAGCATCATTTCCCAGAAGAAGAAGAACAGGAACATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTTAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCa                                 <  1:330997/143‑1 (MQ=255)
                                                                                                            aGGAAGTACATCGGCACCAGCATCATTTCCC‑‑‑AGAAGAAGAACAGGAACATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGa                              >  2:39417/1‑143 (MQ=255)
                                                                                                                   acacCGGCACCAGCATCACTTCCCAGAAGAAGAAGCACAGGACCATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAACAt                          <  1:67039/143‑1 (MQ=255)
                                                                                                                    caccGGCACCAGCATCATTTCCCAGAAGAAGAAGAACAGGAACATGTCGATGGCAAGGAACACGCCAATTACGCCGCCCAGGATCCAAATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAACATa                         <  1:1011842/140‑1 (MQ=255)
                                                                                                                    cATCGGCACCACCATCATTTCCCAGAAGAAGACGAACAGGACCATGTCGATGGCAAGGCACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAACATa                         <  1:849904/143‑1 (MQ=255)
                                                                                                                     aTCGGCACCAGCATCATTTCCCAGAAGAAGAAGAACAGGAACATGTCGATGGCAAGGAACACGCCGATACCGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAACATAg                        >  1:1115503/1‑143 (MQ=255)
                                                                                                                     aTCGGCACCAGCATCACTTCCCAGAAGAAGAAGAACAGGAACATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAACATAg                        <  1:837496/143‑1 (MQ=255)
                                                                                                                           aCCAGCATCATTTCCC‑‑‑AGAAGAAGAACAGGAACATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAACATAGTACCGCCAg               <  1:357101/143‑1 (MQ=255)
                                                                                                                              aGCATCATTTCCCAGAAGAAGAAGAACAGGAACATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAACATAGTACCGCCAg               >  1:863030/1‑143 (MQ=255)
                                                                                                                                 atcatTCCCCAGAAGAAGAAGAACAGGAACATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGCTCCACATCAGGTTGAGGTGGAAGACGGCCTGATATTTTTCGATATGTATCCACGAAAGTAGTACCGCCAGcac            >  2:242372/1‑143 (MQ=255)
                                                                                                                                    atTTCCCAGAAGAAGAAGAACTGGAGCATGTCGATGGCGAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATATCTTTCCACGAGCATAGTACCGCCAGCACAcc         <  2:628824/143‑1 (MQ=255)
                                                                                                                                           agaagaagaagaACAGGAACATGTCGATGGCAAGGCACCCGCCGATAACGCCGCCCAGGATCCACATCACGTCGACGCGGAAGAAGCCCTGATACCCTCCGATCCCTCTCCACCAACAACGTACCCCCACCACACCGacccgc  >  2:587193/1‑138 (MQ=255)
                                                                                                                                             |                                                                                                                                            
ATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCCAGAAGAAGAAGAACAGGAACATGTCGATGGCAAGGAACACGCCGATAACGCCGCCCAGGATCCACATCAGGTTGAGGTGGAAGAAGCCCTGATATTTTTCGATCTCTTTCCACGAACATAGTACCGCCAGCACACCGAGCAGA  >  NC_000913/2392476‑2392757

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