Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,860,6030AT64.7% ‑4.8 / 27.7 34V103E (GTA→GAA) rpoSRNA polymerase sigma factor RpoS
Reads supporting (aligned to +/- strand):  ref base A (2/10);  new base T (22/0);  total (24/10)
Fisher's exact test for biased strand distribution p-value = 5.03e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.51e-05
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTC  >  NZ_CP009273/2860519‑2860691
                                                                                    |                                                                                        
aCCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTaccacc                                                                                     <  2:758675/90‑1 (MQ=255)
 ccGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGGAAATTTTTccccca                                                                                    >  1:431599/1‑90 (MQ=255)
   gcgcgGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGGAATTTTTtccccccaa                                                                                  >  1:342000/1‑87 (MQ=255)
   gcgcgGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGGAATTTTTTCCACCCGa                                                                                  >  2:134912/1‑90 (MQ=255)
   gcgcgGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTtccacccaa                                                                                  >  1:550608/1‑87 (MQ=255)
   gcgcgGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTTCCACCAGa                                                                                  >  2:558857/1‑90 (MQ=255)
   gcgcgGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTTCCACCAGa                                                                                  >  2:51324/1‑90 (MQ=255)
   gcgcgGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTTCCACCAGa                                                                                  >  2:622811/1‑90 (MQ=255)
   gcgcgGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTTCCACCAGa                                                                                  >  2:699969/1‑90 (MQ=255)
     gcgGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTTCCACCAGACg                                                                                >  2:656073/1‑90 (MQ=255)
     gcgGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACg                                                                                <  1:51324/90‑1 (MQ=255)
     gcgGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACg                                                                                <  2:666908/90‑1 (MQ=255)
          tCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGt                                                                           <  1:230581/90‑1 (MQ=255)
          tCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGt                                                                           <  1:897590/90‑1 (MQ=255)
               cccAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTActc                                                                      <  2:610042/90‑1 (MQ=255)
                ccAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAAATTTTTCCACCAGACGCCAATTActct                                                                     >  1:188352/1‑90 (MQ=255)
                  aGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTTCCCCCCAACGCAAGTTACTCTCg                                                                   >  2:755076/1‑90 (MQ=255)
                  aGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTTCCACCAAACGCAAGTTTCTCTCg                                                                   >  2:397092/1‑90 (MQ=255)
                  aGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTTCCACCAAACGCAAAGTACTCTCg                                                                   >  1:387400/1‑90 (MQ=255)
                       gCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGGAATTTTTTCCCCCACACGCAAGTTAATTTCcatcat                                                              >  1:254283/1‑90 (MQ=255)
                       gCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTTCCACCAGACGCAAGtt                                                                          >  2:562881/1‑78 (MQ=255)
                       gCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGtt                                                                          <  1:562881/78‑1 (MQ=255)
                            ctTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGGAATTTTTACCCCCAGACGCAAGTTACTCTCGGTCATcccgg                                                         >  2:252547/1‑87 (MQ=255)
                            ctTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTTCCACCAGACGCAAGTTATTCTCCATCATCcggc                                                         >  2:294113/1‑90 (MQ=255)
                            ctTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTTCCACCAGACGCAAGTTACTCTCGGTCATCccgc                                                         >  2:15098/1‑90 (MQ=255)
                            ctTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTTCCACCAGACGCAAGTTACTCTCGATCATCcggc                                                         >  2:73969/1‑90 (MQ=255)
                            ctTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAACTTTTTCCACCCAACGCAAGTTAGCCCCCCTCATCcggc                                                         >  2:242842/1‑90 (MQ=255)
                             ttCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGGAATTTTTTTCACCAAACGCAAAGTACTCTCGATCATCcggcg                                                        >  1:763480/1‑90 (MQ=255)
                                    aGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTTCCACCCGACGCAAGGTACTCTCGATCATCCGGCGGGGggag                                                 >  1:269350/1‑90 (MQ=255)
                                      gTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTTCCACCCGACGCAAGTTACTCTCGATCATCCGGCGGCGGGAggg                                               >  2:152732/1‑89 (MQ=255)
                                        ccAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAAATTTTTTCACCCAACCCCAATTACTCTCCATCAACCCGCGGCGGGAGGCGa                                             >  1:872608/1‑90 (MQ=255)
                                         cAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGAc                                            <  2:872608/90‑1 (MQ=255)
                                                gCCAGACCACGATTGCCATAACGGCGGGCAATTTTTTCCACCCGACGCAAGTTACTCTCCATCATCCGGCGGCGAGAGGGGAAATCTccc                                     >  2:474868/1‑89 (MQ=255)
                                                   aGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGc                                  <  1:656073/90‑1 (MQ=255)
                                                             tGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCCCGCAGGGCGCGAc                        >  2:583196/1‑90 (MQ=255)
                                                                              aTTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAActtct       <  2:39633/90‑1 (MQ=255)
                                                                                   tACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGTGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACttcttcttc  <  2:170096/90‑1 (MQ=255)
                                                                                    |                                                                                        
ACCGCGCGGATCAGCCCCAGGTTGCCCTCTTCGATAAGGTCCAGCAACGCCAGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTC  >  NZ_CP009273/2860519‑2860691

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