Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,950,292	C→G	10.3%	A77A (GCG→GCC)	aspS ←	aspartate‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,950,292	0	C	G	10.3%	261.0 / 18.9	87	A77A (GCG→GCC)	aspS	aspartate‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base C (28/50);  new base G (5/4);  total (33/54)
Fisher's exact test for biased strand distribution p-value = 2.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.77e-01

TAGCGGTATTTCAGACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGATCA  >  NC_000913/1950146‑1950439                                                                                                                                                   |                                                                                                                                                     tAGCGGTATTTCAGACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCa                                                                                                                                                   <  1:39320/150‑1 (MQ=255)   gCGGTATTTCAGACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑Cg                                                                                                                                                >  1:1029564/1‑150 (MQ=255)   gCGGTATTTCAGACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑Cg                                                                                                                                                >  2:1001153/1‑150 (MQ=255)      gTATTTCAGACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTAc                                                                                                                                             <  2:178666/150‑1 (MQ=255)       tATTTCAGACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCA‑ccgtac                                                                                                                                            >  1:9201/1‑145 (MQ=255)        aTTTCAGACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACgg                                                                                                                                           <  1:957562/150‑1 (MQ=255)        aTTTCAGACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACgg                                                                                                                                           >  1:1230163/1‑150 (MQ=255)            cAGACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGcc                                                                                                                                       <  2:351279/150‑1 (MQ=255)            cAGACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGcc                                                                                                                                       >  1:859454/1‑150 (MQ=255)             aGACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGccc                                                                                                                                      >  1:73796/1‑150 (MQ=255)              gACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCg                                                                                                                                     <  1:1311086/150‑1 (MQ=255)              gACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCg                                                                                                                                     <  2:98398/150‑1 (MQ=255)                cgcgCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTg                                                                                                                                   <  1:56722/150‑1 (MQ=255)                  cgcTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAc                                                                                                                                 <  2:9201/150‑1 (MQ=255)                  cgcTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAc                                                                                                                                 <  1:447102/150‑1 (MQ=255)                   gcTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAcg                                                                                                                                <  2:401147/150‑2 (MQ=255)                   gcTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGAcc                                                                                                                                <  2:939992/150‑1 (MQ=255)                     ttgttcGGTGTTGACGTGGTTAGAGTCTAGCGGCAGTACATGTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGGACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCTCA‑CGTACGGTGCCCGTGACCTg                                                                                                                              <  2:154588/147‑1 (MQ=255)                     ttcttcGGTGTTGACGTGGTTAGTGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTGAGCGAGGACGCAAGCACTTCGATTTAGCCGGTCGGCATATCGCGGTTAATATTTTTTTCGTCACGGGCACCGTAC‑GTGCGCGTGAcgaa                                                                                                                              <  2:1454647/150‑4 (MQ=255)                          cGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGc                                                                                                                         <  1:1164280/150‑1 (MQ=255)                            gtgtTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAg                                                                                                                       >  2:19760/1‑150 (MQ=255)                              gtTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGaa                                                                                                                     >  2:663544/1‑150 (MQ=255)                               ttGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGaat                                                                                                                    <  2:1029564/150‑2 (MQ=255)                               ttGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAAc                                                                                                                    <  1:668378/150‑1 (MQ=255)                                 gACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTc                                                                                                                  >  1:256426/1‑150 (MQ=255)                                 gACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTc                                                                                                                  <  1:1286577/150‑1 (MQ=255)                                 gACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTc                                                                                                                  >  2:57763/1‑150 (MQ=255)                                     tGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTa                                                                                                              <  2:73796/150‑1 (MQ=255)                                     tGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTa                                                                                                              <  1:895543/150‑1 (MQ=255)                                     tGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTa                                                                                                              <  2:1242636/150‑1 (MQ=255)                                     tGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTa                                                                                                              <  1:977977/150‑1 (MQ=255)                                     tGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTa                                                                                                              >  2:237860/1‑150 (MQ=255)                                       gTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCACGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACg                                                                                                            <  2:74934/150‑1 (MQ=255)                                        ttAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGc                                                                                                           <  1:382547/150‑1 (MQ=255)                                         tAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCa                                                                                                          >  1:734717/1‑150 (MQ=255)                                          agagTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAg                                                                                                         >  2:1013384/1‑150 (MQ=255)                                          agagTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAg                                                                                                         >  1:417544/1‑150 (MQ=255)                                                 aGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCagag                                                                                                  >  1:1522178/1‑150 (MQ=255)                                                 aGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCagag                                                                                                  >  1:1527875/1‑150 (MQ=255)                                                 aGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCagag                                                                                                  >  1:649277/1‑150 (MQ=255)                                                 aGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCagag                                                                                                  <  1:1092267/150‑1 (MQ=255)                                                  gCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAgg                                                                                                 <  2:948999/150‑1 (MQ=255)                                                  gCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAgg                                                                                                 <  2:1277875/150‑1 (MQ=255)                                                      cAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGCATGCAGAACTCATTACGCAGTTCAGAGGCCAg                                                                                             >  2:1613936/1‑150 (MQ=255)                                                       aGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGc                                                                                            <  2:1296881/150‑1 (MQ=255)                                                       aGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGc                                                                                            >  1:1319474/1‑150 (MQ=255)                                                          aCATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCttt                                                                                         <  2:1441403/150‑1 (MQ=255)                                                          aCATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCttt                                                                                         <  1:568906/150‑1 (MQ=255)                                                            aTCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTaa                                                                                       <  2:793738/150‑1 (MQ=255)                                                              cTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAcg                                                                                     <  1:19760/150‑1 (MQ=255)                                                              cTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAcg                                                                                     <  1:1515109/150‑1 (MQ=255)                                                                gcgcGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAAcgcg                                                                                   >  2:597730/1‑150 (MQ=255)                                                                 cgcgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGt                                                                                  >  1:83685/1‑150 (MQ=255)                                                                  gcgGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTc                                                                                 <  2:3457/150‑1 (MQ=255)                                                                    ggTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCg                                                                               <  1:165594/150‑1 (MQ=255)                                                                          tgatAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGAt                                                                         <  2:256426/150‑1 (MQ=255)                                                                             tAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCg                                                                      <  1:618401/150‑1 (MQ=255)                                                                               gTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCTCTCGATCCGGa                                                                    <  2:1319474/150‑1 (MQ=255)                                                                               gTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGa                                                                    >  1:513593/1‑150 (MQ=255)                                                                                 cAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATc                                                                  <  1:1041381/150‑1 (MQ=255)                                                                                 cAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATc                                                                  <  1:57763/150‑1 (MQ=255)                                                                                         aCGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATa                                                          <  1:605971/150‑1 (MQ=255)                                                                                         aCGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATa                                                          <  1:597730/150‑1 (MQ=255)                                                                                            ccAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCt                                                       >  1:1278057/1‑150 (MQ=255)                                                                                             cAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTg                                                      <  1:733131/150‑1 (MQ=255)                                                                                                cACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCAc                                                   <  1:1013384/150‑1 (MQ=255)                                                                                                  cTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGa                                                 >  1:1423588/1‑150 (MQ=255)                                                                                                    tCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATa                                               <  1:663544/150‑1 (MQ=255)                                                                                                     cGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATAc                                              >  2:1045820/1‑150 (MQ=255)                                                                                                            gCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTcgcg                                       <  1:78979/150‑1 (MQ=255)                                                                                                             ccGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCgg                                      >  1:1473030/1‑150 (MQ=255)                                                                                                                  cGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTcgcg                                 <  2:83685/150‑1 (MQ=255)                                                                                                                    ccATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCa                               <  2:734717/150‑1 (MQ=255)                                                                                                                     cATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCat                              >  1:1349051/1‑150 (MQ=255)                                                                                                                            cgGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATg                       <  2:513593/150‑1 (MQ=255)                                                                                                                            cgGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATg                       <  1:1001153/150‑1 (MQ=255)                                                                                                                              gTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGaa                     <  1:230237/150‑1 (MQ=255)                                                                                                                                    atTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAg               <  2:246487/150‑1 (MQ=255)                                                                                                                                    atTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAg               >  1:1178140/1‑150 (MQ=255)                                                                                                                                      ttttttCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGc             <  1:237860/150‑1 (MQ=255)                                                                                                                                      cattcaggTCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGc             >  2:37850/8‑150 (MQ=255)                                                                                                                                          ttCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTAcc         <  2:1522178/150‑1 (MQ=255)                                                                                                                                          ttCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTAcc         <  2:1527875/150‑1 (MQ=255)                                                                                                                                            ggTCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCaa       >  2:248905/2‑150 (MQ=255)                                                                                                                                             gTCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAg      >  1:1612344/1‑150 (MQ=255)                                                                                                                                             gTCACGGGCACCGTAC‑GTGCGCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAg      >  2:1423743/1‑150 (MQ=255)                                                                                                                                               cACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGAt    >  1:1562466/1‑150 (MQ=255)                                                                                                                                                 cgcgCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGATCa  >  1:1156685/1‑150 (MQ=255)                                                                                                                                                 cgcgCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGATCa  <  2:1230163/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     TAGCGGTATTTCAGACGCGCTTCTTCGGTGTTGACGTGGTTAGAGTCAAGCGGCAGAACATCTGCGCGGTTGATGATAGTCAGCGAGGACGCCAGCACTTCGATTTCGCCGGTCGCCATATCGCGGTTAATATTTTTTTCGTCACGCGCA‑CGTACGGTGCCCGTGACCTGAATGCAGAACTCATTACGCAGTTCAGAGGCCAGCTTTAACGCGTCCGCACGATCCGGATCGAAAAATACCTGCACGATACCTTCGCGGTCGCGCATATCGATGAAGATCAGGCTACCAAGATCA  >  NC_000913/1950146‑1950439

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