Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,520	A→T	P240P (CCA→CCT)	thrB →	homoserine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,520	0	A	T	94.3%	97.4 / NA	38	P240P (CCA→CCT)	thrB	homoserine kinase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  major base T (17/16);  minor base G (0/2);  total (18/20)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTG‑CCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGG‑TATCTCCGGCTCCGGCCCGACC‑TTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAA  >  NC_000913/3274‑3747                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                      cGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCTGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTG‑CCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCt                                                                                                                                                                                                                                   <  1:75637/251‑1 (MQ=255)   aCATCATCGGCCAGCTAGTGCCATGGTTTGATGATTGGCTGTGGGGGCTGGCGTTTCCGGGTATTAAAGTCTCGACGGTATAAGCCAGTGTTATTTTACCGGCGTAGTATCGCCGTCAGGATTGCCTTTCGCACGAGCGATATCTGGCTGGCTTCAGTCGCGCCTGCTATTCCCGTCCGCCTGAGCTTG‑CCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTc                                                                                                                                                                                                                                 <  2:237543/251‑1 (MQ=255)    catcatCAACCCGCAAGAGCCAGGGTTTTATGAGTTGCTGTGGGAGATGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGAGACATCTGTCAGGCTTCATTACCGCCTGCTATACCCGTAAGCCTGAGCTTG‑CCGCGAAGCTGATGCAAGATTT‑TATCGCTGAACCCTCCCGTAAACGGTTACTGCATGGCTTcc                                                                                                                                                                                                                                >  2:194377/1‑251 (MQ=255)           agccagcAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTG‑CCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCcggcaggc                                                                                                                                                                                                                         >  1:340357/1‑251 (MQ=255)                       ccAGGGTTTGATGAGTGGCTGTGGGTGCTGGAGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATATGGATTGCTTCATTACCGCATGCTATTCCCGTCAGCCTGAGCTTG‑CCGCGCAGCTGTTGAAAGATGT‑TATCGCTGAACCCTACCGTGAACAGTTACTGCCTCGCTTCcggcaggcagggaatgcggg                                                                                                                                                                                                             >  2:353408/1‑239 (MQ=255)                             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)                                              ggTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCTGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTG‑CCGCGACGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCTGGCTTCCGGCAGGCGCGgcgggcgggtggggggataggcgcggtagcgg                                                                                                                                                                                      >  1:370125/1‑227 (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)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                      CGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTGGCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCGCCGCCAGGATTGCATTGCGCACGGGCGACATCTGGCAGGCTTCATTCACGCCTGCTATTCCCGTCAGCCTGAGCTTG‑CCGCGAAGCTGATGAAAGATGT‑TATCGCTGAACCCTACCGTGAACGGTTACTGCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGG‑TATCTCCGGCTCCGGCCCGACC‑TTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAA  >  NC_000913/3274‑3747

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