Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,123,727	C→T	G27R (GGA→AGA)	glcB ←	malate synthase G

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,123,727	0	C	T	97.4%	132.5 / ‑6.5	39	G27R (GGA→AGA)	glcB	malate synthase G
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (18/20);  minor base A (0/1);  total (18/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCTGTCAATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCACTCATCAAGCGCTGCCTGAATGCGATCGCGTTCTGCCAGCAACTGACGATTTTCTGGTGCCAGATCATGAACGATCTCATCAAAATTGCGCCAGAACGCCGCAGCGTCCAGCCCTGTTCCCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAATAACGG  >  NC_000913/3123484‑3123964                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                              gcaggttggcgacgtcggTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCACTCATCAAGCGCTGCCTGAATGCGATCGCGTTCTGCCAGCAACTGACGATTTTCTGGTGCCAGATCATGAACGATCTCATCAAAATTGCGCCAGAACGCCGCAGCGTCCAGCCCTGTTCTCGGTaa                                                                                                                                                                                                                                         >  2:142433/17‑250 (MQ=255)  caggttggcgacgtcggTTTCCCCCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCTCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCACTCATCAAGCGCTGCCTGAATGCGATCGCGTTCTGCCAGCAACTGACGATTTTCTGGTGCCAGATCATGAACGATCTCATCAAAATTGCGCCAGAACGCCGCAGCGTCCAGCCCTGTTCTCGGTaaaa                                                                                                                                                                                                                                       <  1:142433/236‑1 (MQ=255)   cTGTCAATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCACTCATCAAGCGCTGCCTGAATGCGATCGCGTTCTGCCAGCAACTGACGATTTTCTGGTGCCAGATCATGAACGATCTCATCAAAATTGCGCCAGAACGCCGCAGCGTCCAGCCCTGTTCTCGGTAAAAc                                                                                                                                                                                                                                      >  1:88675/1‑251 (MQ=255)                         ccGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCACTCATCAAGCGCTGCCTGAATGCGATCGCGTTCTGCCAGCAACTGACGATTTTCTGGTGCCAGATCATGAACGATCTCATCAAAATTGCGCCAGAACGCCGCAGCGTCCAGCCCTGTTCTCGGTAAAACTTCTTCATCCACAAAACGTTTa                                                                                                                                                                                                                <  1:241886/251‑1 (MQ=255)                           gTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCACTCATCAAGCGCTGCCTGAATGCGATCGCGTTCTGCCAGCAACTGACGATTTGCTGGTGCCAGATCATGAACGATCTCATCAAAATTGCGCCAGAACGCCGCAGCGTCCAGCCCTGTTCTCGGTAAAACTTCTTCATCCACAAAACGTTTaaa                                                                                                                                                                                                              <  1:201474/251‑1 (MQ=255)                                                               tCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCACTCATCAAGCGCTGCCTGAATGCGATCGCGTTCTGCCAGCAACTGACGATTTTCTGGTGCCAGATCATGAACGATCTCATCAAAATTGCGCCAGAACGCCGCAGCGTCCAGCCCTGTTCTCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATgg                                                                                                                                                                           <  1:145960/250‑1 (MQ=255)                                                                     aGGAAAGATTTTTAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCACTCATCAAGCGCTGCCTGAATGCGATCGCGTTCTGCCAGCAACTGACGATTTTCTGGTGCCAGATCATGAACGATCTCATCAAAATTGCGCCAGAACGCCGCAGCGTCCAGCCCTGTTCTCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGCCTCTGGGTTATGGTTTGACt                                                                                                                                                                    >  2:162489/1‑251 (MQ=255)                                                                      ggAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCACTCATCAAGCGCTGCCTGAATGCGATCGCGTTCTGCCAGCAACTGACGATTTTCTGGTGCCAGATCATGAACGATCTCATCAAAATTGCGCCAGAACGCCGCAGCGTCCAGCCCTGTTCTCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTc                                                                                                                                                                   >  2:120403/1‑251 (MQ=255)                                                                       gAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCACTCATCAAGCGCTGCCTGAATGCGATCGCGTTCTGCCAGCAACTGACGATTTTCTGGTGCCAGATCATGAACGATCTCATCAAAATTGCGCCAGAACGCCGCAGCGTCCAGCCCTGTTCTCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCa                                                                                                                                                                  <  1:230304/251‑1 (MQ=255)                                                                         aaGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCACTCATCAAGCGCTGCCTGAATGCGATCGCGTTCTGCCAGCAACTGACGATTTTCTGGTGCCAGATCATGAACGATCTCATCAAAATTGCGCCAGAACGCCGCAGCGTCCAGCCCTGTTCTCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCAtt                                                                                                                                                                <  1:158472/251‑1 (MQ=255)                                                                         aaGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCACTCATCAAGCGCTGCCTGAATGCGATCGCGTTCTGCCAGCAACTGACGATTTTCTGGTGCCAGATCATGAACGATCTCATCAAAATTGCGCCAGAACGCCGCAGCGTCCAGCCCTGTTCTCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCAtt                                                                                                                                                                >  1:98464/1‑251 (MQ=255)                                                                                  aGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCACTCATCAAGCGCTGCCTGAATGCGATCGCGTTCTGCCAGCAACTGACGATTTTCTGGTGCCAGATCATGAACGATCTCATCAAAATTGCGCCAGAACGCCGCAGCGTCCAGCCCTGTTCTCGGTAAAACTTCTTCATCCACAAAGCGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTAtct                                                                                                                                                        <  1:3106/250‑1 (MQ=255)                                                                                          ttATCTTTTACCGGCCCCGGATTGCTGCGATGCCACTCATCAAGCGCTGCCTGAATGCGATCGCGTTCTGCCAGCAACTGACGATTTTCTGGTGCCAGATCATGAACGATCTCATCAAAATTGCGCCAGAACGCCGCAGCGTCCAGCCCTGTTCTCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGtttt                                                                                                                                               >  1:71783/1‑251 (MQ=255)                                                                                                tttACCGGCCCCGGATTGCTGCGATGCCACTCATCAAGCGCTGCCTGAATGCGATCGCGTTCTGCCAGCAACTGACGATTTTCTGGTGCCAGATCATGAACGATCTCATCAAAATTGCGCCAGAACGCCGCAGCGTCCAGCCCTGTTCTCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTTACTCATTGTTTATCTCCTCGTTTTCGCt                                                                                                                                           <  1:108127/249‑1 (MQ=255)                                                                                                    ccGGCCCCGGCTTGCTGCGATGCCACTCATCAAGCGCTGCATGAATGCGATCGCGTTCTGCCAGCAACTGACGATTTTCTGGTGCCAGATCATGAACGATCTCATCAAAATTGCGCCAGAACGCCGCAGCGTCCAGCCCTGTTCTCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCTTCGTTTTCGCTTATTTc                                                                                                                                     >  2:248715/1‑251 (MQ=255)                                                                                                            ggATTGCTGCGATGCCACTCATCAAGCGCTGCCTGAATGCGATCGCGTTCTGCCAGCAACTGACGATTTTCTGGTGCCAGATCATGAACGATCTCATCAAAATTGCGCCAGAACGCCGCAGCGTCCAGCCCTGTTCTCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAacac                                                                                                                             >  1:64427/1‑251 (MQ=255)                                                                                                              aTTGCTGCGATGCCACTCATCAAGCGCTGCCTGAATGCGATCGCGTTCTGCCAGCAACTGACGATTTTCTGGTGCCAGATCATGAACGATCTCATCAAAATTGCGCCAGAACGCCGCAGCGTCCAGCCCTGTTCTCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCg                                                                                                                           >  1:2064/1‑251 (MQ=255)                                                                                                                       aTGCCACTCATCAAGCGCTGCCTGAATGCGATCGCGTTCTGCCAGCAACTGACGATTTTCTGGTGCCAGATCATGAACGATCTCATCAAAATTGCGCCAGAACGCCGCAGCGTCCAGCCCTGTTCTCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCg                                                                                                                  <  2:88675/251‑1 (MQ=255)                                                                                                                          ccACGCATCAAGCGCTGCCTGAATGCGATCGCGTTCTGCCAGCAACTGACGATTTTCTGGTGCCAGATCATGAGCGATCTCATCAAAATTGCGCCAGAACGCCGCAGCGTCCAGCCCTGTTCTCGGTAAAACTTCTTCATCCTCAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACGCATTGTTTATCTCCTCGTTTTCGCTTATTTGGCCAACACCGCTGCTGCCGCtt                                                                                                               <  2:56886/251‑1 (MQ=255)                                                                                                                                            cTGAATGCGATCGCGTTCTGCCAGCAACTGACGATTTTCTGGTGCCAGATCATGAACGATCTCATCAAAATTGCGCCAGAACGCCGCAGCGTCCAGCCCTGTTCTCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGCCCTGCGCATCCt                                                                                             >  2:109215/1‑251 (MQ=255)                                                                                                                                              gAATGCGATCGCGTTCTGCCAGCAACTGACGATTTTCTGGTGCCAGATCATGAACGATCTCATCAAAATTGCGCCAGAACGCCGCAGCGTCCAGCCCTGTTCTCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCtgt                                                                                           <  1:7209/251‑1 (MQ=255)                                                                                                                                                    gATCGCGTTCTGCCAGCAACTGACGATTTTCTGGTGCCAGATCATGAACGATCTCATCAAAATTGCGCCAGAACGCCGCAGCGTCCAGCCCTGTTCTCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCg                                                                                     >  2:198096/1‑251 (MQ=255)                                                                                                                                                     aTCGCGTTCTGCCAGCAACTGACGATTTTCTGGTGCCAGATCATGAACGATCTCATCAAAATTGCGCCAGAACGCCGCAGCGTCCAGCCCTGTTCTCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCgg                                                                                    <  1:262810/251‑1 (MQ=255)                                                                                                                                                        gcgTTCTGCCAGCAACTGACGATTTTCTGGTGCCAGATCATGAACGATCTCATCAAAATTGCGCCAGAACGCCGCAGCGTCCAGCCCTGTTCTCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTa                                                                                 <  1:44441/251‑1 (MQ=255)                                                                                                                                                          gTTCTGCCAGCAACTGACGATTTTCTGGTGCCAGATCATTAACGATCTCATCAAAATTGCGCCAGAACGCCGCAGCGTCCAGCCCTGTTCTCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTaaa                                                                               <  2:98464/251‑1 (MQ=255)                                                                                                                                                               gCCAGCAACTGACGATTTTCTGGTGCCAGATCATGAACGATCTCATCAAAATTGCGCCAGAACGCCGCAGCGTCCAGCCCTGTTCTCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAAcc                                                                             <  1:99565/248‑1 (MQ=255)                                                                                                                                                               gCCAGCAACTGACGATTTTCTGGGTCCAGATCATGAACGATCTCATCAAAATTGCGCCAGAACGCCGCAGCGTCCAGCCCTGTTCTCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAAcc                                                                             <  2:2064/248‑1 (MQ=255)                                                                                                                                                                                  cacgTGCCAGAGCATGAACGATCTCATCAAAATTGCGCCAGAACGCCGCAGCGTCCATCGCTGTTCTCGGAAAAACGTCTTCAACCACAAAGCGTTTAAAATTGGCGGCAATGCGTAAACGGCTCTGGGTTATGGTGTGACTCATGGTTGACCTCCTCGTTTTCGCTTATTTCGCCAACACCGTTGCTGCCGCTTTCGCGACCTGCGCATCCCGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCaat                                                        <  1:161875/247‑1 (MQ=255)                                                                                                                                                                                                   aCGATCTCATCAAAATTGCGCCAGAACGCCGCAGCGTCCAGCCCTGTTCTCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCacaac                                      >  1:113035/1‑251 (MQ=255)                                                                                                                                                                                                      aTCTCATCAAAATTGCGCCAGAACGCCGCAGCGTCCAGCCCTGTTCTCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACaacaac                                   >  2:79154/1‑251 (MQ=255)                                                                                                                                                                                                         tcatcaAAATTGCGCCAGAACGCCGCAGCGTCCAGCCCTGTTCTCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCgg                                >  2:180129/1‑251 (MQ=255)                                                                                                                                                                                                          cagcaAAATGGCGCCAGAACGCCGCAGCGTCCAGCCCTGTTCTCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCGGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCacaaca                                     <  2:168971/242‑1 (MQ=255)                                                                                                                                                                                                               aaaTTGCGCCAGAACGCCGCAGCGTCCAGCCCTGTTCTCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTAcgcc                          >  1:23410/1‑251 (MQ=255)                                                                                                                                                                                                               aaaTTGCGCCAGAACGCCGCAGCGTCCAGCCCTGTTCTCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTAcgcc                          >  1:148645/1‑251 (MQ=255)                                                                                                                                                                                                               aaaTTGCGCCAGAACGCCGCAGCGTCCAGCCCTGTTCTCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTATCGCGACCTGCGCATCCTGTGCGCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTAcgcc                          >  2:93127/1‑251 (MQ=255)                                                                                                                                                                                                                    gcgcCAGAACGCCGCAGGGGCCAGCCCTGTTCACGGTAAAACTTCTTCAGCCACAAAACTTTTAAAATTGGCGTCAATGCGTAAACGGCGCTGGGTTATGGTTTGACTCAGTGTTTATCTCCTCGTTTTCGCTTATTTCGCGAGCACCGCTGCTGCCGCTTTCGCGAGCTGCGCATCCTTTGCTCCGGTTAAACCAGTAACGCCCACGGGACCAATAATTTGCCCATCCACAACAACCGGTTCTCCGCCtt                     <  2:277291/251‑1 (MQ=255)                                                                                                                                                                                                                         aGAACGCCGCAGCGTCCAGCCCTGTTCTCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACGCCGCCTTCCAGc                <  1:108265/251‑1 (MQ=255)                                                                                                                                                                                                                               ccgcAGCGTCCAGCCCTGTTCTCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGtt          <  1:273529/251‑1 (MQ=255)                                                                                                                                                                                                                                       tCCAGCCCTGTTCTCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAACAATTTGCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAATAACgg  >  2:5651/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                              CGCTGTCAATGCCCGTGGTTTCCACCGTCACGCGCTCCGGTTGCGGCACCAGGTAGCCCAGTTCACGCAGGAAAGATTTATAGGCCGCTTTATCTTTTACCGGCCCCGGATTGCTGCGATGCCACTCATCAAGCGCTGCCTGAATGCGATCGCGTTCTGCCAGCAACTGACGATTTTCTGGTGCCAGATCATGAACGATCTCATCAAAATTGCGCCAGAACGCCGCAGCGTCCAGCCCTGTTCCCGGTAAAACTTCTTCATCCACAAAACGTTTAAAATTGGCGTCAATGCGTAAACGGCTCTGGGTTATGGTTTGACTCATTGTTTATCTCCTCGTTTTCGCTTATTTCGCCAACACCGCTGCTGCCGCTTTCGCGACCTGCGCATCCTGTGCTCCGGTTAAACCAGAAACGCCCACGGCACCAATAATTTGCCCATCCACAACAACCGGTACGCCGCCTTCCAGCGACGTTAATAACGG  >  NC_000913/3123484‑3123964

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