Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,797,3730T.64.1% 16.7 / 71.2 39coding (513/960 nt)nrdFclass 1b ribonucleoside‑diphosphate reductase subunit beta
Reads supporting (aligned to +/- strand):  ref base T (6/8);  new base . (6/19);  total (12/27)
Fisher's exact test for biased strand distribution p-value = 2.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.59e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATCGCGGTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAG  >  NZ_CP009273/2797285‑2797458
                                                                                        |                                                                                     
aTCGCGGTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTAtt                                                                                      >  1:558840/1‑90 (MQ=255)
    cgGTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTAttttcc                                                                                  <  2:756394/90‑3 (MQ=255)
        gatgatCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGc                                                                              >  2:427157/1‑90 (MQ=255)
             tCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCgg                                                                         <  1:146578/90‑1 (MQ=255)
              ccGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGaa                                                                       >  2:413139/1‑90 (MQ=255)
               cGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGaaa                                                                      >  2:376746/1‑90 (MQ=255)
                gCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAg                                                                     >  1:320313/1‑90 (MQ=255)
                 cTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGc                                                                    <  2:12669/90‑1 (MQ=255)
                 cTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGc                                                                    <  2:279867/90‑1 (MQ=255)
                 cTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGc                                                                    <  1:376746/90‑1 (MQ=255)
                  tGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGc                                                                    <  2:613488/90‑1 (MQ=255)
                  tGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCt                                                                   >  2:331360/1‑90 (MQ=255)
                        aGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCaa                                                             >  2:574899/1‑90 (MQ=255)
                         gAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAAt                                                            <  1:739812/90‑1 (MQ=255)
                         gAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAAt                                                            <  2:192797/90‑1 (MQ=255)
                           aaaTCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATa                                                           >  2:72735/1‑90 (MQ=255)
                              tCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACcgc                                                       <  1:331360/90‑1 (MQ=255)
                              tCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACcgc                                                       <  2:253141/90‑1 (MQ=255)
                                        gtTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATc                                             <  2:558840/90‑1 (MQ=255)
                                        gtTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATc                                             <  2:274572/90‑1 (MQ=255)
                                           tCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATc                                             >  1:574025/1‑87 (MQ=255)
                                           tCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATc                                             <  2:574025/87‑1 (MQ=255)
                                                              tATTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAg                       <  1:534428/90‑1 (MQ=255)
                                                               aTTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAg                       <  1:821002/90‑1 (MQ=255)
                                                               aTTCCGGTTTCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGc                      <  2:308346/90‑1 (MQ=255)
                                                                     gTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTcc                 <  2:385177/90‑1 (MQ=255)
                                                                      tttCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCAc               <  1:268657/90‑1 (MQ=255)
                                                                       ttCTGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACg              <  1:321099/90‑1 (MQ=255)
                                                                         cTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACgg             >  2:231661/1‑90 (MQ=255)
                                                                          tGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCTCGATGAAGCAGTCCACGGtt           <  2:527808/90‑1 (MQ=255)
                                                                          tGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGtt           <  2:112561/90‑1 (MQ=255)
                                                                          tGGCTGCCGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGtt           <  2:431502/90‑1 (MQ=255)
                                                                              tGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTtact        >  1:322495/1‑90 (MQ=255)
                                                                              tGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTtact        >  1:764107/1‑90 (MQ=255)
                                                                                 cGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACAt    <  1:740501/90‑1 (MQ=255)
                                                                                 cGATGTATTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACAt    <  1:673630/90‑1 (MQ=255)
                                                                                  gATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACAt    <  1:152420/90‑1 (MQ=255)
                                                                                  gATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACAt    <  1:382873/90‑1 (MQ=255)
                                                                                    tGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAg  <  1:103446/90‑1 (MQ=255)
                                                                                     gTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTta          <  1:810702/81‑1 (MQ=255)
                                                                                     gTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTta          >  2:810702/1‑81 (MQ=255)
                                                                                        |                                                                                     
ATCGCGGTGATGATCCGCTGAAAAAGAAAATCGCCAGTGTGTTTCTTGAATCTTTTTTGTTCTATTCCGGTTTCTGGCTGCCGATGTATTTTTCCAGCCGCGGAAAGCTGACCAATACCGCGGACCTGATCCGTCTGATTATCCGCGATGAAGCAGTCCACGGTTACTACATAG  >  NZ_CP009273/2797285‑2797458

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: