Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP009273	4,596,171	G→A	noncoding (48/87 nt)	leuQ ←	tRNA‑Leu

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,596,171	0	G	A	91.8%	115.3 / 1.8	49	noncoding (48/87 nt)	leuQ	tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base A (21/24);  total (23/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 = 4.54e-01

CCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTG  >  CP009273/4596043‑4596315                                                                                                                                 |                                                                                                                                                 ccAAGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAg                                                                                                                              <  1:732781/149‑1 (MQ=255)                                                                             atatcgcctcgtttGGTGCGAGGGGGGGGACTTGAACCCCCATGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGt                                                  <  1:55481/138‑1 (MQ=18)                                                                             atatcgcctcgtttGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGt                                                  >  1:1827545/12‑149 (MQ=17)                                                                              taccgcctcgtttGGGGCGAGGGGGGGGACCTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTg                                                 <  2:99904/139‑1 (MQ=18)                                                                                tcgcctcgtttGGTGCGAGGGGGGTGACTGGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGt                                               <  1:953408/141‑1 (MQ=18)                                                                                tcgcctcgtttGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGt                                               >  1:1703500/9‑149 (MQ=255)                                                                                  gcgtcctttGGTGCGCGGGGGGGGACTTGAAACCCCCCGTTCGTCCGAAAACTAACACCTGAAGCTAGCGCGTCTACCAATCCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGttt                                             <  1:422594/143‑1 (MQ=21)                                                                                   cctcgtttGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCCAGGCGGCGCAGTCTAGTGGtttt                                            >  2:111471/6‑149 (MQ=18)                                                                                   cctcgtttGGTGCGACGGGGGGTACTTCAACCCACACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGtttt                                            <  1:1226838/144‑1 (MQ=21)                                                                                    ctcgtttGGGGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAACAACACGAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTc                                           <  1:1762871/145‑1 (MQ=21)                                                                                     tcgtttGGTGCGAGGTGGGGGCCTTGAACCCCCACGTTCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCa                                          <  2:1837259/146‑1 (MQ=21)                                                                                     tcgtttGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCCACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCACTCTAGTGGGATTCa                                          >  2:660994/4‑149 (MQ=18)                                                                                     gTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTAACACGGAGGCGCATTCTAGTGGTTTTCa                                          >  2:1496255/1‑149 (MQ=255)                                                                                       gattGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGc                                        <  2:1680256/147‑1 (MQ=255)                                                                                          tGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCttt                                     <  2:1640469/149‑1 (MQ=255)                                                                                            gTGCGCGCGGGGGAACTTGAACCCCCACGTCCGTAAGAACACTACCACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTc                                   <  2:1190782/149‑1 (MQ=25)                                                                                            gTGCGAGGGGGGGGACGTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTc                                   <  1:514239/149‑1 (MQ=255)                                                                                             tGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCg                                  >  1:1456272/1‑149 (MQ=255)                                                                                             tGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCg                                  >  1:735124/1‑149 (MQ=255)                                                                                             tGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCg                                  >  2:1570618/1‑149 (MQ=255)                                                                                              gCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGGAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGt                                 >  1:394958/1‑149 (MQ=255)                                                                                              gCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCTACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGt                                 <  2:1598089/149‑1 (MQ=255)                                                                                              gCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGt                                 <  1:1679348/149‑1 (MQ=255)                                                                                              gCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGt                                 <  2:579286/149‑1 (MQ=255)                                                                                                   gggggggACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATa                            <  1:125718/149‑1 (MQ=255)                                                                                                   ggggcggACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATa                            <  2:1348945/149‑1 (MQ=255)                                                                                                   ggggccgACTGGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCACCTTTTCGTCAATa                            >  2:17582/7‑149 (MQ=25)                                                                                                     gggggACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTTAAATTATTTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTATTCGTCTATGGt                          >  2:498850/1‑149 (MQ=21)                                                                                                       gggACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTTGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTa                        <  1:1365393/149‑1 (MQ=255)                                                                                                       gggACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCCCCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCACTAGTTa                        >  1:863674/1‑149 (MQ=255)                                                                                                       gggACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTa                        >  1:212796/1‑149 (MQ=255)                                                                                                         gtcTTGAACCCCCACGTCCGTAAGAACACTGACTCCTGTCGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATGCCATCAATTCTTATTTAGTATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCTATAGTTAAt                      <  2:254451/147‑1 (MQ=21)                                                                                                            ttGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTAt                   >  1:1543151/1‑149 (MQ=255)                                                                                                            ttGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCCAGTGGTTTTCAGCTTTTCGTCAATAGTTAAtaac                   >  2:1654991/1‑146 (MQ=255)                                                                                                              gAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCg                 >  1:242234/1‑149 (MQ=255)                                                                                                               aaCCCCCCCGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTCAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGa                <  1:1215458/149‑1 (MQ=255)                                                                                                               aaCCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGa                >  1:154901/1‑149 (MQ=255)                                                                                                               aaCCCCCACGTCCGTAAAAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGa                <  1:718848/149‑1 (MQ=255)                                                                                                                aCCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGATTTCAGCTTTTCGTCAATAGTTAATTATCGAc               >  1:1108423/1‑149 (MQ=255)                                                                                                                    ccACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACGATCAATTCTTAAAAACAATTGCTACCACGGAGGCGCATCCTAGTGGCTTTCAGCTTTTCGTCAATAGTTAATTATCGAcgcag           >  2:601259/1‑145 (MQ=37)                                                                                                                     cACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAgg          <  1:543151/149‑1 (MQ=255)                                                                                                                     cACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAgg          <  2:987471/149‑1 (MQ=255)                                                                                                                      aCGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGgt         <  1:1935063/149‑1 (MQ=255)                                                                                                                         tCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTa      <  1:1689070/149‑1 (MQ=255)                                                                                                                          ccGTAAGAACACTACCACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCTATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTaa     <  2:1341736/149‑1 (MQ=255)                                                                                                                           cGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAAt    >  1:1842396/1‑149 (MQ=255)                                                                                                                            gTAAGGACACTAACAACTTAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGCTTTCAGCTTTTCGTCAATTGTTAATTATCGACCGAGGTGTAAtc   >  1:1476027/1‑148 (MQ=38)                                                                                                                            gTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCCCGGAGGCGCATCCTAGTGGTTCTCAGCTTTTCGCCAATAGTTAAAAACCGACCGAGGAGTAAtt   >  1:261643/1‑149 (MQ=25)                                                                                                                             tAAGAACACTAACACCTGAAGCTAGCGAGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTg  >  2:1797453/1‑149 (MQ=255)                                                                                                                                 |                                                                                                                                                 CCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTG  >  CP009273/4596043‑4596315

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