Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	388,743	G→A	intergenic (+97/+10)	tauD → / ← hemB	taurine dioxygenase/porphobilinogen synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	388,743	0	G	A	94.9%	126.5 / NA	39	intergenic (+97/+10)	tauD/hemB	taurine dioxygenase/porphobilinogen synthase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (18/19);  minor base C (2/0);  total (20/19)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACCGCGTGACCCAGCACTATGCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGATGTCGAGGTACGCTCCAGCAGGTTT  >  NC_000913/388497‑388989                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       cAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACCGCGTGACCCAGCACTATGCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAgaagaa                                                                                                                                                                                                                                                     >  2:204096/1‑250 (MQ=255) cAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACCGCGTGACCCAGCACTATGCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAAGaaa                                                                                                                                                                                                                                                    <  2:45072/251‑1 (MQ=255)  aGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACCGCGTGACCCAGCACTATGCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAAGAAAc                                                                                                                                                                                                                                                   >  1:63537/1‑251 (MQ=255)              gctggcAACCAAATGATATTGCGATTTGGGATAACCGCGTGACCCAGCACTATGCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAAGAAACAAAATTAACGCa                                                                                                                                                                                                                                       <  1:98885/251‑1 (MQ=255)                   cAACCAAATGATATTGCGATTTGGGATAACCGTGTGACCCAGCACTATGCCAATGCCGATTACCTGCCAAAGCGACGGATAATGCATCGGGCGACGATCATTGGGGATCAACCGTTTTATCGGGCGGGGTAATACGAGAGCGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGACGAAACAAAATTAACGCAAAAtc                                                                                                                                                                                                                                  >  2:211911/1‑251 (MQ=255)                   cAACCAAATGATATTGCGATTTGGGATAACCGCGTGACCCAGCACTATGCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGCGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAAGAAACAAAATTAACGCAGAAtc                                                                                                                                                                                                                                  >  2:309708/1‑251 (MQ=255)                     aCCAAATGATATTGCGATTTGGGATAACCGCGTGACCCAGCACTATGCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAAGAAACAAAATTAACGCAGAAtctt                                                                                                                                                                                                                                >  1:124570/1‑251 (MQ=255)                        aaaTGATAGTGCGATTTGGGATAACCGCGTGACCCAGCACTATGCCAAGGTCGATTCACTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGTGGTAATACGAGAGTGGACGGTCCCCTCGCCACCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTGCAACAGGTTAACTCCCCCTTTCTGAGAAGATACAAAATTAACGCAGAAtcttctt                                                                                                                                                                                                                             <  2:312127/251‑1 (MQ=255)                           tGATACTGCGATTTGGGATAACCGCGTGACCCAGCACTATGCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAAGAAACAAAATTAACGCAGAATCTTCTtctc                                                                                                                                                                                                                          >  2:151935/1‑251 (MQ=255)                                  gCGATTTGGGATAACCGCGTGACCCAGCACTATGCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCTTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAAGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCaa                                                                                                                                                                                                                    <  1:20999/250‑1 (MQ=255)                                   cGATTTGGGATAACCGCGGGACCCAGCACTATGCCAATGCCGTTTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGGTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGTGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAAGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAAt                                                                                                                                                                                                                  <  2:307753/251‑1 (MQ=255)                                              gaCCGCTTGACCCAGCACTCTGCCAATGCCGATTCCCTGACACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAAGAAACAAAATTAACGCATAATCTTCTTCTCAGCCAAATCCAGCGCAAAg                                                                                                                                                                                                       <  2:281389/250‑1 (MQ=255)                                                 cgcgTGACCCAGCACTATGCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAATAGGTTAACTCCCCCTTTCTGAGAAGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAg                                                                                                                                                                                                    <  2:192109/251‑1 (MQ=255)                                                   cgTGACCCAGCACTATGCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAAGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCt                                                                                                                                                                                                  >  1:113235/1‑251 (MQ=255)                                                          cAGCACTATGCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAAGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATc                                                                                                                                                                                           >  1:146478/1‑251 (MQ=255)                                                                   gCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTGTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAAGAAACAAAATTAACGCAGAATCTTCTTCTCAGGCAAATCCAGCGCAAAGTAGCTGAAAATCAGATccgca                                                                                                                                                                                  <  2:53758/251‑1 (MQ=255)                                                                    ccAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGTTTAACTCCCCCTTTCTGAGAAGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATccgcac                                                                                                                                                                                 >  1:198184/1‑251 (MQ=255)                                                                       aTGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGCGAGTGGACGGTCCCCTCGCCACCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAAGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGAAGCTGAAAATCAGATccgcacccg                                                                                                                                                                              >  2:302019/1‑251 (MQ=255)                                                                           cGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAAGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACg                                                                                                                                                                          >  1:128035/1‑251 (MQ=255)                                                                                             aCGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAAGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCt                                                                                                                                                        <  2:42294/251‑1 (MQ=255)                                                                                               ggATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATAAGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAAGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCAACCGCACGCTTAATCGAACCTAAGCttt                                                                                                                                                      >  1:348392/1‑251 (MQ=255)                                                                                                         gcGGGCGACGAGCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAAGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGAc                                                                                                                                            <  1:194007/250‑1 (MQ=255)                                                                                                              cgacgaTCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAAGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTtct                                                                                                                                       <  2:67865/251‑1 (MQ=255)                                                                                                                                                      aTACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAAGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCACCCGCCAGCGCGGCGAACTTAatcat                                                                                               >  2:109616/1‑251 (MQ=255)                                                                                                                                                           agagTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGTTTAACTCCCCCTTTCTGAGAAGAAACAAAAGTAACGCAGAATCTTCTTCCCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGTTCCGCCCCCCCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCAt                                                                                          <  2:176243/251‑1 (MQ=255)                                                                                                                                                           agagTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAAGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCAt                                                                                          >  1:168445/1‑251 (MQ=255)                                                                                                                                                                      tCCCCTCGACCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAAGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTc                                                                               >  1:277857/1‑251 (MQ=255)                                                                                                                                                                         ccTCGCCCCCTTGGGTAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAAGAAACAAAATGAACGCAGAACCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGACCCGCCCCCGGACGTTTAGTCGAACCTAAGCTTGCGAGCACGACTTTCCCTTCATCGATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCAcc                                                                            <  2:312425/251‑1 (MQ=255)                                                                                                                                                                            cacccccTTGGGGAGAGGGGTAGGGTGAGGGGGCGTTCACCTTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAAGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGAAAAGTAGCTGAAAATCAGATCCGCCCCCGCACGCTAAATCGAACCTAAGCTTTCGAGCACAACTTTCTCTTCATCTATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGa                                                                         <  2:210456/249‑1 (MQ=255)                                                                                                                                                                                 ccTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAAGAAACAAAATTACCGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATCGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATAcgc                                                                    >  2:252063/1‑251 (MQ=255)                                                                                                                                                                                                             cGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAAGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCCCCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGc                                        <  1:193406/251‑1 (MQ=255)                                                                                                                                                                                                                      tACTTTCAACAGGTTAACTCCCCCTTTCTGAGAAGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTcacgc                               >  1:22701/1‑251 (MQ=255)                                                                                                                                                                                                                         tttCAACAGGTTAACTCCCCCTTTCTGAGACGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAACTCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCACCCGCCAGCGCGGCGCACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTcacgcacg                            >  1:244783/1‑251 (MQ=255)                                                                                                                                                                                                                          ttCAACAGGTTAACTCCCCCTTTCTGAGAAGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTCCATCTATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTcacgcac                             <  1:188756/249‑1 (MQ=255)                                                                                                                                                                                                                                  gTTAACTCCCCCTTTCTGAGAAGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGATGTCGAgg                   <  2:128035/251‑1 (MQ=255)                                                                                                                                                                                                                                  gTTAACGCCCCCTTTCGGAGAAGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAACCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAACCTTTCGAGCACGACTTTCTCTTCATCTATGGTCCCCGCCAGCGCGGCGAACTTAATCAGCGCATACTCACCGCGCACCGGATTCGCGCCAAGCGGCAATTCATTACGTTCACGCAGCTCACGCACGATGTCGAgg                   <  1:217939/251‑1 (MQ=255)                                                                                                                                                                                                                                     aaCTCCCCCTTTCTGAGAAGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGATGTCGAGGTAc                >  1:41047/1‑251 (MQ=255)                                                                                                                                                                                                                                        tCCCCCTTTCTGAGAAGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGATGTCGAGGTACGCt             >  2:176265/1‑251 (MQ=255)                                                                                                                                                                                                                                                  agagaAGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGAAAAGTAGCTGAAAATCAGATCCGCCCCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGATGTCGAGGTACGCTCCAGCAGGtt   <  2:124570/250‑1 (MQ=255)                                                                                                                                                                                                                                                   gagaAGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGATGTCGAGGTACGCTCCAGCAGGttt  >  2:303915/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       CAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACCGCGTGACCCAGCACTATGCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGGGTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGGTTAACTCCCCCTTTCTGAGAGGAAACAAAATTAACGCAGAATCTTCTTCTCAGCCAAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGCACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCTATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGATGTCGAGGTACGCTCCAGCAGGTTT  >  NC_000913/388497‑388989

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