Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,889,662	A→G	E311E (GAA→GAG)	tnaA →	tryptophanase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,889,662	0	A	G	97.7%	306.2 / ‑2.4	87	E311E (GAA→GAG)	tnaA	tryptophanase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (40/45);  total (41/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

CCGCCAAGAAAGATGCGATGGTGCCGATGGGCGGCCTGCTGTGCATGAAAGACGACAGCTTCTTTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGG  >  NC_000913/3889529‑3889796                                                                                                                                      |                                                                                                                                       ccGCCAAGAAAGATGCGATGGTGCCGATGGGCGGCCTGCTGTGCATGAAAGACGACAGCTTCTTTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGggcggc                                                                                                                                  >  2:266194/1‑140 (MQ=255)    ccAAGAAAGATGCGATGGTGCCGATGGGCGGCCTGCTGTGCATGAAAGACGACAGCTTCTTTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGcgcg                                                                                                                               >  1:897181/1‑140 (MQ=255)    ccAAGAAAGATGCGATGGTGCCGATGGGCGGCCTGCTGTGCATGAAAGACGACAGCTTCTTTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGcgcg                                                                                                                               >  2:797403/1‑140 (MQ=255)         aaagaTGCGATGGTGCCGATGGGCGGCCTGCTGTGCATGAAAGACGACAGCTTCTTTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGa                                                                                                                          <  2:746705/140‑1 (MQ=255)          aagaTGCGATGGTGCCGATGGGCGGCCTGCTGTGCATGAAAGACGACAGCTTCTTTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAg                                                                                                                         <  1:916465/140‑1 (MQ=255)             aTGCGATGGTGCCGATGGGCGGCCTGCTGTGCATGAAAGACGACAGCTTCTTTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGt                                                                                                                      <  1:433771/140‑1 (MQ=255)               gCGATGGTGCCGATGGGCGGCCTGCTGTGCATGAAAGACGACAGCTTCTTTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCt                                                                                                                    <  1:797403/140‑1 (MQ=255)                 gATGGTGCCGATGGGCGGCCTGCTGTGCATGAAAGACGACAGCTTCTTTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGt                                                                                                                      >  1:610328/1‑136 (MQ=255)                 gATGGTGCCGATGGGCGGCCTGCTGTGCATGAAAGACGACAGCTTCTTTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGt                                                                                                                      <  2:610328/136‑1 (MQ=255)                   tGGTGCCGATGGGCGGCCTGCTGTGCATGAAAGACGACAGCTTCTTTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCg                                                                                                                <  2:146212/140‑1 (MQ=255)                   tGGTGCCGATGGGCGGCCTGCTGTGCATGAAAGACGACAGCTTCTTTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCg                                                                                                                >  2:211005/1‑140 (MQ=255)                    ggTGCCGATGGGCGGCCTGCTGTGCATGAAAGACGACAGCTTCTTTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCgg                                                                                                               <  2:427331/140‑1 (MQ=255)                    ggTGCCGATGGGCGGCCTGCTGTGCATGAAAGACGACAGCTTCTTTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCgg                                                                                                               <  1:753859/140‑1 (MQ=255)                    ggTGCCGATGGGCGGCCTGCTGTGCATGAAAGACGACAGCTTCTTTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCgg                                                                                                               <  2:897181/140‑1 (MQ=255)                      tGCCGATGGGCGGCCTGCTGTGCATGAAAGACGACAGCTTCTTTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTa                                                                                                             >  1:900568/1‑140 (MQ=255)                       gCCGATGGGCGGCCTGCTGTGCATGAAAGACGACAGCTTCTTTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAg                                                                                                            <  2:986744/140‑1 (MQ=255)                        ccGATGGGCGGCCTGCTGTGCATGAAAGACGACAGCTTCTTTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAgg                                                                                                           >  2:937230/1‑140 (MQ=255)                               gcggcCTGCTGTGCATGAAAGACGACAGCTTCTTTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTAt                                                                                                    >  1:120711/1‑140 (MQ=255)                                cggcCTGCTGTGCATGAAAGACGACAGCTTCTTTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATg                                                                                                   >  2:275775/1‑140 (MQ=255)                                cggcCTGCTGTGCATGAAAGACGACAGCTTCTTTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATg                                                                                                   <  2:724482/140‑1 (MQ=255)                                  gcCTGCTGTGCATGAAAGACGACAGCTTCTTTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGAc                                                                                                 >  2:593434/1‑140 (MQ=255)                                  gcCTGCTGTGCATGAAAGACGACAGCTTCTTTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGAc                                                                                                 >  2:860708/1‑140 (MQ=255)                                      gctgTGCATGAAAGACGACAGCTTCTTTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCa                                                                                             >  1:716016/1‑140 (MQ=255)                                      gctgTGCATGAAAGACGACAGCTTCTTTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCa                                                                                             <  1:101169/140‑1 (MQ=255)                                        tgtgCATGAAAGACGACAGCTTCTTTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATg                                                                                           >  2:565424/1‑140 (MQ=255)                                            cATGAAAGACGACAGCTTCTTTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAAtc                                                                                       >  1:67272/1‑140 (MQ=255)                                               gAAAGACGACAGCTTCTTTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCg                                                                                    >  2:169787/1‑140 (MQ=255)                                                aaaGACGACAGCTTCTTTGATGTGTACACCGAGTGCAGAGCCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGa                                                                                   <  1:211839/140‑1 (MQ=255)                                                aaaGACGACAGCTTCTTTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGa                                                                                   <  2:158771/140‑1 (MQ=255)                                                 aaGACGACAGCTTCTTTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGAc                                                                                  >  1:1043045/1‑140 (MQ=255)                                                  aGACGACAGCTTCTTTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGG‑AGGGGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGActg                                                                                >  2:504576/1‑140 (MQ=255)                                                   gacgacAGCTTCTTTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGActg                                                                                <  2:67272/140‑1 (MQ=255)                                                    acgacAGCTTCTTTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGActgg                                                                               >  2:214478/1‑140 (MQ=255)                                                         aGCTTCTTTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACtggctggc                                                                          <  1:721495/140‑1 (MQ=255)                                                            ttcttTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTta                                                                       <  1:290644/140‑1 (MQ=255)                                                                  gATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTtatcgtat                                                                 <  2:891736/140‑1 (MQ=255)                                                                  gATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTtatcgtat                                                                 <  2:716115/140‑1 (MQ=255)                                                                  gATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTtatcgtat                                                                 <  2:716016/140‑1 (MQ=255)                                                                  gATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTtatcgtat                                                                 <  2:900568/140‑1 (MQ=255)                                                                    tgtgTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTtatcgtatcg                                                               <  1:214478/140‑1 (MQ=255)                                                                      tgtACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATcgcg                                                             <  2:120711/140‑1 (MQ=255)                                                                      tgtACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATcgcg                                                             <  1:266194/140‑1 (MQ=255)                                                                       gtACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCgcgc                                                            <  2:899795/140‑1 (MQ=255)                                                                       gtACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCgcgc                                                            <  2:1043045/140‑1 (MQ=255)                                                                       gtACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCgcgc                                                            <  2:225515/140‑1 (MQ=255)                                                                       gtACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCgcgc                                                            <  1:860708/140‑1 (MQ=255)                                                                        tACACCGAGTGAAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCa                                                           >  1:598137/1‑140 (MQ=255)                                                                         acacCGAGTGAAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAg                                                          <  2:598137/140‑1 (MQ=255)                                                                           acCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGt                                                        <  1:26591/140‑1 (MQ=255)                                                                            ccGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTa                                                       <  1:593434/140‑1 (MQ=255)                                                                            ccGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTa                                                       <  1:275775/140‑1 (MQ=255)                                                                             cGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTta                                                                       >  1:89801/1‑123 (MQ=255)                                                                             cGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTta                                                                       <  2:89801/123‑1 (MQ=255)                                                                             cGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATcgcg                                                             >  1:137286/1‑133 (MQ=255)                                                                             cGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATcgcg                                                             <  2:137286/133‑1 (MQ=255)                                                                              gAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTa                                                       >  2:238361/1‑138 (MQ=255)                                                                              gAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTa                                                       <  1:238361/138‑1 (MQ=255)                                                                               aGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCa                                                           <  1:980015/133‑1 (MQ=255)                                                                               aGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCa                                                           >  2:980015/1‑133 (MQ=255)                                                                               aGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAg                                                    >  2:820268/1‑140 (MQ=255)                                                                               aGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAg                                                    >  2:744369/1‑140 (MQ=255)                                                                               aGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAg                                                    >  1:425640/1‑140 (MQ=255)                                                                               aGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAg                                                    >  2:956646/1‑140 (MQ=255)                                                                                gTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGt                                                   >  1:893771/1‑140 (MQ=255)                                                                                  gCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTAt                                                 >  2:69955/1‑140 (MQ=255)                                                                                        cccTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTc                                           >  1:570234/1‑140 (MQ=255)                                                                                         ccTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCg                                          <  1:956646/140‑1 (MQ=255)                                                                                          cTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGa                                         <  1:172628/140‑1 (MQ=255)                                                                                          cTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGa                                         <  1:820268/140‑1 (MQ=255)                                                                                          cTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGa                                         <  1:69955/140‑1 (MQ=255)                                                                                           tttGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGAt                                        >  1:220690/1‑140 (MQ=255)                                                                                               cGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTc                                    <  1:937230/140‑1 (MQ=255)                                                                                                           aGGCTTCCCGACATATGGCGGCCTGG‑AGGGGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTgg                       <  1:504576/140‑1 (MQ=255)                                                                                                              cTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATg                                       >  2:56673/1‑122 (MQ=255)                                                                                                              cTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATg                                       <  1:56673/122‑1 (MQ=255)                                                                                                              cTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCg                     >  2:522601/1‑140 (MQ=255)                                                                                                              cTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCg                     >  2:103373/1‑140 (MQ=255)                                                                                                              cTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCg                     >  1:587389/1‑140 (MQ=255)                                                                                                              cTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCg                     >  2:368960/1‑140 (MQ=255)                                                                                                              cTTCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCg                     >  2:491121/1‑140 (MQ=255)                                                                                                              cTTCCCGACATATGGCAGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCg                     >  2:304199/1‑140 (MQ=255)                                                                                                                tCCCGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGtt                   <  2:425640/140‑1 (MQ=255)                                                                                                                 cccGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTg                  >  1:490603/1‑140 (MQ=255)                                                                                                                  ccGACATATGGCGGCCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGt                 >  2:366846/1‑140 (MQ=255)                                                                                                                             cggcCTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAgg      <  1:882572/140‑1 (MQ=255)                                                                                                                                ccTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAggcgg   >  2:810790/1‑140 (MQ=255)                                                                                                                                 cTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAggcggg  >  1:704883/1‑140 (MQ=255)                                                                                                                                 cTGGAGGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAggcggg  <  1:810790/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                                       CCGCCAAGAAAGATGCGATGGTGCCGATGGGCGGCCTGCTGTGCATGAAAGACGACAGCTTCTTTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGG  >  NC_000913/3889529‑3889796

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑