Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,390,406	T→G	5.4%	intergenic (‑46/+51)	psd ← / ← rsgA	phosphatidylserine decarboxylase/ribosome small subunit‑dependent GTPase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,390,406	0	T	G	5.4%	261.6 / 6.5	74	intergenic (‑46/+51)	psd/rsgA	phosphatidylserine decarboxylase/ribosome small subunit‑dependent GTPase A
Reads supporting (aligned to +/- strand):  ref base T (33/37);  new base G (2/2);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.39e-01

TCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCCGCGAT  >  NC_000913/4390261‑4390549                                                                                                                                                  |                                                                                                                                                 tCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGcc                                                                                                                                               <  3:233469/149‑1 (MQ=255) tCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGcc                                                                                                                                               >  3:369720/1‑149 (MQ=255)   aGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCgg                                                                                                                                             >  4:98895/1‑149 (MQ=255)   aGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCgg                                                                                                                                             >  3:48740/1‑149 (MQ=255)     ccATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGAt                                                                                                                                           >  1:49476/1‑149 (MQ=255)      cATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATc                                                                                                                                          <  4:48740/149‑1 (MQ=255)      cATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGCAAGCCATAGTTTCGGCAGAACGTACTGTAGCGAAAGTTTAAGTGAATTTCACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCCGCCTACATGCCGGATc                                                                                                                                          <  3:430419/149‑1 (MQ=255)           tGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCATGTAGGCCGGATc                                                                                                                                          >  1:403337/1‑144 (MQ=255)           tGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCATGTAGGCCGGATc                                                                                                                                          <  2:403337/144‑1 (MQ=255)              ccGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑aaaaa                                                                                                                                 >  2:169712/1‑149 (MQ=255)                gcttgcttGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑aaaaaaa                                                                                                                               <  2:310916/149‑1 (MQ=255)                gcttgcttGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑aaaaaaa                                                                                                                               >  4:262973/1‑149 (MQ=255)                      ttgcGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGAc                                                                                                                         >  2:190116/1‑149 (MQ=255)                      ttgcGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGAc                                                                                                                         >  2:362107/1‑149 (MQ=255)                        gcgcCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCATAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGGCTCCAGGCCATTGTTTTGTCGTTCCTGGTCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGa                                                                                                                       <  3:580781/149‑1 (MQ=255)                           ccccAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTc                                                                                                                    >  1:73027/1‑149 (MQ=255)                           ccccAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTc                                                                                                                    <  2:283585/149‑1 (MQ=255)                            cccAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGAtt                                                                                                                     <  4:292941/147‑1 (MQ=255)                            cccAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGAtt                                                                                                                     >  3:292941/1‑147 (MQ=255)                             ccAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCGGATCCGGCATGTAGGCCGGATcaggaa                                                                                                                                     <  2:377487/131‑6 (MQ=255)                             ccAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACCGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCGGATCCGGCATGTAGGCCGGATcaggaa                                                                                                                                     >  1:377487/1‑126 (MQ=255)                                     ccAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGCTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGt                                                                                                          <  3:674363/149‑1 (MQ=255)                                       aGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGATTTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGCCGATTCTAACGACGGTTa                                                                                                        <  4:164536/149‑1 (MQ=255)                                             gTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTATCCTCCAGGCCATTGTTTTGTCGTTGCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAac                                                                                                                >  3:36301/1‑135 (MQ=255)                                             gTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAac                                                                                                                <  4:36301/135‑1 (MQ=255)                                               aaGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAAtt                                                                                                <  1:169712/149‑1 (MQ=255)                                               aaGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAAtt                                                                                                <  3:411621/149‑1 (MQ=255)                                                aGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTg                                                                                               >  2:359823/1‑149 (MQ=255)                                                aGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTg                                                                                               >  2:224539/1‑149 (MQ=255)                                                   cATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTgtca                                                                                            <  3:262973/149‑1 (MQ=255)                                                    aTAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTgtcag                                                                                           <  3:26070/149‑1 (MQ=255)                                                         ttCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATc                                                                                      <  2:180595/149‑1 (MQ=255)                                                         ttCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATc                                                                                      <  3:218200/149‑1 (MQ=255)                                                          tCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATcc                                                                                     <  4:369720/149‑1 (MQ=255)                                                              cAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTAt                                                                                 >  2:484213/1‑149 (MQ=255)                                                              cAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTAt                                                                                 <  2:73027/149‑1 (MQ=255)                                                                   tGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGaa                                                                            >  3:316207/1‑149 (MQ=255)                                                                        tGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGtt                                                                       >  1:657306/1‑149 (MQ=255)                                                                         gTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGttt                                                                      <  1:224539/149‑1 (MQ=255)                                                                         gTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGttt                                                                      <  1:359823/149‑1 (MQ=255)                                                                         gTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGttt                                                                      <  1:484213/149‑1 (MQ=255)                                                                          tAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGttt                                                                      <  1:489500/148‑1 (MQ=255)                                                                          tAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGttt                                                                      >  2:489500/1‑148 (MQ=255)                                                                             cGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGttt                                                                      >  2:285293/1‑145 (MQ=255)                                                                             cGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGttt                                                                      <  1:285293/145‑1 (MQ=255)                                                                             cGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGtt                                                                       <  2:558379/144‑1 (MQ=255)                                                                             cGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGtt                                                                       >  1:558379/1‑144 (MQ=255)                                                                               aaaGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTAcgc                                                                >  2:645963/1‑149 (MQ=255)                                                                                   tttAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAgg                                                                                                                             <  1:147325/84‑1 (MQ=255)                                                                                   tttAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAgg                                                                                                                             >  2:147325/1‑84 (MQ=255)                                                                                    ttAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGtttt                                                           <  1:5560/149‑1 (MQ=255)                                                                                       aaTGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGaaaa                                                                          <  3:449856/131‑1 (MQ=255)                                                                                       aaTGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAATCATCCGTATCAGaaaa                                                                          >  4:449856/1‑131 (MQ=255)                                                                                        aTGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTAcc                                                       >  4:271617/1‑149 (MQ=255)                                                                                           aaTTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTgc                                                    <  2:321154/149‑1 (MQ=255)                                                                                             tttAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTgcgc                                                  >  1:68001/1‑149 (MQ=255)                                                                                                aaCAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCAt                                               >  2:509231/1‑149 (MQ=255)                                                                                                       tAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTcc                                        <  3:98895/149‑1 (MQ=255)                                                                                                         gCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAg                                      >  3:112848/1‑149 (MQ=255)                                                                                                              cAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATa                                  >  1:418307/1‑148 (MQ=255)                                                                                                              cAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATa                                  <  2:418307/148‑1 (MQ=255)                                                                                                                 gCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGt                              <  4:466591/149‑1 (MQ=255)                                                                                                                               gTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGaaa                <  3:271617/149‑1 (MQ=255)                                                                                                                               gTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGaaa                >  3:72054/1‑149 (MQ=255)                                                                                                                                ttCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATCGTTTTCGAAAc               <  1:362107/149‑1 (MQ=255)                                                                                                                                  ccTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACgg             <  4:112848/149‑1 (MQ=255)                                                                                                                                  ccTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACgg             <  1:509231/149‑1 (MQ=255)                                                                                                                                    tGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGt           >  3:628070/1‑149 (MQ=255)                                                                                                                                       tCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTc        >  1:509366/1‑149 (MQ=255)                                                                                                                                        ccGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTcc       <  4:628070/149‑1 (MQ=255)                                                                                                                                          ggCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCcgc     <  1:267070/149‑1 (MQ=255)                                                                                                                                           gCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCcgcg    >  4:674688/1‑149 (MQ=255)                                                                                                                                           gCCTACATGCCGGATCCCGAAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCcgc     >  3:343893/1‑149 (MQ=255)                                                                                                                                             cTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCCGCGAt  >  2:503151/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                 TCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCCGCGAT  >  NC_000913/4390261‑4390549

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