Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,953,554	G→A	G26E (GGA→GAA)	ilvD →	dihydroxy‑acid dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,953,554	0	G	A	94.6%	118.7 / ‑2.9	38	G26E (GGA→GAA)	ilvD	dihydroxy‑acid dehydratase
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (3/32);  total (3/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

ACCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCGGCTGGCGGCGTTGCCAAAGAGTTCAACACCATTGCGGTGGATGATGGGATTGCCATGGGCCACGGGGGGATGCTTTATTCACTGCCATCTCGCGAACTGATCGCTGATTCCGTT  >  NC_000913/3953327‑3953789                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            aCCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCg                                                                                                                                                                                                                         <  2:342494/248‑1 (MQ=255)     aaCGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGAAATGACCGACGCCGATTTCGGTAAGcc                                                                                                                                                                                                                  <  1:204933/251‑1 (MQ=255)                           cGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGAAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCg                                                                                                                                                                                            <  1:310622/251‑1 (MQ=255)                              ccTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGAAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTc                                                                                                                                                                                         <  1:201109/251‑1 (MQ=255)                                             aaCCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGAAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCg                                                                                                                                                                          <  1:116911/251‑1 (MQ=255)                                                 gAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGAAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTc                                                                                                                                                                      <  2:83537/251‑1 (MQ=255)                                                       aaaTGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAAGATGGCGGGTGCTCGTGCGCTGTGGCGCTCCACCGAAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTcc                                                                                                                                                                <  1:33286/251‑1 (MQ=255)                                                       aaaTGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCATCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGAAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTcc                                                                                                                                                                <  1:86949/251‑1 (MQ=255)                                                         aTGNGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACGCTGGGAGTAAATAAAGTATGCCTAAGTGCCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCTCGTGTGCTGTGGCGCGCCACCGGAATGACCGATGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCATCCAATTTGTACCGGGTCACGTCCAt                                                                                                                                                              <  2:379641/251‑1 (MQ=255)                                                                  gTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGAAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATCTGCGCGAt                                                                                                                                                     <  2:143093/251‑1 (MQ=255)                                                                         cAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGAAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTa                                                                                                                                              <  1:23779/251‑1 (MQ=255)                                                                                aTCAATAAATACAAAAAATGGGACGTCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGAAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGt                                                                                                                                       <  2:230600/251‑1 (MQ=255)                                                                                aTCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGAAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGt                                                                                                                                       <  1:104600/251‑1 (MQ=255)                                                                                                aaTGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGAAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGaa                                                                                                                       <  2:219224/251‑1 (MQ=255)                                                                                                          cacgcacCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGAAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCGGCTggcg                                                                                                             <  1:108685/251‑1 (MQ=255)                                                                                                             gcacCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGAAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACATCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCGGCTGgcggcg                                                                                                          <  1:408698/251‑1 (MQ=255)                                                                                                                  gTCCCATTTACGAGACAGCCACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGAAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCGGCTGGCGGCGTTGcc                                                                                                     <  1:340497/251‑1 (MQ=255)                                                                                                                     ccATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGAAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCGGCTGGCGGCGTTGCCaaa                                                                                                  <  1:109330/251‑1 (MQ=255)                                                                                                                              agacagacACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGAAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCGGCTGGCGGCGTTGCCAAAGAGTTCAac                                                                                         <  1:179117/251‑1 (MQ=255)                                                                                                                                   gacaCTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGAAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCGGCTGGCGGCGTTGCCAAAGAGTTCAACACCAt                                                                                    <  1:282105/251‑1 (MQ=255)                                                                                                                                    acacTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGAAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCGGCTGGCGGCGTTGCCAAAGAGTTCAACACCAtt                                                                                   <  2:12889/251‑1 (MQ=255)                                                                                                                                     cacTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGAAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCGGCTGGCGGCGTTGCCAAAGAGTTCAACACCATTg                                                                                  <  1:89998/251‑1 (MQ=255)                                                                                                                                                   aaaGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGAAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCGGCTGGCGGCGTTGCCAAAGAGTTCAACACCATTGCGGTGGATGATggg                                                                    <  1:346633/251‑1 (MQ=255)                                                                                                                                                      gTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGAAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCGGCTGGCGGCGTTGCCAAAGAGTTCAACACCATTGCGGTGGATGATGGGAtt                                                                 <  1:161889/251‑1 (MQ=255)                                                                                                                                                              aattaccgtaccgCCACCACCACTCATTGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACAGAAATGACCGACGCCGATTTCGGTAATCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCGGCTGGCGGCGTTGCCAAAGAGTTCAACACCATTGCGGTGGATGATGGGATTGCCATgg                                                          <  2:365163/247‑1 (MQ=255)                                                                                                                                                                  aCCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGAAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCGGCTGGCGGCGTTGCCAAAGAGTTCAACACCATTGCGGTGGATGATGGGATTGCCATGGGCCa                                                      <  2:215250/250‑1 (MQ=255)                                                                                                                                                                                      cATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGAAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCGGCTGGCGGCGTTGCCAAAGAGTTCAACACCATTGCGGTGGATGATGGGATTGCCATGGGCCACGGGGGGATGCTTTATTCACt                                 <  2:202079/251‑1 (MQ=255)                                                                                                                                                                                      cATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGAAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCGGCTGGCGGCGTTGCCAAAGAGTTCAACACCATTGCGGTGGATGATGGGATTGCCATGGGCCACGGGGGGATGCTTTATTCACt                                 <  2:360842/251‑1 (MQ=255)                                                                                                                                                                                       aTGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGAAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCGGCTGGCGGCGTTGCCAAAGAGTTCAACACCATTGCGGTGGATGATGGGATTGCCATGGGCCACGGGGGGATGCTTTATTCACTg                                <  1:212687/251‑1 (MQ=255)                                                                                                                                                                                       aTGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGAAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCGGCTGGCGGCGTTGCCAAAGAGTTCAACACCATTGCGGTGGATGATGGGATTGCCATGGGCCACGGGGGGATGCTTTATTCACTg                                <  1:211651/251‑1 (MQ=255)                                                                                                                                                                                            cGTNATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGAAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCGGCTGGCGGCGTTGCCAAAGAGTTCAACACCATTGCGGTGGATGATGGGAGTGCCATGGGCCACGGGGGGATGCTTTATTCACTGCCAtc                           <  2:372212/251‑1 (MQ=255)                                                                                                                                                                                            cGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGAAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCGGCTGGCGGCGTTGCCAAAGAGTTCAACACCATTGCGGTGGATGATGGGATTGCCATGGGCCACGGGGGGATGCTTTATTCACTGCCAt                            <  1:136684/250‑1 (MQ=255)                                                                                                                                                                                               aaTATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGAAATGACCGACGCCGATTTCGGGAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCGGCTGGCGGCGTTGCCAAAGAGTTCAACACCATGGCGCTGGATGATGGGATTGCCATGGGCCGCGGGGGGATGCTTTATTCACTGCCAtctc                         <  1:174262/250‑1 (MQ=255)                                                                                                                                                                                                          tGCTCGTGCGCTGTGGCGCGCCACCGAAATGACCGACGCCGCTTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCGGCTGGCGGCGTTGCCAAAGAGTTCAACACCATTGCGGTGGATGATGGGATTGCCATGGGCCACGGGGGGATGCTTTATTCACTGCCATCTCGCGAACTGNTCg             >  2:371537/1‑251 (MQ=255)                                                                                                                                                                                                           gCTCGTGCGCTGTGGCGCGCCACCG‑AATGACCGGCGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGGAAACTGGTCGCCGAACAAATTGAAGCGGCTGGCGGCGTTGCCAAAGAGTTCAACACCATTGCGGTGGATGATGGGATTGCCATGGGCCACGGGGGGATGCTTTATTCACTGCCATCTCGCGAACTGATCGCt           <  1:131636/251‑1 (MQ=255)                                                                                                                                                                                                              cGTGCGCTGTGGCGCGCCACCGAAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCGGCTGTCGGCGTTGCCAAAGAGTTCAACACCATTGCGGTGGATGATGGGATTGCCATGGGCCACGGGGGGATGCTTTATTCACTGCCATCTCGCGAACTGATCGCTg          >  1:415776/1‑250 (MQ=255)                                                                                                                                                                                                              cGTGCGCTGTGGCGCGCCACCGAAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCGGCTGGCGGCGTTGCCAAAGAGTTCAACACCATTGCGGTGGATGATGGGATTGCCATGGGCCACGGGGGGATGCTTTATTCACTGCCATCTCGCGAACTGATCGCt           <  2:314095/249‑1 (MQ=255)                                                                                                                                                                                                                     tgtgGCGCGCCACCGAAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCGGCTGGCGGCGTTGCCAAAGAGTTCAACACCATTGCGGTGGATGATGGGATTGCCATGGGCCACGGGGGGCTGCTTTATTCACTGCCATCTCGCGGACTGATCGCTGATTCCGtt  >  2:212303/1‑251 (MQ=255)                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            ACCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATTATCGCGGTTGTGAACTCGTTCACCCAATTTGTACCGGGTCACGTCCATCTGCGCGATCTCGGTAAACTGGTCGCCGAACAAATTGAAGCGGCTGGCGGCGTTGCCAAAGAGTTCAACACCATTGCGGTGGATGATGGGATTGCCATGGGCCACGGGGGGATGCTTTATTCACTGCCATCTCGCGAACTGATCGCTGATTCCGTT  >  NC_000913/3953327‑3953789

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