Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,102	T→C	G123G (GGT→GGC)	thrC →	threonine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,102	0	T	C	97.7%	125.4 / ‑6.5	43	G123G (GGT→GGC)	thrC	threonine synthase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (25/17);  total (26/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

TGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCACAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTG‑GCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGTGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCG‑AATGTGAAAGTGGTTATCCTCTATCCAC‑GAGGCAAAATCAGTCCACTGCAAGAAAAAC‑TGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAACTGAA  >  NC_000913/3873‑4341                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   tGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTG‑GCATTTAAAGATTTCGGCGTTCGCTTTCTAGACCCAATGCTGACCCATATTGCGGGCGAtgatgcagtggtcattctt                                                                                                                                                                                                                               >  2:363197/1‑233 (MQ=255)            tttGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTG‑GCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCt                                                                                                                                                                                                                     <  1:165296/250‑1 (MQ=255)              gggcaaccGCAGTGCGAAGATCCTCGCGGCTGTTATTGGTGATGCAATCCCGCATGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGGTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTG‑GCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTcg                                                                                                                                                                                                                   <  2:114657/244‑2 (MQ=255)                atccccgcagtgatgatatcttcgtggCGTTTATTGGTTATGAACTCCCGCAGGAAAGCCTGGAAGAGCTCGGTCGCTGGGCGGTAGCGTTACCGGCTCTGGGCGCCAATGTTGAAAGCGAGGTCGGGTGTCGGGGATTGTTCCAGGGGCCAACGCTG‑GCATTTAAAGATTTCGGCGGTCGCTTTATGGCCCAAATGCTGACCCAAATTTGGGGCGATAAGCAAGTGATCATTCTGATCGCGACATCCGGt                                                                                                                                                                                                                <  2:202828/226‑1 (MQ=255)                   cccGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTG‑GCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGAACCATATTGCGGGCGAGAAGCCCGTGACCATTCTGACCGCGACCTCCGGGGAt                                                                                                                                                                                                             >  2:393458/1‑251 (MQ=255)                     cGCAGTGCGAAGATACTCTCGGCGTTTATTGGTGATAAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGACTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTG‑GCATTTAAAGATTTCGGCGGACGCTTTATTGCACAAATGCTGACCCATATTTCGGGCGATTAGCCAGTGACACTACTGACCGCGACCTCCGGTgctcc                                                                                                                                                                                                           >  1:293726/1‑247 (MQ=255)                                         ggCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTG‑GCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTC‑GGTGATACCGGAGCGGCAGTGGCTCATGc                                                                                                                                                                                      <  1:74165/251‑1 (MQ=255)                                           ggTTTATTGGTGAGGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTG‑GCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCt                                                                                                                                                                                     <  2:9145/250‑1 (MQ=255)                                                   ggTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGGGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTG‑GCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACgg                                                                                                                                                                             <  2:282959/251‑1 (MQ=255)                                                       atgaAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTG‑GCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAACCAGTGACCATTCTGACCGCAACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGtttt                                                                                                                                                                         >  1:415541/1‑250 (MQ=255)                                                                        aTCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTG‑GCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGGGGCTCATGCTTTCTACGGTTTACCG‑AATGGGAAAGTggg                                                                                                                                                       >  1:225938/1‑250 (MQ=255)                                                                            tGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTG‑GCATTTAAAGATTTCGGCGGTCGCTTTATGGCCCAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCG‑AATGTGAAAGTGGttttc                                                                                                                                                   >  1:103804/1‑248 (MQ=255)                                                                                 gagCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTG‑GCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGGGCGGCAGTGGCTCATGCTTTCTACGGTTTACCG‑AATGTGAAAGGGGTTATCctct                                                                                                                                               >  1:441470/1‑250 (MQ=255)                                                                                            cggggcgTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTTTCGGGAATTGTTCCACGGGCCAACGCTG‑GCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCG‑AATGTGAAAGTGGTTATCCTCTATCCAC‑GAGGCa                                                                                                                                  <  2:197126/248‑1 (MQ=255)                                                                                                  gTTTGCCTTCCCGGCTCCGGTCGCCAATGGTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTG‑GCATTTAAAGAGTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTTACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCG‑AATGTGAAAGTGGTTATCCTCTATCCAC‑GAGGCAAAATCa                                                                                                                            <  2:433245/251‑1 (MQ=255)                                                                                                     tGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTG‑GCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCG‑AATGTGAAAGTGGTTATCCTCTATCCAC‑GAGGCAAAATCAGt                                                                                                                          <  1:358519/250‑1 (MQ=255)                                                                                                      gCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTG‑GCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCG‑AATGTGAAAGTGGTTATCCTCTATCCAC‑GAGGCAAAATCAGTcc                                                                                                                        <  2:317127/251‑1 (MQ=255)                                                                                                            ccGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTG‑GCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGAACCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCCGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCG‑AATGTGAAAGTGGTTATCCTCTATCCAC‑GAGGCAAAATCAGTCCCCTGc                                                                                                                   >  1:67633/1‑250 (MQ=255)                                                                                                                cTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTG‑GCATTTAAAGATTTCGGCGGTCGCTTTGTGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGGGCGGCAGTGGCTCATGCTTTCTACGGTTTACCG‑AATGTGGAAGGGGTTATCCTCTATCCAC‑GAGGCCAAATCAGTCCACTGCAAGa                                                                                                               >  1:69713/1‑250 (MQ=255)                                                                                                                 tCCGGTCGCAAATGTTTAAAGCGATGGCGGTTGTATGGAATTGTTCCACGGGCCAACGCTG‑GCATTTAAAGATTTCGGCGGTTGCTTTATGGCACAAATGCTGGCCCATATTGCGGCCGATAAGCCAGTGACCATTCTGACCCCGACCTCCGGTAGTACCTGAACGGAAGTGGCTGATGCTTTTTACGGTTTACCG‑AATGTGGAAGTGGTGATCCCTTATCCAC‑TGGACACAATTAGTCCACTgaaaaaa                                                                                                              >  2:414253/1‑244 (MQ=255)                                                                                                                       cGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTG‑GCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCG‑GATGTGAAAGTGGTTATCCTCTATCCAC‑GAGGCAAAATCAGTCCACTGCAAGAAAAAC‑TGt                                                                                                      >  1:474000/1‑251 (MQ=255)                                                                                                                                          gTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTG‑GCATTTAAAGATTTCGGCGGTCGCTTTTTGGCACAAATGCTGGCCCATATTGCGGGCGATAAGCCAGTTACCATTCTGACCGCCACCTCCGGTGATAACGGAGCGGCAGGTGCTCATTCCTTTTACTGTTTACCGGATTTTGAAGGTG‑TTATCCTCTTTCCCC‑GAGGCAAAATCAGTACCCTTCAAGAAAAAC‑ATTTTTTTACATTGGGCGGCa                                                                                    >  2:111927/1‑250 (MQ=255)                                                                                                                                            cggtggttgtgAATTGTTCCACGTGCCACCG‑TGTGCATTTAAAGAGTTCGGCGGGCGCTTTTTGGCACAAACTCTGACCCATATTGCGGGCGTTAAGCCAGTGACCATTGTGACCGCGACCTCCGGTGAGACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCT‑AATGTGAAAGTGGTTAGCCTCTATCCAC‑GAGGCAAAAACAGTCCACTGCAAGAAAAAC‑TGTTCTGTACATTGGGCGGCAat                                                                                  <  2:273976/240‑1 (MQ=255)                                                                                                                                               ttGTCTGGAATTGTTCCACGGGCCAACGCTG‑GCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCG‑AATGTGAAAGTGGTTATCCTCTATCCAC‑GAGGCAAAATCAGTCCACTGCAAGAAAAAC‑TGTTCTGTACATTGGGCGGCAATATCg                                                                              >  1:483733/1‑251 (MQ=255)                                                                                                                                                   cTGGAATTGTTCCACGGGCCAACGCTG‑GCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCG‑AATGTGAAAGTGGTTATCCTCTATCCCC‑GAGGCAAAATCAGTCCACTGCAAGAAAAAc                                                                                                          >  1:56662/1‑220 (MQ=255)                                                                                                                                                   cTGGAATTGTTCCACGGGCCAACGCTG‑GCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCG‑AATGTGAAAGTGGTTATCCTCTATCCAC‑GAGGCAAAATCAGTCCACTGCAAGAAAAAc                                                                                                          <  2:56662/220‑1 (MQ=255)                                                                                                                                                            gtCCACGGGCCAACGCTG‑GCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCG‑AATGTGAAAGTGGTTATCCTCTATCCAC‑GAGGCAAAATCAGTCCACTGCAAGAAAAAC‑TGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCAt                                                                 <  2:103804/250‑1 (MQ=255)                                                                                                                                                                      cAACGCTG‑GCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGGCCCATATTGCGGGCGGTAAGCCCGTGACCATTCTGAACGCGACCTCCGGTGATGCCGGAGCCGCAGGGGCCCCTGCTTTCTACGGTTTACCG‑AATGGGAAAAGGGGTATCCCCTATCCACCGAGCCAAAA‑CACTCCACTGCCAGAAAAACCAGTT‑TTTAAATTTGGCGGCAATATCCGAACTTTTGCCATCCACCGCGAt                                                       >  2:174879/1‑251 (MQ=255)                                                                                                                                                                         cGCTG‑GCATTTAAAGATTTCGGCGGTCGCTTTATGGCGCAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCG‑AATGTGAAAGTGGTTATCCTCTATCCAC‑GAGGCAAAATCAGTCCACTGCAAGAAAAAC‑TGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTc                                                    <  1:386420/251‑1 (MQ=255)                                                                                                                                                                                   ttAAAGATTTCGGCGGTTGCTTTATGGCACAAATGCTGGCCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCG‑AATGTGAAAGTGGTTATCCTCTATCCAC‑GAGGCAAAAGCAGTCCACTGCAAGAAAAAC‑TGTTCTGTACAGTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTGCGATGCCTg                                            <  2:69713/250‑1 (MQ=255)                                                                                                                                                                                       aGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCG‑AATGTGAAAGTGGTTATCCTCTATCCAC‑GAGGCAAAATCAGTCCACTGCAAGAAAAAC‑TGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAgg                                       >  1:65541/1‑251 (MQ=255)                                                                                                                                                                                           ttCGGCGGTCGCTGTATGGCCCAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCCTTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGGGGCTCATTCTTTTTAAGGGTTACCG‑AATGTGGAAGGGGGTATCCTCTATCCAC‑GGGGAAAAATCAGTCCACTGCAAGAAAAAC‑TGTTCTGGACATTGGGCGGAAATCTCGAAACTGTTGCCATCGACGGCGAGTTCGATGCCTGGCCGgcgcc                                   >  2:49830/1‑250 (MQ=255)                                                                                                                                                                                            tCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCG‑AATGTGAAAGTGGTTATCCTCTATCCAC‑GAGGCAAAATCAGTCCACTGCAAGAAAAAC‑TGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCt                                   <  2:65541/250‑1 (MQ=255)                                                                                                                                                                                                 ggTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCG‑AATGTGAACGTGGTTATCCTCTATCCAC‑GAGGCAAAATCAGTCCACTGCAAGAAAAAC‑TGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGa                              >  1:28287/1‑250 (MQ=255)                                                                                                                                                                                                     gCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCG‑AATGTGAAAGTGGTTATCCTCTATCCAC‑GAGGCAAAATCAGTCCACTGCAAGAAAAAC‑TGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGc                           >  1:115643/1‑249 (MQ=255)                                                                                                                                                                                                      cTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGAGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCG‑AATGTGAAAGTGGTTATCCTCTATCCAC‑GAGGCAAAATCAGTCCAATGCAAGAAAAAC‑TGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGcggg                        >  2:463741/1‑248 (MQ=255)                                                                                                                                                                                                             gCACAAATGCTGACCCATATTGCGGGCGATAAGCAAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCG‑AATGTGAAAGTGGTTCTCCTCTATCCAC‑GAGGCAAAATCAGTCCACTGCAAGAAAAAC‑TGCTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGCCGGCGATTTCGATGCCTGTCAGGCGATGGTGAAGCAGGCGTTtg                  >  2:315519/1‑250 (MQ=255)                                                                                                                                                                                                               acaAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTCCGGTTTACCG‑AATGTGAAAGTGGTTATCCTCTATCCAC‑GAGGCAAAATCAGTCCACTGCAAGAAAAAC‑TGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTTAAGCAGGCGTTtggtg               >  2:234917/1‑248 (MQ=255)                                                                                                                                                                                                                 aaaTGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCG‑AATGTGAAAGTGGTTATCCTCTATCCAC‑GAGGCAAAATCAGTCCACTGCAAGAAAAAC‑TGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAGGCAGGCGTTtgatgat             <  2:483733/251‑1 (MQ=255)                                                                                                                                                                                                                   aTGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCG‑AATGTGAAAGTGGTTATCCTCTATCCAC‑GAGGCAAAATCAGTCCACTGCAAGAAAAAC‑TGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGGTtgatgatgg           >  2:393032/1‑250 (MQ=255)                                                                                                                                                                                                                       tGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCG‑AATGTGAAAGTGGTTATCCTCTATCCAC‑GAGGCAAAATCAGTCCACTGCAAGAAAAAC‑TGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGCTTTCGATGCCTGTCAGGCGCTGGTGACGCAGGCGTTTGTTGATgaag         >  2:205177/1‑249 (MQ=255)                                                                                                                                                                                                                         aCCCATATTGCGGGTGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCG‑AATGTGAAAGTGGTTATCCTCTATCCAC‑GAGGCAAACTCAGTCCACTGCAAGAAAAAC‑TGTTCTGTACATTGGGCGGCAATATCGAAAATGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGGTGAAGAACt     >  2:222595/1‑251 (MQ=255)                                                                                                                                                                                                                             atatTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCG‑AATGTGAAAGTGGTTATCCTCTATCCAC‑GAGGCAAAATCAGTCCACTGCAAGAAAAAC‑TGTTCTGTACATTGGGCGGCAATATCGCAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCTGTTGATGATGAAGAGCTGaa  >  2:378608/1‑250 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   TGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCACAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTG‑GCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGTGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCG‑AATGTGAAAGTGGTTATCCTCTATCCAC‑GAGGCAAAATCAGTCCACTGCAAGAAAAAC‑TGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAACTGAA  >  NC_000913/3873‑4341

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