Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,285,866	Δ1 bp	6.7%	coding (1506/1959 nt)	acs ←	acetyl‑CoA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,285,866	0	C	.	6.7%	312.8 / 15.9	90	coding (1506/1959 nt)	acs	acetyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base C (46/38);  new base . (3/3);  total (49/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 = 7.00e-01

AATACCTACTACGGCGGCTTCGGCAATCTTCGGATGCGCCACCAGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTA  >  NC_000913/4285722‑4286013                                                                                                                                                 |                                                                                                                                                     aaTACCTACTACGGCGGCTTCGGCAATCTTCGGATGCGCCACCAGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGAcgcgcg                                                                                                                                                   >  2:498387/1‑148 (MQ=255)        actacGGCGGCTTCGGCAATCTTCGGATGCGCCACCAGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACG‑GCGCGCGTcgc                                                                                                                                           >  1:702005/1‑148 (MQ=255)          tacGGCGGCTTCGGCAATCTTCGGATGCGCCACCAGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACG‑GCGCGCGTcgcga                                                                                                                                         <  2:702005/148‑3 (MQ=255)                ggcTTCGGCAATCTTCGGATGCGCCACCAGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGc‑c                                                                                                                                                >  4:395696/1‑135 (MQ=255)                ggcTTCGGCAATCTTCGGATGCGCCACCAGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGc‑c                                                                                                                                                <  3:395696/135‑1 (MQ=255)                   ttCGGCAATCTTCGGATGCGCCACCAGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATAc                                                                                                                                <  1:445768/148‑1 (MQ=255)                   ttCGGCAATCTTCGGATGCGCCACCAGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACa                                                                                                                               <  3:409151/149‑1 (MQ=255)                     cGGCAATCTTCGGATGCGCCACCAGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACata                                                                                                                             >  4:244268/1‑149 (MQ=255)                       gCAATCTTCGGATGCGCCACCAGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATAtt                                                                                                                           >  4:270425/1‑149 (MQ=255)                             ttCGGATGCGCCACCAGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGaa                                                                                                                     <  3:809797/149‑1 (MQ=255)                             ttCGGATGCGCCACCAGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGaa                                                                                                                     <  4:230173/149‑1 (MQ=255)                             ttCGGATGCGCCACCAGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGaa                                                                                                                     <  4:506637/149‑1 (MQ=255)                              tCGGATGCGCCACCAGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGCAATACATATTTTTGaa                                                                                                                     >  1:588217/1‑148 (MQ=255)                              tCGGATGCGCCACCAGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGaa                                                                                                                     >  1:562915/1‑148 (MQ=255)                              tCGGATGCGCCACCAGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAg                                                                                                                    <  3:676786/149‑1 (MQ=255)                                 gATGCGCCACCAGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGt                                                                                                                  <  2:389893/148‑1 (MQ=255)                                 gATGCGCCACCAGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGt                                                                                                                  <  1:759910/148‑1 (MQ=255)                                 gATGCGCCACCAGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTg                                                                                                                 <  3:250243/149‑1 (MQ=255)                                 gATGCGCCACCAGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACCCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACG‑GCGCGCGTcgcgatg                                                                                                                                       >  4:590012/1‑123 (MQ=255)                                 gATCCGCCACCAGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACG‑GCGCGCGTcgcgatg                                                                                                                                       <  3:590012/127‑5 (MQ=255)                                        caccaGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGt                                                                                                          >  4:696745/1‑149 (MQ=255)                                        caccaGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTGATCCAGTAATAGCCATCTTCATCGCGACGCACGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGt                                                                                                          >  4:728322/1‑149 (MQ=255)                                          ccaGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTa                                                                                                         <  1:39031/148‑1 (MQ=255)                                          ccaGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAg                                                                                                        <  4:412059/149‑1 (MQ=255)                                          ccaGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAg                                                                                                        >  4:655682/1‑149 (MQ=255)                                            aGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGt                                                                                                          >  1:633780/1‑145 (MQ=255)                                            aGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGt                                                                                                          <  2:633780/145‑1 (MQ=255)                                            aGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAgg                                                                                                       <  1:498387/148‑1 (MQ=255)                                            aGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAgg                                                                                                       <  2:588217/148‑1 (MQ=255)                                                    cTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCa                                                                                               >  2:763284/1‑148 (MQ=255)                                                     tCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGt                                                                                                  <  3:589113/144‑1 (MQ=255)                                                     tCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGt                                                                                                  >  4:589113/1‑144 (MQ=255)                                                     tCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCaaa                                                                                             >  3:536564/1‑149 (MQ=255)                                                      cAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCaaa                                                                                             >  2:494878/1‑148 (MQ=255)                                                      cAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCaaa                                                                                             <  2:562915/148‑1 (MQ=255)                                                       aaTCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCaaaa                                                                                            <  1:763284/148‑1 (MQ=255)                                                        aTCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACG‑GCGCGCGTcgcgat                                                                                                                                        <  4:256319/103‑4 (MQ=255)                                                        aTCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACG‑GCGCGCGTcgcgat                                                                                                                                        >  3:256319/1‑100 (MQ=255)                                                           tctGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCTGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGtt                                                                                        <  1:405811/148‑1 (MQ=255)                                                           tctGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTc                                                                                       <  4:536564/149‑1 (MQ=255)                                                           tctGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTc                                                                                       <  3:240403/149‑1 (MQ=255)                                                            ctGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGGGAAGTAGGTCTGTTCAAAACGTTc                                                                                       >  1:89086/1‑148 (MQ=255)                                                            ctGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTc                                                                                       >  1:247426/1‑148 (MQ=255)                                                            ctGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTc                                                                                       >  2:728403/1‑148 (MQ=255)                                                            ctGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCg                                                                                      <  3:69171/149‑1 (MQ=255)                                                                 gTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGAt                                                                                  <  2:642351/148‑1 (MQ=255)                                                                 gTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATc                                                                                 >  3:225101/1‑149 (MQ=255)                                                                    cccAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGAt                                                                                  <  2:603487/145‑1 (MQ=255)                                                                    cccAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGAt                                                                                  >  1:603487/1‑145 (MQ=255)                                                                    cccAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGcc                                                                              >  4:443275/1‑149 (MQ=255)                                                                     ccAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGcc                                                                              >  2:532020/1‑148 (MQ=255)                                                                       aGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCaa                                                                            >  2:116142/1‑148 (MQ=255)                                                                            gTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGc                                                                       >  2:711603/1‑148 (MQ=255)                                                                             tGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCg                                                                      <  1:711603/148‑1 (MQ=255)                                                                              gACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGt                                                                     >  2:20148/1‑148 (MQ=255)                                                                                   gagaCGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTAcgcg                                                                >  2:765117/1‑148 (MQ=255)                                                                                     gaCGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTAcgc                                                                 <  4:807362/145‑1 (MQ=255)                                                                                     gaCGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTAcgc                                                                 >  3:807362/1‑145 (MQ=255)                                                                                     gaCGTTCAGCACGTCGTCCACACGCCCGGTTACCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGcc                                                              >  2:542635/1‑148 (MQ=255)                                                                                        gTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATc                                                                                 <  3:282419/126‑1 (MQ=255)                                                                                        gTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATc                                                                                 >  4:282419/1‑126 (MQ=255)                                                                                           cAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCTTGAccc                                                        >  1:644696/1‑148 (MQ=255)                                                                                                  tcgtcCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCTTGACCCGGCCAgg                                                 <  2:644696/148‑1 (MQ=255)                                                                                                  tcgtcCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAgg                                                 >  1:546525/1‑148 (MQ=255)                                                                                                  tcgtcCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAgg                                                 >  1:401572/1‑148 (MQ=255)                                                                                                    gtcCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCCGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGt                                              >  3:854189/1‑149 (MQ=255)                                                                                                     tcCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGt                                              <  2:247426/148‑1 (MQ=255)                                                                                                      ccACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCg                                            <  4:722577/149‑1 (MQ=255)                                                                                                        acacGCCCGGTTATCCAGTAATAGCCGTCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATAGTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAGCGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCACGAGTCGGt                                          >  3:560811/1‑149 (MQ=255)                                                                                                        acacGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGGGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCgg                                           <  2:89086/148‑1 (MQ=255)                                                                                                        acacGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGt                                          >  3:833598/1‑149 (MQ=255)                                                                                                          acGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGa                                        <  4:435000/149‑1 (MQ=255)                                                                                                             cccGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGAtt                                      >  1:251424/1‑148 (MQ=255)                                                                                                                 gTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCa                                  <  2:25665/148‑1 (MQ=255)                                                                                                                  ttATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAg                                 >  2:155775/1‑148 (MQ=255)                                                                                                                   tATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGc                               <  4:280276/149‑1 (MQ=255)                                                                                                                     tCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTa                             >  4:27253/1‑149 (MQ=255)                                                                                                                       cAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTAcc                           >  4:611765/1‑149 (MQ=255)                                                                                                                        aGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTAcc                           >  2:632741/1‑148 (MQ=255)                                                                                                                          tAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTc                        >  4:885882/1‑149 (MQ=255)                                                                                                                                ccATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTGTGTTCAAAACGTTCGTGATAGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGc                  >  3:556342/1‑149 (MQ=255)                                                                                                                                ccATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGc                  >  4:681573/1‑149 (MQ=255)                                                                                                                                 cATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGc                  >  2:856223/1‑148 (MQ=255)                                                                                                                                 cATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGc                  >  1:110915/1‑148 (MQ=255)                                                                                                                                  aTCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGccc                <  3:443275/149‑1 (MQ=255)                                                                                                                                  aTCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGccc                >  4:667896/1‑149 (MQ=255)                                                                                                                                      tCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCt             <  1:116142/148‑1 (MQ=255)                                                                                                                                       cATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTcc           <  3:270425/149‑1 (MQ=255)                                                                                                                                        aTCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTcc           <  1:155775/148‑1 (MQ=255)                                                                                                                                        aTCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCa          >  3:4886/1‑149 (MQ=255)                                                                                                                                              aCGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCggg     >  1:627930/1‑148 (MQ=255)                                                                                                                                                gcgcgc‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTa  <  3:27253/149‑1 (MQ=255)                                                                                                                                                 cgcgc‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTa  >  1:156500/1‑148 (MQ=255)                                                                                                                                                 |                                                                                                                                                     AATACCTACTACGGCGGCTTCGGCAATCTTCGGATGCGCCACCAGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGC‑CGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTA  >  NC_000913/4285722‑4286013

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