Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,013,589	G→A	11.7%	noncoding (346/385 nt)	PP_mr44 →	CrcZ ncRNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,013,589	0	G	A	11.7%	87.2 / 2.7	34	noncoding (346/385 nt)	PP_mr44	CrcZ ncRNA
Reads supporting (aligned to +/- strand):  ref base G (13/17);  new base A (2/2);  total (15/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.66e-01

AAAGCAGGTCAGCGCTCAAAAAAATACGTTTGCTCTTGACCCCGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCTCCACCAGCAACGGCATCGTGCCCAGTACCAGCAGGGCCAGTACCGTGCTGATCAGCGCCGTGGCTTCGCGCCCCATGC  >  NC_002947/4013443‑4013733                                                                                                                                                   |                                                                                                                                                 aaaGCTGGGGAGCGCTCAACATAATACGTTTGCTCTTGACCCCGCCTGGGGGTCGCCCGCAACACCTGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAAAACAAGAGCAGGCTACAAGGCTTTGAGGGGAGCTTCgg                                                                                                                                                 <  2:244697/148‑1 (MQ=255) aaaGCAGCTCAGCGCTCAAAAAAATACGTTTGCTCTTGAGCCCGCCTGGGGGTCGCCCGCATCACCACGACAAGCCGACAAAAACAACAATACGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCgg                                                                                                                                                 <  4:15852/148‑1 (MQ=255)         tCAGCGCTCAAAAAAATACGTTTGCTCTTGACCCCGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCt                                                                                                                                          >  8:172905/1‑147 (MQ=255)         tCAGCGCTCAAAAAAATACGTTTGCTCTTGACCCCGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCa                                                                                                                                        <  4:63830/149‑1 (MQ=255)           aGCTCTCAAAAAAATACGTTTGCTCTTGACCCCGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCCGAAGCTCCCCTCaaa                                                                                                                                      <  3:247846/149‑3 (MQ=255)                    aaaaaTACGTTTGCTCTTGACCCCGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGc                                                                                                                             <  7:167377/149‑1 (MQ=255)                       aaTACGTTTGCTCTTGACCCCGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCct                                                                                                                          <  3:55645/149‑1 (MQ=255)                         tACGTTTGCTCTTGACCCCGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCctct                                                                                                                        <  3:292640/149‑1 (MQ=255)                            gTTTGCTCTTGACCCCGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTtcc                                                                                                                     >  2:195047/1‑149 (MQ=255)                             tttGCTCTTGACCCCGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTtcct                                                                                                                    >  8:245507/1‑149 (MQ=255)                                gCTCTTGACCCCGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTtcctcc                                                                                                                  <  7:172905/148‑1 (MQ=255)                                 ctctTGACCCCGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAAAATAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCtt                                                                                                                <  4:66935/149‑1 (MQ=255)                                 ctctTGACCCCGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCtt                                                                                                                <  7:103654/149‑1 (MQ=255)                                  tctTGACCCCGCCTGGGGGTCGCTAGCAAAACCAGGACAAGCAGAAAAAAACAACAAAAGGCCGGTACAAATAATAAAAAAAGAGCAGGCAAAAACGATTTGAGGGGAGCTTCGGCTCCCATCAGTGCTTCCTGCCCTCTTCCTACTTc                                                                                                               >  7:48346/1‑149 (MQ=255)                                     tGACCCCGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGc                                                                                                            >  6:251423/1‑149 (MQ=255)                                         cccGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTtct                                                                                                        >  4:287191/1‑149 (MQ=255)                                                                        aaGCCGCCAAAAACAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCGTCCTCCTTCTGCTTCACCTCATCCTGCCCGTCGACCTTCTCCACCAg                                                                         >  6:140121/1‑149 (MQ=255)                                                                        aaGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCTCCACCAg                                                                         >  4:39565/1‑149 (MQ=255)                                                                         aGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAAACAAAGAGCAGGCTACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCACCACCAGc                                                                        <  6:307792/149‑1 (MQ=255)                                                                             gACAAAAACAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCTCCACCAGCAACg                                                                    <  7:245507/149‑1 (MQ=255)                                                                              aCAAAAACAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCTCCACCAGCAACgg                                                                   >  2:43744/1‑149 (MQ=255)                                                                                      aacaaTAGGCCGGTCCCAATAATAAAAAAAGAGCAGGAAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTACCATCTTCATCCTTCTGCTTATCCTCATCCTGCCCAACGACCATATCCACAAGCAACGACATCGAGc                                                           >  2:293584/1‑149 (MQ=255)                                                                                          aTAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCTCCACCAGCAACGGCATCGTGCCCAg                                                       <  2:14598/149‑1 (MQ=255)                                                                                                     ccAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCTCCACCAGCAACGGCATCGTGCCCa                                                        <  6:213516/137‑1 (MQ=255)                                                                                                       aataatAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCTCCACCAGCAACGGCATCGTGCCCAGTACCAGCAGGGc                                           >  5:166384/1‑148 (MQ=255)                                                                                                              aaaaaaGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCTCCACCAGCAACGGCATCGTGCCCAGTACCAGCAGGGCCAGTAcc                                    <  3:287191/148‑1 (MQ=255)                                                                                                                aaaaGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCTCCACCAGCAACGGCATCGTGCCCAGTACCAGCAGGGCCAGTACCGTg                                 <  1:43744/149‑1 (MQ=255)                                                                                                                  aaGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCTCCACCAGCAACGGCATCGTGCCCAGTACCAGCAGGGCCAGTACCGTGCt                               <  1:170518/149‑1 (MQ=255)                                                                                                                   atagCAGTCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCTCCACCAGCAACGGCATCGTGCCCAGTACCAGCAGGGCCAGTACCGTGCTg                              <  5:317692/147‑1 (MQ=255)                                                                                                                       cAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCTCCACCAGCAACGGCATCGTGCCCAGTACCAGCAGGGCCAGTACCGTGCTGAt                            >  7:175207/1‑147 (MQ=255)                                                                                                                            aacaacGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCTCCACCAGCAACGGCATCGTGCCCAGTACCAGCAGGGCCAGTACCGTGCTGATCAGCGcc                     >  3:279681/1‑149 (MQ=255)                                                                                                                                    actGAGGGGAGCCGAAGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCTCCACCAGCAACGGCATc                                                               <  2:114762/97‑1 (MQ=38)                                                                                                                                    actGAGGGGAGCCGAAGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCTCCACCAGCAACGGCATc                                                               >  1:114762/3‑99 (MQ=38)                                                                                                                                      tGAGGGGAGCCGAAGCTCCCCTCAGTGCTTCCTGCCCTCGTCCTCCTTCTGCTTCTCCTCCTCCTGCCCGTCGCCCTTCTCCACCAGCAACGCCATCGTGCCCAGTACCAGCAGGGCCAGTACCGTGCTGATCAGCGCCGTGGCTTcgc           >  5:328556/1‑149 (MQ=255)                                                                                                                                            gAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCTCCACCAGCAACGGCATCGTGCCCAGTACCAGCAGGGCCAGTACCGTGCTGATCAGCGCCGTGGCTTCGCGCCCCa     <  5:251423/149‑1 (MQ=255)                                                                                                                                               cTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCTCCACCAGCAACGGCATCGTGCCCAGTACCAGCAGGGCCAGTACCGTGCTGATCAGCGCCGTGGCTTCGCGCCCCATGc  <  6:166384/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                 AAAGCAGGTCAGCGCTCAAAAAAATACGTTTGCTCTTGACCCCGCCTGGGGGTCGCCCGCAACACCAGGACAAGCCGACAAAAACAACAATAGGCCGGTCCCAATAATAAAAAAAGAGCAGGCAACAACGCTTTGAGGGGAGCTTCGGCTCCCCTCAGTGCTTCCTGCCCTCTTCCTCCTTCTGCTTCTCCTCATCCTGCCCTTCGACCTTCTCCACCAGCAACGGCATCGTGCCCAGTACCAGCAGGGCCAGTACCGTGCTGATCAGCGCCGTGGCTTCGCGCCCCATGC  >  NC_002947/4013443‑4013733

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