Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,350,618	0	T	C	40.5%	16.7 / 47.9	42	noncoding (42/121 nt)	arcZ	small regulatory RNA ArcZ
Reads supporting (aligned to +/- strand):  ref base T (14/11);  new base C (8/9);  total (22/20)
Fisher's exact test for biased strand distribution p-value = 7.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TGACTGAAACACATTAATTTTTTTAATAAAAATGGTACGCATCACACATTTAACTGATTCATGTAACAAATCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCAT  >  NC_000913/3350480‑3350757                                                                                                                                           |                                                                                                                                            tGACTGAAACACATTAATTTTTTTAATAAAAATGGTACGCATCACACATTTAACTGATTCATGTAACAAATCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCata                                                                                                                                            <  1:260099/140‑1 (MQ=255)     tGAAACACATTAATTTTTTTAATAAAAATGGTACGCATCACACATTTAACTGATTCATGTAACAAATCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCAtaataa                                                                                                                                        >  1:225111/1‑140 (MQ=255)             aTTAATTTTTTTAATAAAAATGGTACGCATCACACATTTAACTGATTCATGTAACAAATCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGgc                                                                                                                                >  1:163089/1‑140 (MQ=255)                     ttttAATAAAAATGGTACGCATCACACATTTAACTGATTCATGTAACAAATCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCACAATAATTTACGGCGCAGCCaa                                                                                                                        <  1:12906/140‑1 (MQ=255)                             aaaaTGGTACGCATCACACATTTAACTGATTCATGTAACAAATCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTccc                                                                                                                <  2:414188/140‑1 (MQ=255)                                        cATCACACATTTAACTGATTCATGTAACAAATCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGgcg                                                                                                     <  1:780898/140‑1 (MQ=255)                                        cATCACACATTTAACTGATTCATGTAACAAATCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCACAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGgcg                                                                                                     <  2:45054/140‑1 (MQ=255)                                         aTCACACATTTAACTGATTCATGTAACAAATCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCACAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGgcgc                                                                                                    >  2:675395/1‑140 (MQ=255)                                                aTTTAACTGATTCATGTAACAAATCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTc                                                                                             >  1:279857/1‑140 (MQ=255)                                                  ttAACTGATTCATGTAACAAATCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTcgc                                                                                           >  1:420234/1‑140 (MQ=255)                                                   tAACTGATTCATGTAACAAATCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTcgcg                                                                                          <  2:420234/140‑1 (MQ=255)                                                   tAACTGATTCATGTAACAAATCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCACAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTcgcg                                                                                          <  2:596104/140‑1 (MQ=255)                                                     aCTGATTCATGTAACAAATCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCa                                                                                        >  2:51733/1‑140 (MQ=255)                                                     aCTGATTCATGTAACAAATCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCa                                                                                        >  1:801396/1‑140 (MQ=255)                                                        gATTCATGTAACAAATCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCACAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCAccc                                                                                     <  1:675395/140‑1 (MQ=255)                                                                   cAAATCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCACAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCg                                                                          >  1:547044/1‑140 (MQ=255)                                                                           ttAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCACAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCAtt                                                                  >  2:192352/1‑140 (MQ=255)                                                                                   tGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTg                                                          >  2:679947/1‑140 (MQ=255)                                                                                         cTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGCTTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTgg                                                         >  2:566458/1‑135 (MQ=255)                                                                                         cTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTgg                                                         <  1:566458/135‑1 (MQ=255)                                                                                          ttAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCACAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTg                                                   >  1:596474/1‑140 (MQ=255)                                                                                             aCTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCACAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGcca                                                >  2:309666/1‑140 (MQ=255)                                                                                                gcgtgcgGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGccaccc                                             >  2:692099/1‑140 (MQ=255)                                                                                                gcgtgcgGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGccaccc                                             <  1:233469/140‑1 (MQ=255)                                                                                                 cgtgcgGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCACAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGccaccca                                            >  1:290358/1‑140 (MQ=255)                                                                                                  gtgcgGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGccacccac                                           >  2:769220/1‑140 (MQ=255)                                                                                                  gtgcgGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGccacccac                                           >  2:207490/1‑140 (MQ=255)                                                                                                  gtgcgGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCACAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGccacccac                                           >  2:59932/1‑140 (MQ=255)                                                                                                   tgcgGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACg                                          >  1:481460/1‑140 (MQ=255)                                                                                                   tgcgGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACg                                          >  1:450464/1‑140 (MQ=255)                                                                                                   tgcgGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACg                                          >  2:761949/1‑140 (MQ=255)                                                                                                   tgcgGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCACAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACg                                          >  1:494285/1‑140 (MQ=255)                                                                                                       gCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCACAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCttt                                      <  1:192352/140‑1 (MQ=255)                                                                                                          tGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAg                                   <  1:51733/140‑1 (MQ=255)                                                                                                          tGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGCGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAg                                   <  2:801396/140‑1 (MQ=255)                                                                                                          tGAAAAACAGTGCTGTGCCCTTGTAACTCATCACAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAg                                   <  2:290358/140‑1 (MQ=255)                                                                                                                          gCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGa                   <  2:481460/140‑1 (MQ=255)                                                                                                                             cTTGTAACTCATCACAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGc                <  2:537183/140‑1 (MQ=255)                                                                                                                             cTTGTAACTCATCACAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGc                <  1:309666/140‑1 (MQ=255)                                                                                                                              ttGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGcc               <  2:163089/140‑1 (MQ=255)                                                                                                                               tGTAACTCATCACAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCa              <  1:59932/140‑1 (MQ=255)                                                                                                                                gTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCAc             >  1:129100/1‑140 (MQ=255)                                                                                                                                 tAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTTTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCAct            <  1:207490/140‑1 (MQ=255)                                                                                                                                 tAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCAct            <  1:769220/140‑1 (MQ=255)                                                                                                                                         catAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTAcgtcg                      >  1:625028/1‑122 (MQ=255)                                                                                                                                         catAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTAcgtcg                      <  2:625028/122‑1 (MQ=255)                                                                                                                                         catAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTc    >  1:3982/1‑140 (MQ=255)                                                                                                                                           taataaTTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCAt  <  1:662791/140‑1 (MQ=255)                                                                                                                                           taataaTTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCAt  >  1:330871/1‑140 (MQ=255)                                                                                                                                           taataaTTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCAt  <  1:203880/140‑1 (MQ=255)                                                                                                                                           |                                                                                                                                            TGACTGAAACACATTAATTTTTTTAATAAAAATGGTACGCATCACACATTTAACTGATTCATGTAACAAATCATTTAAGTTTTGCTATCTTAACTGCGTGCGGCCTGAAAAACAGTGCTGTGCCCTTGTAACTCATCATAATAATTTACGGCGCAGCCAAGATTTCCCTGGTGTTGGCGCAGTATTCGCGCACCCCGGTCTAGCCGGGGTCATTTTTTAGTGGCTTTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCAT  >  NC_000913/3350480‑3350757

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