Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,024,6090TG20.8% 43.1 / 11.0 24intergenic (+18/+26)rcsA/dsrBtranscriptional regulator of colanic acid capsular biosynthesis/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (7/12);  new base G (3/2);  total (10/14)
Fisher's exact test for biased strand distribution p-value = 6.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CAAGACCGTTTCATCGCATAAAGGTAATATTAAACGTAAGATCAAAACGCATAATAAACAGGTTATCTACCATGTCGTCCGACTGACGGATAATGTGACTAATGGTATTTTTGTCAACATGCGCTAACACATTCTGACTGGTGGTTTCCCACCAGTCAGGCTGAATAAGATTACTCTGCTTTCTCCACAAAGATACCGTCCTGATGCCCTGCTTCATTAAAGAACCAGATGCCGAGCGGGTAGTCTTCCAGCGAAACCAGGTACATTGTGCCTTCACTAAACT  >  NC_000913/2024465‑2024747
                                                                                                                                                |                                                                                                                                          
cAAGACCGTTTCATCGCATAAAGGTAATATTAAACATAAGATCAAAACGCATAATAAACAGGTTATCTACCATGTCGTCCGACTGACGGATAATGTGACTAATGGTATTTTTGTCAACATGCGCTAACACATTCTGACTGGTGGTTTcc                                                                                                                                        <  1:16310/149‑1 (MQ=255)
   gACCGTTTCATCGCATAAAGGTAATATTAAACGTAAGATCAAAACGCATAATAAACAGGTTATCTACCATGTCGTCCGACTGACGGATAATGTGACTAATGGTATTTTTGTCAACATGCGCTAACACATTCTGACTGGTGGTTTCccac                                                                                                                                     <  4:158351/149‑1 (MQ=255)
        tttCATCGCATAAAGGTAATATTAAACGTAAGATCAAAACGCATAATAAACAGGTTATCTACCATGTCGTCCGACTGACGGATAATGTGACTAATGGTATTTTTGTCAACATGCGCTAACACATTCTGACTGGTGGGAAACCACCAGTc                                                                                                                                <  7:84203/149‑1 (MQ=255)
                cATAAAGGTAATATTAAACGTAAGATCAAAACGCATAATAAACAGGTTATCTACCATGTCGTCCGACTGACGGATAATGTGACTAATGGTATTTTTGTCAACATGCGCTAACACATTCTGACTGGTGGTTTCCCACCAGTCAGGCTGaa                                                                                                                        >  4:51107/1‑149 (MQ=255)
                        tAATATTAAACGTAAGAGCAAAACGCATAATAAACAGGTTATCTACCATGTCGTCCGACTGACGGATAATGTGACTAATGGTATTTTTGTCAACATGCGCTAACACATTCTGACTGGTGGTTTCCCACCAGTCAGGCTGAATAAGATTa                                                                                                                <  3:51107/149‑1 (MQ=255)
                                    tAAGATCAAAACGCATAATAAACAGGTTATCTACCATGTCGTCCGACTGACGGATAATGTGACTAATGGTATTTTTGTCAACATGCGCTAACACATTCTGACTGGTGGTTTCCCACCAGTCAGGCTGAATAAGATTACTCTGCTTtctc                                                                                                    <  4:97814/149‑1 (MQ=255)
                                           aaaaCGCATAATAAACAGGTTATCTACCATGTCGTCCGACTGACGGATAATGTGACTAATGGTATTTTTGTCAACATGCGCTAACACATTCTGACTGGTGGTTTCCCACCAGTCAGGCTGAATAAGATTACTCTGCTTTCTCCACaaa                                                                                              <  5:24931/148‑1 (MQ=255)
                                                    aataaACAGGTTATCTACCATGTCGTCCGACTGACGGATAATGTGACTAATGGTATTTTTGTCAACATGCGCTAACACATTCTGACTGGTGGTTTc                                                                                                                                         >  1:95764/1‑96 (MQ=255)
                                                    aataaACAGGTTATCTACCATGTCGTCCGACTGACGGATAATGTGACTAATGGTATTTTTGTCAACATGCGCTAACACATTCTGACTGGTGGTTTc                                                                                                                                         <  2:95764/96‑1 (MQ=255)
                                                    aataaACAGGTTATCTACCATGTCGTCCGACTGACGGATAATGTGACTAATGGTATTTTTGTCAACATGCGCTAACACATTCTGACTGGTGGTTTCCCACCAGTCAGGCTGAATAAGATTACTCTGCt                                                                                                         >  8:98904/1‑128 (MQ=255)
                                                    aataaACAGGTTATCTACCATGTCGTCCGACTGACGGATAATGTGACTAATGGTATTTTTGTCAACATGCGCTAACACATTCTGACTGGTGGTTTCCCACCAGTCAGGCTGAATAAGATTACTCTGCt                                                                                                         <  7:98904/128‑1 (MQ=255)
                                                       aaaCAGGTTATCTACCATGTCGTCCGACTGACGGATAATGTGACTAATGGTATTTTTGTCAACATGCGCTAACACATCCTGACTGGTGGGAAACCACCAGTCAgaa                                                                                                                            >  2:157477/1‑104 (MQ=38)
                                                       aaaCAGGTTATCTACCATGTCGTCCGACTGACGGATAATGTGACTAATGGTATTTTTGTCAACATGCGCTAACACATCCTGACTGGTGGGAAACCACCAGTCAgaa                                                                                                                            <  1:157477/106‑3 (MQ=38)
                                                             gTTATCTACCATGTCGTCCGACTGACGGATAATGTGACTAATGGTATTTTTGTCAACATGCGCTAACACATTCTGACTGGTGGTTTCCCACCAGTCAGGCTGAATAAGATTACTCTGCTTTCTCCACAAAGATACCGTCCTGATGCCCt                                                                           >  7:63186/1‑149 (MQ=255)
                                                                aTCTACCATGTCGTCCGACTGACGGATAATGTGACTAATGGTATTTTTGTCAACATGCGCTAACACATTCTGACTGGTGGTTTCCCACCAGTCAGGCTGAATAAGATTACTCTGCTTTCTCCACAAAGATACCGTCCTGATGCCCTGc                                                                         <  8:63186/148‑1 (MQ=255)
                                                                            gtcCGACTGACGGATAATGTGACTAATGGTATTTTTGTCAACATGCGCTAACACATTCTGACTGGTGGTTTCCCACCAGTCAGGCTGAATAAGATTACTCTGCTTTCTCCACAAAGATACCGTCCTGATGCCCTGCTTCATTAAAGaa                                                             <  3:60871/148‑1 (MQ=255)
                                                                            gtcCGACTGACGGATAATGTGACTAATGGTATTTTTGTCAACATGCGCTAACACATTCTGACTGGTGGTTTCCCACCAGTCAGGCTGAATAAGATTACTCTGCTTTCTCCACAAAGATACCGTCCTGATGCCCTGCTTCATTAAAGAAc                                                            >  8:19684/1‑149 (MQ=255)
                                                                                                  cTAATGGTATTTTTGTCAACATGCGCTAACACATTCTGACTGGTGGTTTCCCACCAGTCAGGCTGAATAAGATTACTCTGCTTTCTCCACAAAGATACCGTCCTGATGCCCTGCTTCATTAAAGAACCAGATGCCGAGCGGGTAGTCtt                                      <  6:55152/149‑1 (MQ=255)
                                                                                                     aTGGTATTTTTGTCAACATGCGCTAACACATTCTGACTGGTGGTTTCCCACCAGTCAGGCTGAATAAGATTACTCTGCTTTCTCCACAAAGATACCGTCCTGATGCCCTGCTTCATTAAAGAACCAGAt                                                       <  7:11595/129‑1 (MQ=255)
                                                                                                     aTGGTATTTTTGTCAACATGCGCTAACACATTCTGACTGGTGGTTTCCCACCAGTCAGGCTGAATAAGATTACTCTGCTTTCTCCACAAAGATACCGTCCTGATGCCCTGCTTCATTAAAGAACCAGAt                                                       >  8:11595/1‑129 (MQ=255)
                                                                                                                                acaTTCTGACTGGTGGTTTCCCACCAGTCAGGCTGAATAAGATTACTCTGCTTTCTCCACAAAGATACCGTCCTGATGCCCTGCTTCATTAAAGAACCAGATGCCGAGCGGGTAGTCTTCCAGCGAAACCAgg                        >  3:84812/1‑133 (MQ=255)
                                                                                                                                acaTTCTGACTGGTGGTTTCCCACCAGTCAGGCTGAATAAGATTACTCTGCTTTCTCCACAAAGATACCGTCCTGATGCCCTGCTTCATTAAAGAACCAGATGCCGAGCGGGTAGTCTTCCAGCGAAACCAgg                        <  4:84812/133‑1 (MQ=255)
                                                                                                                                      tGACTGGTGGGAAACCACCAGTCAGGCTGAATAAGATTACTCTGCTTTCTCCACAAAGATACCGTCCTGATGCCCTGCTTCATTAAAGAACCAGATGCCGAGCGGGTAGTCTTCCAGCGAAACCAGGTACATTGTGCCTTCACTAAAct  >  1:86667/1‑149 (MQ=255)
                                                                                                                                      tGACTGGTGGGAAACCACCAGTCAGGCTGAATAAGATTACTATGCTTTCTCCACAAAGATACCGTCCTGATGCCCTGCTTCATTAAAGCACCAGATGCCGAGCGGGTAGTCTTCCAGCGAAACCAGGTACATTGTGCCTTCACTAAAct  >  4:90669/1‑149 (MQ=255)
                                                                                                                                            gtggtTTCCCACCAGTCAGGCTGAATAAGATTACTCTGCTTTCTCCACAAAGATACCGTCCTGATGCCCTGCTTCATTAAAGAACCAGATGCCGAGCGGGTAGTCTTCCAGCGAAACCAGGTACATtgt                <  5:74577/129‑1 (MQ=255)
                                                                                                                                             tggtTTCCCACCAGTCAGGCTGAATAAGATTACTCTGCTTTCTCCACAAAGATACCGTCCTGATGCCCTGCTTCATTAAAGAACCAGATGCCGAGCGGGTAGTCTTCCAGCGAAACCAGGTACAtt                  >  2:5381/1‑126 (MQ=255)
                                                                                                                                             tggtTTCCCACCAGTCAGGCTGAATAAGATTACTCTGCTTTCTCCACAAAGATACCGTCCTGATGCCCTGCTTCATTAAAGAACCAGATGCCGAGCGGGTAGTCTTCCAGCGAAACCAGGTACAtt                  <  1:5381/126‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                          
CAAGACCGTTTCATCGCATAAAGGTAATATTAAACGTAAGATCAAAACGCATAATAAACAGGTTATCTACCATGTCGTCCGACTGACGGATAATGTGACTAATGGTATTTTTGTCAACATGCGCTAACACATTCTGACTGGTGGTTTCCCACCAGTCAGGCTGAATAAGATTACTCTGCTTTCTCCACAAAGATACCGTCCTGATGCCCTGCTTCATTAAAGAACCAGATGCCGAGCGGGTAGTCTTCCAGCGAAACCAGGTACATTGTGCCTTCACTAAACT  >  NC_000913/2024465‑2024747

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