Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,792,007	C→G	G182A (GGT→GCT)	kbl ←	2‑amino‑3‑ketobutyrate CoA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,792,007	0	C	G	97.8%	148.3 / NA	46	G182A (GGT→GCT)	kbl	2‑amino‑3‑ketobutyrate CoA ligase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (18/27);  minor base A (1/0);  total (19/27)
Fisher's exact test for biased strand distribution p-value = 4.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

AGCCACTCAACCACTTCTTTGCGCGCCGCGGTATAACCACCAGAAGCCCCGCCCAGCGCTTTACCAAGCGTACCGGTGATAATATCGACCCGGCCCATCACATCGCAGTATTCATGGGAACCACGACCATTTTCACCGACAAAACCGACCGCGTGGGAGTCGTCTACCATCACCA‑GGGCATCATATTTATCTGCCAGATCGCAAACGCCCTTCAGGTTGGCAATCACGCCGTCCATTGAGAACACACCATCGGTGGCGATCAGCACATGACGCGCACCGGCTTCACGCGCTTCTTTCAGACGTGCTTCCAGCTCCTGCATATCGTTGTTGGCATAGCGATAGCGTTTAGCTTTGCACAGACGCACACCATCAATAATAGACGCGTGGTTCAGTGCGTCGGAGATAATGGCGTCTTCCGCACCCAGAAGCGTTTCAAACAGGCCACCGTTAGCATCAAAGCAGGAAGAGTAGAGAATCGCATC  >  NC_000913/3791761‑3792242                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                            cacccacccaaccctgcttTGCGCGCCGCGGTATAACCACCAGAAGCCCCGCCCAGCGCTTTACCAAGCGTACCGGTGATAATATCGACCCGGCCCATCACATCGCAGTATTCATGGGAACCACGACCATTTTCACCGACAAAACCGACCGCGTGGGAGTCGTCGACCATCACCA‑GGGCATCATATTTATCTGCCAGATCGCAAACGCCCGTCAGGTTGGCAATCACGCCGTCCATTGAGAACACAGCATc                                                                                                                                                                                                                                         <  2:222301/235‑1 (MQ=255)           ccACTTCTTTGCGCGCCGCGGTATAACCACCAGAAGCCCCGCCCAGCGCTTTACCAAGCGTACCGGTGATAATATCGACCCGGCCCATCACATCGCAGTATTCATGGGAACCACGACCATTTTCACCGACAAAACCGACCGCGTGGGAGTCGTCTACCATCACCA‑GGGCATCATATTTATCTGCCAGATCGCAAACGCCCTTCAGGTTGGCAATCACGCCGTCCATTGAGAACACAGCATCGGTGGCGATc                                                                                                                                                                                                                               <  1:112357/251‑1 (MQ=255)                        cgccgcGGTATAACCACCAGAAGCCCCGCCCAGCGCTTTACCAAGCGTACCGGTGATAATATCGACCCGGCCCATCACATCGCAGTATTCATGGGAACCACGACCATTTTCACCGACAAAACCGACCGCGTGGGAGTCGTCTACCATCACCA‑GGGCATCATATTTATCTGCCAGATCGCAAACGCCCTTCAGGTTGGCAATCACGCCGTCCATTGAGAACACAGCATCGGTGGCGATCAGCACATGAcgcg                                                                                                                                                                                                                  >  1:313285/1‑251 (MQ=255)                                                cccGCCCAGCGCTTTACCCAGAGTACCGGTGATCATATCGACCCGGCACATCACCTCGCAGTATTCATGGGAACCACGAACATTTTCACCGACAAAACCGACCGCGAGGGAGTCGTCTACCATCACCA‑GGGCATCATATTTATCTGCCAGATCGCAAACGCCCTTCAGGTTGGCAATCACGCCGTCCCTTGAGAACACAGCATCGGTg                                                                                                                                                                                                                                     >  2:146155/1‑208 (MQ=255)                                                cccGCCCAGCGCTTTACCAAGCGTACCGGTGATAATATCGACCCGGCCCATCACATCGCAGTATTCATGGGAACCACGACCATTTTCACCGACAAAACCGACCGCGTGGGAGTCGTCTACCATCACCA‑GGGCATCATATTGATCTGCCAGATCGCAAACGCCCTGCAGGGTGGCAATCACGCCGTCCATTGAGAACACAGCATCGGTg                                                                                                                                                                                                                                     <  1:146155/208‑1 (MQ=255)                                                          gcTTTACCAAGCGTACCGGTGATAATATCGACCCGGCCCATCACATCGCAGTATTCATGGGAACCACGACCATTTTCACCGACAAAACCGACCGCGTGGGAGTCGTCTACCATCACCA‑GGGCATCATATTTATCTGCCAGATCGCAAACGCCCTTCAGGTTGGCAATCACGCCGTCCATTGAGAACACAGCATCGGTGGCGATCAGCACATGACGCGCACCGGCTTCACGCGCTTCTTTCAGACGTGCTTc                                                                                                                                                                                >  1:176707/1‑251 (MQ=255)                                                          gcTTTACCAAGCGTACCGGTGATAATATCGACCCGGCCCATCACATCGCAGTATTCATGGGAACCACGACCATTTTCACCGACAAAACCGACCGCGTGGGAGTCGTCTACCATCACCA‑GGGCATCATATTTATCTGCCAGATCGCAAACGCCCTTCAGGTTGGCAATCACGCCGTCCATTGAGAACACAGCATCGGTGGCGATCAGCACATGACGCGCACCGGCTTCACGCGCTTCTTGCAGACGTGCTTc                                                                                                                                                                                <  2:160325/251‑1 (MQ=255)                                                             ttACCAAGCTAACCGGTGATAAGACCGACCCGGCCCATCACATCGCAGTATTCATGGGAACCACGACCATTTTCACCGACAAAACCGACCGCGTGGGAGTCGTCCACCATCGCCA‑GGGCATCATATTTGTCTGCCAGATCGCAAACGCCCTTCAGGTTGGCAATCACGCCGTCCATTGAGAACACAGCATCGGTGGCGATCAGCACATGACGCGCACCGGCTTCACGCGCTTCTTTCAGACGTGCTTcc                                                                                                                                                                               <  2:4610/249‑1 (MQ=255)                                                                ccAAGCGTACCGGTGATAATATCGACCCGGCCCATCACATCGCAGTATTCATGGGAACCACGACCATTTTCACCGACAAAACCGACCGCGTGGGAGTCGTCTACCATCACCA‑GGGCATCATATTTATCTGCCAGATCGCAAACGCCCTTCAGGTTGGCAATCACGCCGTCCATTGAGAACACAGCATCGGTGGCGATCAGCACATGACGCGCACCGGCTTCACGCGCTTCTTTCAGACGTGCTTCCAGCTc                                                                                                                                                                          >  1:3004/1‑251 (MQ=255)                                                                     cGTACCGGGGATAATATCGACCCGGCCCATCACATCGCAGTATTCATGGGAACCACGACCATTTTCACCGACAAAACCGACCGCGTGGGAGTCGTCTACCATCACCA‑GGGCATCATATTTATCTGCCAGATCGCAAACGCCCTTCAGGTTGGCAATCACGCCGTCCATTGAGAACACAGCATCGGTGGCGATCAGCACATGACGCGCACCGGCTTCACGCGCTTCTTTCAGACGTGCTTCCAGCTCCTGCa                                                                                                                                                                     <  2:138870/251‑1 (MQ=255)                                                                                       gACCCGGCCCATCACATCGCAGTATTCATGGGAACCACGACCATTTTCACCGACAAAACCGACCGCGTGGGAGTCGTCTACCATCACCA‑GGGCATCATCTTTATCTGCCAGATCGCAAACGCCCTTCAGGTTGGCAATCACGCCGTCCATTGAGAACACAGCATCGGTGGCGATCAGCACATGACGCGCACCGGCTTCACGCGCTTCTTTc                                                                                                                                                                                           >  2:348815/1‑211 (MQ=255)                                                                                       gACCCGGCCCATCACATCGCAGTATTCATGGGAACCACGACCATTTTCACCGACAAAACCGACCGCGTGGGAGTCGTCTACCATCACCA‑GGGCATCATATTTATCTGCCAGATCGCAAACGCCCTTCAGGTTGGCAATCACGCCGTCCATTGAGAACACAGCATCGGTGGCGATCAGCACATGACGCGCACCGGCTTCACGCGCTTCTTTc                                                                                                                                                                                           <  1:348815/211‑1 (MQ=255)                                                                                        aCCCGGCCCATCACATCGCAGTATTCATGGGAACCACGACCATTTTCACCGACAAAACCGACCGCGTGGGAGTCGTCTACCATCACCA‑GGGCATCATATTTATCTGCCAGATCGCAAACGCCCTTCAGGTTGGCAATCACGCCGTCCATTGAGAACACAGCATCGGTGGCGATCAGCACATGACGCGCACCGGCTTCACGCGCTTCTTTCAGACGTGCTTCCAGCTCCTGCATATCGTTGTTGGCatagcg                                                                                                                                                  <  1:193407/251‑1 (MQ=255)                                                                                           cGGCCCATCACATCGCAGTATTCATGGGAACCACGACCATTTTCACCGACACAACCGACCGCGTGGGAGTCGTCTACCATCACCA‑GGGCATCATATTTATCTGCCAGATCGCAAACGCCCTTCAGGTTGGCAATCACGCCGTCCATTGAGAACACAGCATCGGTGGCGATCAGCACATGACGCGCACCGGCTTCACGCGCTTCTTTCAGACGTGCTTCCAGCTCCTGCATATCGTTGTTGGCatagcgat                                                                                                                                                >  2:37466/1‑250 (MQ=255)                                                                                             gCCCATCACATCGCAGTATTCATGGGAACCACGACCATTTTCACCGACAAAACCGACCGCGTGGGAGTCGTCTACCATCACCA‑GGGCATCATATTTATCTGCCAGATCGCAAACGCCCTTCAGGTTGGCAATCACGCCGTCCATTGAGAACACAGCATCGGTGGCGATCAGCACATGACGCGCACCGGCTTCACGCGCTTCTTTCAGACGTGCTTCCAGCTCCTGCATATCGTTGTTGGCatagcgatagc                                                                                                                                             >  1:16795/1‑251 (MQ=255)                                                                                             gCCCATCACATCGCAGTATTCATGGGAACCACGACCATTTTCACCGACAAAACCGACCGCGTGGGAGTCGTCTACCATCACCA‑GGGCATCATATTTATCTGCCAGATCGCAAACGCCCTTCAGGTTGGCAATCACGCCGTCCATTGAGAACACAGCATCGGTGGCGATCAGCACATGACGCGCACCGGCTTCACGCGCTTCTTTCAGACGTGCTTCCAGCTCCTGCATATCGTTGTTGGCatagcgatagc                                                                                                                                             <  1:157771/251‑1 (MQ=255)                                                                                                       tCGCAGTATTCATGGGAACCACGACCATTTTCACCGACAAAACCGACCGCGTGGGAGTCGTCTACCATCACCA‑GGGCATCATATTTATCTGCCAGATCGCAAACGCCCTTAAGGTTGGCAATCACGCCGTCCATTGAGAACACAGCATCGGTGGCGATCAGCACATGACGCGCACCGGCTTCACGCGCTTCTTTCAGACGTGCTTCCAGCTCCTGCATATCGTTGTTGGCATAGCGATAGCGTTTAGCttt                                                                                                                                   >  1:138982/1‑251 (MQ=255)                                                                                                              aTTCATGGGAACCACGACCATTTTCACCGACAAAACCGACCGCGTGGGAGTCGTCTACCATCACCA‑GGGCATCATATTTATCTGCCAGATCGCAAACGCCCTTCAGGTTGGCAATCACGCCGTCCATTGAGAACACAGCATCGGTGGCGATCAGCACATGACGCGCACCGGCTTCACGCGCTTCTTGCAGACGGGCTTCCAGCTCCTGCATATCGTTGTTGGCATAGCGATAGCGTTTAGCTTTGCACAg                                                                                                                             <  2:175136/250‑1 (MQ=255)                                                                                                                          cACGACCATTTTCACCGACAAAACCGACCGCGTGGGAGTCGTCTACCATCACCA‑GGGCATCATATTTATCTGCCAGATCGCAAACGCCCTTCAGGTTGGCAATCACGCCGTCCATTGAGAACACAGCATCGGTGGCGATCAGCACATGACGCGCACCGGCTTCACGCGCTTCTTTCAGACGTGCTTCCAGCTCCTGCATATCGTTGTTGGCATAGCGATAGCGTTTAGCTTTGCACAGACGCACACCATCa                                                                                                                <  1:213593/251‑1 (MQ=255)                                                                                                                             gACCATTTACACCGACAAAACCGACCGCGTGGGAGTCGTCTACCATCACCA‑GGGCATCATATTTATCTGCCAGATCGCAAACGCCCTTCAGGTTGGCAATCACGCCGTCCATTGAGAACACAGCATCGGTGGCGATCAGCACATGACGCGCACCGGCTTCACGCGCTTCTTTCAGACGTGCTTCCAGCTCCTGCATATCGTTGTTGGCATAGCGATAGCGTTTAGCTTTGCACAGACGCACACCATCaata                                                                                                             <  1:55205/251‑1 (MQ=255)                                                                                                                                 aTTTTCACCGACAAAACCGACCGCGTGGGAGTCGCCTACCATC‑CCACGGGAATCATATTTATCTGCCAGCTCGCAAACGCCCTTCAGGTTGGCAATCACGCCGTCCCTTGAGAACACAGCATCGGTGGCGATCAGCACATGACGCGCACCGGCTTCACGCGCTTCTTTCAGACGTGCTTCCAGCTCCTGAATATCGTTGTTGGCATAGCGCTAGCGTTTAGCTTTGCACAGACGCACACCATCaataat                                                                                                           >  1:327273/1‑249 (MQ=255)                                                                                                                                    ttCACCGACAAAACCGACCGCGTGGGAGTCGTCTACCATCACCA‑GGGCATCATATTTATCTGCCAGATCGCAAACGCCCTTCAGGTTGGCAATCACGCCGTCCATTGAGAACACAGCATCGGTGGCGATCAGCACATGACGCGCACCGGCTTCACGCGCTTCTTTCAGACGTGCTTCCAGCTCCTGCATATCGTTGTTGGCATAGCGATAGCGTTTAGCTTTGCACAGACGCACACCATCAATAATAGAcg                                                                                                      >  1:240400/1‑251 (MQ=255)                                                                                                                                     tCACCGACAAAACCGACCGCGGGGGAGTCGTCTACCACCACCA‑GGGCATCATATTTATCTGCCAGATCGCAAACGCCCTTCAGGTTGGCAATCACGCCGTCCATTGAGAACACAGCAGCGGTGGCGATCAGCACATGACGCGCACCGGCTTCACGCGCTTCTTTCAGACGTGCTTCCAGCTCCTGCATATCGTTGTTGGCATAGCGATAGCGTTTAGCTTTGCACAGACGCACACCATCAATAATAGAcg                                                                                                      <  2:118849/250‑1 (MQ=255)                                                                                                                                        ccGACAAAACCGACCGCGTGGGAGTCGTCTACCATCACCA‑GGGCATCATATTTATCTGCCAGAGCGCAAACGCCCTTCAGGTTGGCAATCACGCCGTCCAGTGAGAACACAGCAGCGGTGGCGATCAGCACAGGACGCGCACCGGCTTCACGCGCTTCTTTCAGACGTGCTGCCAGCTCCTGCATATCGTTGTTGGCATAGCGATAGCGTTTAGCTTTGCACAGACGCACACCATCAATAATAGACGCGTg                                                                                                  <  2:138982/251‑1 (MQ=255)                                                                                                                                         cGACAAAACCGACCGCGTGGGAGTCGTCTACCATCACCA‑GGGCATCATATTTATCTGCCAGATCGCAAACGCCCTTCAGGTTGGCAATCACGCCGTCCATTGAGAACACAGCATCGGTGGCGATCAGCACATGACGCGCACCGGCTTCACGCGCTTCTTTCAGACGTGCTTCCAGCTCCTGCATATCGTTGTTGGCATAGCGATAGCGTTTAGCTTTGCACAGACGCACACCATCAATAATAGACGCGTgg                                                                                                 <  1:70533/251‑1 (MQ=255)                                                                                                                                                                      aCCATCACCA‑GGGCATCATATTTATCTGCCAGATCGCAAACGCCCTTCAGGTTGGCAATCACGCCGTCCATTGAGAACACAGCATCGGTGGCGATCAGCACATGACGCGCACCGGCTTCACGCGCTTCTTTCAGACGTGCTTCCAGCTCCTGCATATCGTTGTTGGCATAGCGATAGCGTTTAGCTTTGCACAGACGCACACCATCAATAATAGACGCGTGGTTCAGTGCGTCGGAGATAATGGCGTCTTc                                                                    <  1:364673/251‑1 (MQ=255)                                                                                                                                                                      aCCATCACCA‑GGGCATCATATTTATCGGCCAGATCGCAATCGCCCTTCAGGTTGGCAATCACGCCGGCCATTGATAACACAGCATCGGTGGCGATCAGCACATGACGCGCACCGGCATCACGCGCTTCTTGCAGACGTGCTTCCAGCTCCTGCATATCGTTGTTGGCATAGCGATAGCGTTTAGCTTTGCACAGACGCACACCATCAATAATAGACGCGTGGTTCAGTGCGTCTGAGATAATGGCGTCt                                                                      <  2:132037/249‑1 (MQ=255)                                                                                                                                                                            aCCA‑GGGCATCATATTTATCTGCCAGATCGCAAACGCCCTTCAGGTTGGCAATCACGCCGTCCATTGAGAACACAGCATCGGTGGCGATCAGCACATGACGCGCACCGGCTTCACGCGCTTCTTTCAGACGTGCTTCCAGCTCCTGCATATCGTTGTTGGCATAGCGATAGCGTTTAGCTTTGCACAGACGCACACCATCAATAATAGACCCGTGGTTCAGTGCGTCGGAGATAATGGCGTCTTCCGCAcc                                                              >  2:258277/1‑251 (MQ=255)                                                                                                                                                                                        atatTTATCTGCCAGATCGCAAACGCCCTTCAGGTTGGCAATCACGCCGTCCATTGAGAACACAGCATCGGTGGCGATCAGCACATGACGCGCACCGGCTTc                                                                                                                                                                                                        >  1:205744/1‑102 (MQ=255)                                                                                                                                                                                        atatTTATCTGCCAGATCGCAAACGCCCTTCAGGTTGGCAATCACGCCGTCCATTGAGAACACAGCATCGGTGGCGATCAGCACATGACGCGCACCGGCTTc                                                                                                                                                                                                        <  2:205744/102‑1 (MQ=255)                                                                                                                                                                                             tATCTGCCAGATCGCAAACGCCCTTCAGGTTGGCAATCACGCCGTCCATTGAGAACACAGCATCGGTGGCGATCAGCACATGACGCGCACCGGCTTCACGCGCTTCTTTCAGACGTGCTTCCAGCGCCTGCATATCGTTGTTGGCATAGCGATAGCGTTTAGCTTTGCACAGACGCACACCATCAATAATAGACGCGTGGTTCAGTGCGTCGGAGATAATGGCGTCTTCCGCACCCAGAAGCGTTTCaaa                                               <  2:135098/250‑1 (MQ=255)                                                                                                                                                                                              aTCTGCCAGATCGCAAGCGCCCTTCAGGTTGGCAATCACGCCGTCCATTGAGAACACAGCATCGGTGGCGATCAGCACATGACGCGCACCGGCTTCACGCGCTTCTTTCAGACGTGCTTCCAGCTCCTGCATATCGTTGTTGGCATAGCGATAGCGTTTAGCTTTGCACAGACGCACACCATCAATAATAGACGCGTGGTTCAGTGCGTCGGTGATAATGGCGTCTTCCGCACCCAGAAGCGTTTCAAAc                                              <  2:204846/250‑1 (MQ=255)                                                                                                                                                                                                      gATCGCAAACGCCCTTCAGGTTGGCAATCACGCCGTCCATTGAGAACACAGCATCGGTGGCGATCAGCACATGACGCGCACCGGCTTCACGCGCTTCTTTCAGACGTGCTTCCAGCTCCTGCATATCGTTGTTGGCATAGCGATAGCGTTTAGCTTTGCACAGACGCACACCATCAATAATAGACGCGTGGTTCAGTGCGTCGGAGATAATGGCGTCTTCCGCACCCAGAAGCGTTTCAAACAGGCCACCg                                     <  2:143679/251‑1 (MQ=255)                                                                                                                                                                                                          gCAAACGCCCTTCAGGTTGGCAATCACGCCGTCCATTGAGAACACAGCATCGGTGGCGATCAGCACATGACGCGCACCGGCTTCACGCGCTTCTTTCAGACGTGCTTCCAGCTCCTGCATATCGTTGTTGGCa                                                                                                                                                       >  1:293461/1‑133 (MQ=255)                                                                                                                                                                                                          gCAAACGCCCTTCAGGTTGGCAATCACGCCGTCCATTGAGAACACAGCATCGGTGGCGATCAGCACATGACGCGCACCGGCTTCACGCGCTTCTTTCAGACGTGCTTCCAGCTCCTGCATATCGTTGTTGGCa                                                                                                                                                       <  2:293461/133‑1 (MQ=255)                                                                                                                                                                                                                   cttcaactttgcactcccGCCGTCCACTGCGAAACCAGCATCCGTGACGACCAGCACCTGACGCGCCCCGGCTTCACCACCTTCTCTCAGACGTGCTTCCAGCTCCTGCATATCATTGTTCGCATAGCGATAGCGTTCAGCTTTGCGCAAACGAACACCATCAATAATCGACGCGTGGTTCAGTGCGTCGGAGATAATGGCGTCTTCCGCACCCAGAAGCGTTTCATACAGGCCACCGTTAGCATCaaatc                        >  2:163872/18‑249 (MQ=255)                                                                                                                                                                                                                   cTTCAGGTTGGCAATCACGCCGTCCATTGAGAACACAGCATCGGTGGCGATCAGCACATGACGCGCACCGGCTTCACGCGCTTCTTTCAGACGTGCTTCCAGCTCCTGCATATCGTTGTTGGCATAGCGATAGCGTTTAGCTTTGCACAGACGCACACCATCAATAATAGACGCGTGGTTCAGTGCGTCGGAGATAATGGCGTCTTCCGCACCCAGAAGCGTTTCAAACAGGCCACCGTTAGCATCAAAGc                        >  1:310659/1‑251 (MQ=255)                                                                                                                                                                                                                        gggtGGCAATCACGCCGTCAAGTGAGAACACAGCATCGGTGGCGATCAGCACATGCAGCGCACCGGCTTCACGCGCTTCTGTCAGACGGGCTTCCAGCTCCTGCATATCGTTGTTGGCATAGCGATAGCGTTTAGCTTTGCACAGACGCACACCATCAATAATAGACGCGTGGTTCAGTGCGTCGGAGATAATGGCGTCTTCCGCACCCAGAGGCGGTTCAAACAGGCCACCGTTAGCATCAAAGCAGGaa                   <  2:148285/248‑1 (MQ=255)                                                                                                                                                                                                                           tGGCATTCACGCCGTCCATTGAGAACACAGCATCGGTGGCGATCAGCTCATGACGCGCACCGGCTTCACGCGCTTCTTTCAGACGTGCTTCCAGCTCCTGCATATCGTTGTTGGCATAGCGATAGCGTTTAGCTTTGCACAGACGCACACCATCAATAATAGACGCGTGGTTCAGTGCGTCGGAGATAATGGCGTCTTCCGCCCCCAGAAGCGTTTCAAACAGGCCACCGTTAGCATCAAAGCAGGAaga                 <  1:127309/250‑1 (MQ=255)                                                                                                                                                                                                                            ggCAATCCCGCCGTCAATTGAGAACACAGCATCGGTGGCGATCAGCACAGGACGCGCACCGGCTTCACGCGCTTCTTTCAGACGGGTTTCCAGCTCCTGCATATCGTTGTTGGCatagcgatagc                                                                                                                                             <  1:247280/125‑1 (MQ=255)                                                                                                                                                                                                                            ggCAATCACGCCGTCCATTGAGAACACAACATCGGTGGCGATCAGCACATAACGCGCACCAGCTTCACGAGCTTCTTTCAGACGTGCTTCCAGCTACTGCATATCGTTGTTGGCatagcgatagc                                                                                                                                             >  2:247280/1‑125 (MQ=255)                                                                                                                                                                                                                             gCAATCACGCCGTCCATTGAGAACACAGCATCGGTGGCGATCAGCACATGACGCGCACCGGCTTCACGCGCTTCTTTCAGACGTGCTTCCAGCTCCTGCATATCGTTGTTGGCATAGCGATAGCGTTTAGCTTTGCACAGACGCACACCATCAATAATAGACGCGTGGTTCAGTGCGTCGGAGATAATGGCGTCTTCCGCACCCAGAAGCGTTTCAAACAGGCCACCGTTAGCATCAAAGCAGGAAGAGTa              <  2:254345/251‑1 (MQ=255)                                                                                                                                                                                                                                  cACGCCGTCCATTGAGAACACAGCATCGGTGGCGATCAGCACATGACGCGCACCGGCTTCACGCGCTTCTTTCAGACGTGCTTCCAGCTCCTGCATATCGTTGTTGGCATAGCGATAGCGTTTAGCTTTGCACAGACGCACACCATCAATAATAGACGCGTGGTTCAGTGCGTCGGAGATAATGGCGTCTTCCGCACCCAGAAGCGTTTCAAACAGGCCACCGTTAGCATCAAAGCAGGAAGAGTAgaga          >  2:257177/1‑250 (MQ=255)                                                                                                                                                                                                                                     gCCGTCCATTGAGAACACAGCATCGGTGGCGATCAGCACATGACGCGCACCGGCTTCACGCGCTTCTTTCAGACGTGCTTCCAGCTCCTGCATATCGTTGTTGGCATAGCGATAGCGTTTAGCTTTGCACAGACGCACACCATCAATAATAGACGCGTGGTTCAGTGCGTCGGAGATAATGGCGTCTTCCGCACCCAGAAGCGTTTCAAACAGGCCACCGTTAGCATCAAAGCAGGAAGAGTAGAGAATCg      <  1:358669/251‑1 (MQ=255)                                                                                                                                                                                                                                        gTCCATTGAGAACACAGCATCGGTGGCGATCAGCACATGACGCGCACCGGCTTCACGCGCTTCTTTCAGACGTGCTTCCAGCTCCTGCATATCGTTGTTGGCATAGCGATAGCGTTTAGCTTTGCACAGACGCACACCATCAATAATAGACGCGTGGTTCAGTGCGTCGGAGATAATGGCGTCTTCCGCACCCAGAAGCGTTTCAAACAGGCCACCGTTAGCATCAAAGCAGGAAGAGTAGAGAATCGCAt   >  1:320758/1‑251 (MQ=255)                                                                                                                                                                                                                                         tCCATTGAGAACACAGCATCGGTGGCGATCAGCACATGACGCGCACCGGCTTCACGCGCTTCTTTCAGACGTGCTTCCAGCTCCTGCATATCGTTGTTGGCATAGCGATAGCGTTTAGCTTTGCACAGACGCACACCCTCAATAATAGACGCGTGGTTCAGTGCGTCGGAGATAATGGCGTCTTCCGCACCCAGAAGCGTTTCAAACAGGCCACCGTTAGCATCAAAGCAGGAAGAGAAGAGAATCGCATc  >  1:239417/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                            AGCCACTCAACCACTTCTTTGCGCGCCGCGGTATAACCACCAGAAGCCCCGCCCAGCGCTTTACCAAGCGTACCGGTGATAATATCGACCCGGCCCATCACATCGCAGTATTCATGGGAACCACGACCATTTTCACCGACAAAACCGACCGCGTGGGAGTCGTCTACCATCACCA‑GGGCATCATATTTATCTGCCAGATCGCAAACGCCCTTCAGGTTGGCAATCACGCCGTCCATTGAGAACACACCATCGGTGGCGATCAGCACATGACGCGCACCGGCTTCACGCGCTTCTTTCAGACGTGCTTCCAGCTCCTGCATATCGTTGTTGGCATAGCGATAGCGTTTAGCTTTGCACAGACGCACACCATCAATAATAGACGCGTGGTTCAGTGCGTCGGAGATAATGGCGTCTTCCGCACCCAGAAGCGTTTCAAACAGGCCACCGTTAGCATCAAAGCAGGAAGAGTAGAGAATCGCATC  >  NC_000913/3791761‑3792242

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: ‑