Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,555,796	(T)5→4	coding (525/831 nt)	yjiC ←	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,555,792	0	T	.	93.1%	227.4 / 11.1	58	coding (529/831 nt)	yjiC	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base . (30/24);  total (32/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

TGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAA  >  NC_000913/4555697‑4555887                                                                                                |                                                                                                tGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTTGTCCATCttttg                                                                                             <  1:382263/100‑2 (MQ=255)   cGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATCttttgcc                                                                                           <  2:1199608/100‑4 (MQ=255)   cGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATCttttgcc                                                                                           >  1:1170620/1‑97 (MQ=255)    ggCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATCttttgcct                                                                                          >  2:42097/1‑96 (MQ=255)      cAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATC‑TTTTGCCTga                                                                                       >  1:1204586/1‑100 (MQ=255)               aTCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATC‑TTTTGCCTGAGCAGCTGGa                                                                              <  1:685605/100‑1 (MQ=255)                tCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGaa                                                                             <  1:1045786/100‑1 (MQ=255)                 ccGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAg                                                                            <  1:646663/100‑1 (MQ=255)                 ccGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAg                                                                            >  2:657100/1‑100 (MQ=255)                  cGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAgg                                                                           >  1:271449/1‑100 (MQ=255)                  cGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAgg                                                                           >  1:499037/1‑100 (MQ=255)                   ggAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAggg                                                                          >  1:1262869/1‑100 (MQ=255)                   ggAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAggg                                                                          <  2:597755/100‑1 (MQ=255)                    gaagaaTTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAGGGt                                                                         <  2:407516/100‑1 (MQ=255)                      agaaTTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAGGGTTg                                                                       <  2:553118/100‑1 (MQ=255)                      agaaTTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAGGGTTg                                                                       >  1:522746/1‑100 (MQ=255)                        aaTTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAGGGTTGaa                                                                     >  2:412351/1‑100 (MQ=255)                        aaTTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAGGGTTGaa                                                                     <  2:138533/100‑1 (MQ=255)                         aTTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAGGGTTGaa                                                                     >  2:800203/1‑99 (MQ=255)                          tttGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAGGGTTGAAat                                                                   <  2:522742/100‑1 (MQ=255)                           ttGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAGGGTTGAAata                                                                  <  2:1059847/100‑1 (MQ=255)                           ttGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCATGGCTGCTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAGGGTTGAAata                                                                  <  2:1176489/100‑1 (MQ=255)                              cTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAGGGTTGAAatatat                                                               <  2:509214/100‑1 (MQ=255)                                aCCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATcc                                                             >  1:275673/1‑100 (MQ=255)                                    cTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATc                                                         >  2:348289/1‑100 (MQ=255)                                     tGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCg                                                        <  1:412355/100‑1 (MQ=255)                                       ttACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGt                                                       >  1:576193/1‑100 (MQ=255)                                        tACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTgg                                                     >  2:946733/1‑100 (MQ=255)                                            gTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTcc                                                 <  1:305653/100‑1 (MQ=255)                                               aTTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGc                                              <  1:194513/100‑1 (MQ=255)                                                ttATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCa                                             >  1:625249/1‑100 (MQ=255)                                                  aTCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATCTTTTTGCCTGAGAAGCTGGCAGGGTTGAAATATATCCTATCGTGGTTCCGGCAt                                            >  2:640028/1‑100 (MQ=255)                                                    cTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCGATC‑TTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAAc                                         <  2:675613/100‑1 (MQ=255)                                                      gATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAAc                                         >  2:139869/1‑98 (MQ=255)                                                      gATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCa                                       >  1:1048068/1‑100 (MQ=255)                                                        tgagtataGATATCAGTCCTTGGCGGGTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAgg                                     <  2:1032762/96‑1 (MQ=255)                                                        tCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAgg                                     <  1:42098/100‑1 (MQ=255)                                                          aaTATAGATATCAATCCTTGGCTGCTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGaa                                   <  2:1048058/100‑1 (MQ=255)                                                              taGATATCAATCCTTGGCTGCTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCTGGCATAACCAGGAACAtt                               >  2:1206976/1‑100 (MQ=255)                                                              taGATATCAATCCTTGGCTGCTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACAtt                               <  2:161807/100‑1 (MQ=255)                                                              taGATATCAATCCTTGGCTGCTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGCCATAACCAGGAACAtt                               >  2:761844/1‑100 (MQ=255)                                                               aGATATCAATCCTTGGCTGCTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTg                              <  1:1117942/100‑1 (MQ=255)                                                               aGATATCAATCCTTGGCTGCTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTg                              <  1:1220883/100‑1 (MQ=255)                                                                gATATCAATCCTTGGCTGCTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGc                             >  1:512240/1‑100 (MQ=255)                                                                 atatCAATCCTTGGCTGCTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCa                            <  1:348292/100‑1 (MQ=255)                                                                 atatCAATCCTTGGCTGCTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCa                            >  2:277530/1‑100 (MQ=255)                                                                   atCAATCCTTGGCTGCTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAAc                          >  2:933506/1‑100 (MQ=255)                                                                      aaTCCTTGGCTGCTCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCa                        <  2:257165/100‑1 (MQ=255)                                                                      aaTCCTTGGCTGCTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAg                       >  1:410844/1‑100 (MQ=255)                                                                      aaTCCTTGGCTGCTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAg                       >  2:309075/1‑100 (MQ=255)                                                                      aaTCCTTGGCTGCTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAg                       >  1:926580/1‑100 (MQ=255)                                                                        tCCTTGGCTGCTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAggg                     >  1:500502/1‑100 (MQ=255)                                                                        tCCTTGGCTGCTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAggg                     >  2:653684/1‑100 (MQ=255)                                                                        tCCTTGGCTGCTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAggg                     >  1:8542/1‑100 (MQ=255)                                                                             ggCTGCTCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTaa                >  1:1063840/1‑100 (MQ=255)                                                                               ctgctCATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATg              <  2:1136495/100‑1 (MQ=255)                                                                                   tCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTa           <  2:576189/100‑1 (MQ=255)                                                                                    cATTGGTCCATC‑TTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTAtt         >  1:228179/1‑100 (MQ=255)                                                                                       tGGTCCATC‑TTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCa      <  1:933514/100‑1 (MQ=255)                                                                                        ggTCCATC‑TTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCAt     >  2:370488/1‑100 (MQ=255)                                                                                           ccATC‑TTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGaa  >  2:571259/1‑100 (MQ=255)                                                                                           ccATC‑TTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCAGAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGaa  >  2:1155125/1‑100 (MQ=255)                                                                                                |                                                                                                TGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAA  >  NC_000913/4555697‑4555887

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