Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,390,401	C→A	5.2%	intergenic (‑41/+56)	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,401	0	C	A	5.2%	271.5 / 2.7	76	intergenic (‑41/+56)	psd/rsgA	phosphatidylserine decarboxylase/ribosome small subunit‑dependent GTPase A
Reads supporting (aligned to +/- strand):  ref base C (33/39);  new base A (2/2);  total (35/41)
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.99e-01

AGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCCGCGAT  >  NC_000913/4390254‑4390549                                                                                                                                                    |                                                                                                                                                      aGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCt                                                                                                                                                      <  3:122029/149‑1 (MQ=255)        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                                                                                                                               >  4:262973/1‑149 (MQ=255)                       gcttgcttGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑aaaaaaa                                                                                                                               <  2:310916/149‑1 (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                                                                                                                    <  2:283585/149‑1 (MQ=255)                                  ccccAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTc                                                                                                                    >  1:73027/1‑149 (MQ=255)                                   cccAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGAtt                                                                                                                     >  3:292941/1‑147 (MQ=255)                                   cccAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGAtt                                                                                                                     <  4:292941/147‑1 (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)                                                     tAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTAc                                                                                                                                                    >  1:153152/1‑99 (MQ=255)                                                     tAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTAc                                                                                                                                                    <  2:153152/99‑1 (MQ=255)                                                      aaGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAAtt                                                                                                <  1:169712/149‑1 (MQ=255)                                                      aaGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAAtt                                                                                                <  3:411621/149‑1 (MQ=255)                                                       aGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTg                                                                                               >  2:224539/1‑149 (MQ=255)                                                       aGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTg                                                                                               >  2:359823/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:359823/149‑1 (MQ=255)                                                                                gTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGttt                                                                      <  1:484213/149‑1 (MQ=255)                                                                                gTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGttt                                                                      <  1:224539/149‑1 (MQ=255)                                                                                 tAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGttt                                                                      >  2:489500/1‑148 (MQ=255)                                                                                 tAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGttt                                                                      <  1:489500/148‑1 (MQ=255)                                                                                    cGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGttt                                                                      <  1:285293/145‑1 (MQ=255)                                                                                    cGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGttt                                                                      >  2:285293/1‑145 (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                                  <  2:418307/148‑1 (MQ=255)                                                                                                                     cAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATa                                  >  1:418307/1‑148 (MQ=255)                                                                                                                        gCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGt                              <  4:466591/149‑1 (MQ=255)                                                                                                                                      gTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGaaa                >  3:72054/1‑149 (MQ=255)                                                                                                                                      gTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGaaa                <  3:271617/149‑1 (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)                                                                                                                                                    |                                                                                                                                                      AGTTTTGTCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCCGCGAT  >  NC_000913/4390254‑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: ‑