Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,547	C→T	V249V (GTC→GTT)	thrB →	homoserine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,547	0	C	T	94.1%	92.1 / ‑3.0	35	V249V (GTC→GTT)	thrB	homoserine kinase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (13/19);  total (14/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 = 6.41e-01

TTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTG‑CCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGG‑TATCTCCGGCTCCGGCCCGACC‑TTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCG  >  NC_000913/3302‑3781                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                             tttGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTG‑CCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCggg                                                                                                                                                                                                                                         >  1:258904/1‑250 (MQ=255)           ggCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTG‑CCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGcgcg                                                                                                                                                                                                                               <  1:280007/251‑1 (MQ=255)               gtgGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTG‑CCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGGAATCGGCGCGGt                                                                                                                                                                                                                             >  1:154166/1‑249 (MQ=255)                             aTCCGGGGATTAAAGGCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTG‑CCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGG‑TAtctc                                                                                                                                                                                                             <  1:35068/250‑1 (MQ=255)                                       tAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTG‑CCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATAGGCGCGGGAGCGAGCGG‑GATCCCCGGCTCCGGcc                                                                                                                                                                                                  >  1:207984/1‑251 (MQ=255)                                                         aGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTG‑CCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGG‑TATCTCCGGCTCCGGCCCGAAC‑TTGTTCGCTCtgt                                                                                                                                                                               >  1:144861/1‑251 (MQ=255)                                                                     tgtcccggtgcagttccgccgcCATGATTGAATTGCGCACGGGTGCCATCTGGCTGGCTTCATTCTCGCCTGCTATCCCCGTCAGCCTGAGCTTG‑CCGCGAAGCTGATGAAAGATGT‑TAGCGCTGATCCTTACCGTGAACGGTTACTGCCTGGTTTCCGGCAGGAGCGGCCGTCGGTTGCGGAAATCGGCGCGGTAGCGAGCGG‑TAGCTCCGGCTCCGGCCCGACC‑TTGTTCGCGATGTGTGACAAgacg                                                                                                                                                                    <  2:5604/234‑4 (MQ=255)                                                                        aCCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTG‑CCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGATAGCGAGCGG‑TATCTCCGGCTCCGGCCCGACC‑TTGTTCGCTCTGTGTGACAAGCCGGat                                                                                                                                                                 <  1:207840/250‑2 (MQ=255)                                                                             cgcAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTG‑CCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGAGGGCGG‑GATCTCCGGATCCGGCCCGACC‑TTGGTCGCTCTGGGTGTCAAGCCGGGAACCGcc                                                                                                                                                           >  1:265378/1‑251 (MQ=255)                                                                               cAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTG‑CCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGG‑GATCTCCGGCCCCGGCCCGACC‑TTGTTCGCTCTGTGTGACAAGCCGGATAACGcccc                                                                                                                                                         >  1:338857/1‑250 (MQ=255)                                                                                   aTCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTG‑CCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTTCGGAAATCGGCGCGGGAGCGAGCGG‑TATCTCCGGCTCCGGCCCGGCC‑TTGTTCGCTCTGTGTGACAAGCCGGAATACGCCCAgcgc                                                                                                                                                     >  1:196253/1‑251 (MQ=255)                                                                                   aTCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTG‑CCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGGAATAGGGGCGGGAACGAGAGG‑AATCTCCGGCTCCGGCCCGACC‑ATGTTCGCTCTGTGTCACAAGACGGATACCGCCCAgcgc                                                                                                                                                     >  1:261877/1‑251 (MQ=255)                                                                                      gccgccAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTG‑CCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGG‑TATCTCCGGCTCCGGCCCGACC‑TTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGtt                                                                                                                                                  <  1:178858/251‑1 (MQ=255)                                                                                       cctccAGGATTGCATTGCGCTCGGGCGACAGTTGGCTGGCTTCATTCACGCCGGCGATTCCCGTCAGCCTGAGCTTG‑CCGGGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCGGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCTGTAGCGAGCGG‑TATCTCCGGCTCCGTCCCGACC‑TTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTg                                                                                                                                                 <  2:368578/248‑1 (MQ=255)                                                                                        cgccAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCAGCTATTCCCGTCAGCCTGAGCTTG‑CCGCGAAGCTGATGAAAGATGT‑TCTCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGGGAACGGTGCGGTAGAGCGCGGGCATC‑CCGGCTCCGGGCCGACCATTGTT‑GGTCCGGGTGACAAGCCGGagcacgcccaacgccgtgg                                                                                                                                                >  2:77237/1‑233 (MQ=255)                                                                                               attgcattgcGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTG‑CCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGG‑TATCTCCGGCTCCCGCCCGACC‑TTGTTCGCTCTGTGTGACAAGCCGGGTACCGCCCAGCGCGTTGCCGGCTgg                                                                                                                                         >  1:122527/1‑251 (MQ=255)                                                                                                  gcatggcGCACGGGTGACCTCTGGCGGGCTTCAGTCACGCCTGCTTTTCCCGGCCGCCTGAG‑TTGTCCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGG‑TATCTCCGGCGCCGGCCCGACC‑TGGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGtt                                                                                                                                       <  2:338857/250‑1 (MQ=255)                                                                                                           aCGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTG‑CCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGG‑TATCTCCGGCTCCGGCCCGACC‑tt                                                                                                                                                                                          <  1:60431/190‑1 (MQ=255)                                                                                                                          ctggCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTG‑CCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGGAGCCAGCGG‑GATCTCCGGCTCCGGCCCGACC‑TTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATc                                                                                                                >  1:52712/3‑249 (MQ=255)                                                                                                                             gCTTCATGCACGCCTGATATTCTCGTCAGCCTGTGCTTG‑CCGAGTAGCTGATGAAAGA‑GTGTATCGTTGAACCCTCCCGGGAACGGTTAATGCCGGGCTTCCGGCAGGCGCTGCAGGCGGTGGCGGAAAACGGCGCGGGAGTGAGCGG‑TATCTCCGTCTCCGGCCCGACC‑TTGTTCGCTCTGTGTGACAAGCCGGATACCGCGCATCGCGTTGCCGACTGGTTGGGTAAGAACTTCCTGCAAATTCAGGaa                                                                                                           <  2:251038/251‑1 (MQ=255)                                                                                                                                     caCGCCTGCTATTCCCGTCAGCCTGAGCTTG‑CCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGTGGTTGCGGAAATCGGCGCGGTAGTGAGCGG‑TATCGCCTGCTCCGGCCCGACC‑TTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGt                                                                                                   <  1:48111/251‑1 (MQ=255)                                                                                                                                       cGCCTGCTATTCCCGTCAGCCTGAGCTTG‑CCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGGAATCGGCGCGGTAGCGAGCGG‑GATCTCCGGCTCCGGCCCGACC‑TTGGTCGCTCTGTGTGACAAGCCGGGTACCGCCCAGCGCGATGCCGACTGGGTGGGGAAGAACTACCTGCACAAACAGGAAGGTTTTGttt                                                                                                 >  1:137976/1‑250 (MQ=255)                                                                                                                                                 tCCCGTCAGCCTGAGCTTG‑CCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGG‑TATCTCCGGCTCCGGCCCGACC‑TTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGCCTGGGTGGGGAAGAACTACCTGCAAAATCAGGAAGGGTTTGTTCATATTTTCCg                                                                                       >  1:78528/1‑251 (MQ=255)                                                                                                                                                       cAGCCTGAGCTTG‑CCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGG‑TATCTCCGGCTCCGGCCCGACC‑TTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTgg                                                                                  <  1:110133/250‑1 (MQ=255)                                                                                                                                                             gtgcgtG‑CCGCGAAGCTGATGAAGGATGT‑TATCGCTGATCCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGG‑TATCTCCGGCTCCGGCCCGACC‑TTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATAcg                                                                             <  2:78528/244‑1 (MQ=255)                                                                                                                                                                       gcgAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGG‑TATCTCCGGCTCCGGCCCGACC‑TTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACggcggggggccg                                                                  >  1:167627/1‑245 (MQ=255)                                                                                                                                                                               gatgaAAGATGG‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCTGAAATCGGCGCGGTAGCGAGCGG‑TATCTCCGGCTCCGGCCCGACC‑TTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTgg                                                          <  2:167627/251‑1 (MQ=255)                                                                                                                                                                                  gaAAGATGT‑TATCGCTGAACACTACCGTGAACGGGTACTGCCTGGCTTCCGGCATGCGCGGCCGGCGGTTGCGGAAATCGGCTCGGTAGCGCTCGG‑TATCTCCGGCTCCGGCCCTACC‑TTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACGGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCGGGATACGGCGGGCGCACGAGTACTGGaaa                                                       <  2:9166/251‑1 (MQ=255)                                                                                                                                                                                   aaaGTTGT‑TATCGCTGAACCCTACCGTGAACGGGTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAAGCGGCGCGGTAGCGAGCGG‑TATCTCCGGCTCCGGCCCGACC‑TTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGaaaa                                                      <  1:12004/251‑1 (MQ=255)                                                                                                                                                                                                   gAACCCTACCGTGAACGTTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCTG‑TATCTCCGGCTCCGGCCCGACC‑TTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTTCTGGAAAACTAAATGAAACTCTa                                       <  1:330107/251‑1 (MQ=255)                                                                                                                                                                                                      cccTACCGTGAACTGTTACTTCCCGGCTTCCGGCATGCGCGGCATGCGGTTGCGTAAAGCGGCGCGGTAGCGAGACG‑TATCTCCGGCTCCGGCCCGACC‑TCGTTCGCTCTGTGTGACAAGCCGGATTCCGCCCAGCGCGTCGCCGACTTGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACaa                                    <  2:154166/251‑1 (MQ=255)                                                                                                                                                                                                        cTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGG‑TATCTCCGGCTCCGGCCCGACC‑TTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCaaaa                                                                                                                  <  1:150559/171‑1 (MQ=255)                                                                                                                                                                                                                          gcctggCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGG‑TATCTCCGGCTCCGGCCCGACC‑TTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAg                 <  1:18553/246‑1 (MQ=255)                                                                                                                                                                                                                              ggttcccGGCAGGCGCGGCAGGCGGTTGCTGAAATCGGCGCGGTGGCGAGCGG‑TATCTCCGGCGCCGGCCCGACC‑TTGTTCGCTCTGTGTGACAAGCCGGATACCGTCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAAGCAGGAAGGTTTTGGTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGTAAAACTAAATGAAACTCTACAATCTGAATGATCACAATGAGCAGGt            <  2:196253/246‑1 (MQ=255)                                                                                                                                                                                                                                        cggcgcggcagtagGTTGCTGAAATCGGCGCGGGAGCGATCGG‑TATCGCCGTCTCCGGCCCGACC‑TTGGTCGCTCTGTGTGCCAAGCCGGATACCGCCCATTGCGTTTCCGACTGGTTGGGTAAGAACTACCTTCAACATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTAATGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTgcg  <  2:302839/250‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                             TTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTG‑CCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGG‑TATCTCCGGCTCCGGCCCGACC‑TTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCG  >  NC_000913/3302‑3781

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