Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,367,700	G→T	5.2%	L210M (CTG→ATG)	aspA ←	aspartate ammonia‑lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,367,700	0	G	T	5.2%	317.1 / 3.3	96	L210M (CTG→ATG)	aspA	aspartate ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base G (46/44);  new base T (2/4);  total (48/48)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCATCTAC  >  NC_000913/4367553‑4367838                                                                                                                                                    |                                                                                                                                           aGCAAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTcatcag                                                                                                                                            >  2:161439/1‑148 (MQ=255)    cAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATg                                                                                                                                        >  2:370807/1‑149 (MQ=255)    cAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATg                                                                                                                                        >  8:33121/1‑149 (MQ=255)     aGTAACTTCAGACAGTTTTTTCACTGCCAGAGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTAACTTCTTCTTTCAGCAGGATGc                                                                                                                                       >  3:232139/1‑149 (MQ=255)      gTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCt                                                                                                                                      >  1:179086/1‑149 (MQ=255)      gTAACTTCAGCCAGTTTTTTCACTGCCAGAGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCt                                                                                                                                      >  1:69871/1‑149 (MQ=255)         aCTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGa                                                                                                                                    <  7:63032/148‑1 (MQ=255)          cTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGaaa                                                                                                                                  <  2:190871/149‑1 (MQ=255)             cAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCATCCTGCTGAAAg                                                                                                                                 <  6:95120/147‑1 (MQ=255)             cAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAgc                                                                                                                                <  7:33121/148‑1 (MQ=255)              aGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAgcgc                                                                                                                              <  3:315456/149‑1 (MQ=255)               gCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGcgcg                                                                                                                             <  1:132100/149‑1 (MQ=255)                ccAGTTTTTTCACTGCCAGATGAGAGTACTCTTTCGGCTTGTTCAGAACAGTACCGATTGCTGTGGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGagcgt                                                                                                                            <  4:224139/149‑2 (MQ=255)                   gTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAAt                                                                                                                         >  6:326621/1‑149 (MQ=255)                       attCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTcc                                                                                                                      <  1:29243/147‑1 (MQ=255)                        ttCACTGCCAGAGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCATGATGCTGAAAGCGCGGAATTCCt                                                                                                                     <  2:69871/148‑1 (MQ=255)                         tCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGa                                                                                                                   >  1:151497/1‑149 (MQ=255)                               ccAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAgg                                                                                                             >  5:298243/1‑149 (MQ=255)                               ccAGCGGAGAGGACTCTTTCGGCGTTGTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTAAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAgg                                                                                                             <  2:124821/149‑1 (MQ=255)                                cAGCGGATAGTACTCTTTCGGCGTGTTCAGTCCAGTACCTATTGATGTTGCTCCAAGGTTAACTTCCAGCAGCAGTGCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCATGATGCTGAAAGCGCGGAATTCCTTACCGAggg                                                                                                            <  4:232139/149‑1 (MQ=255)                                         gTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTa                                                                                                   <  8:18433/149‑1 (MQ=255)                                         gTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTa                                                                                                   <  4:69564/149‑1 (MQ=255)                                         gTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTa                                                                                                   >  4:307568/1‑149 (MQ=255)                                          tACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGa                                                                                                                   <  6:290554/132‑1 (MQ=255)                                          tACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGa                                                                                                                   >  5:290554/1‑132 (MQ=255)                                            ctctTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGAGTCATCGGTACt                                                                                                 <  2:151497/148‑1 (MQ=255)                                             tctTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTa                                                                                                   >  4:364003/1‑145 (MQ=255)                                             tctTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTa                                                                                                   <  3:364003/145‑1 (MQ=255)                                              ctTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTcagcag                                                                                                                                            >  6:227503/1‑103 (MQ=255)                                               tttCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGc                                                                                               <  2:178239/147‑1 (MQ=255)                                               tttCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGt                                                                                             >  5:195383/1‑149 (MQ=255)                                               tttCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGt                                                                                             >  1:53428/1‑149 (MQ=255)                                                 tCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACATTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCAGTACTGCGTc                                                                                            <  2:147622/148‑1 (MQ=255)                                                    gCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGAAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTg                                                                                         >  2:224565/1‑148 (MQ=255)                                                         ttCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCt                                                                                    <  6:195383/148‑1 (MQ=255)                                                            aGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCATCCTGCTGAAAgaag                                                                                                                              >  3:344012/1‑100 (MQ=255)                                                            aGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCATCCTGCTGAAAgaag                                                                                                                              <  4:344012/103‑4 (MQ=255)                                                            aGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGt                                                                                <  4:221064/149‑1 (MQ=255)                                                             gACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGATTTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTa                                                                               >  1:44253/1‑149 (MQ=255)                                                                cAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATAGGTACTGCGTCCTGCAGCTGGGTACgc                                                                            >  7:199820/1‑148 (MQ=255)                                                                      cGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCAtt                                                                      >  4:88498/1‑149 (MQ=255)                                                                      cGATTGCTGTTGCACCAAGGTTAACTTCCAGAAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGAGCGGAATTCCTGAACGAGGGTCATCGGTACTGCGTCCTGCAGCTTGGTACTACCCAtt                                                                      >  7:182809/1‑149 (MQ=255)                                                                        aTTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGATAAGCGGTACGTAGGATGTTTTTCACTTCTTCTATCAGCAGCATCCTGAAAGCGCCGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCAtttt                                                                    >  3:171639/1‑149 (MQ=255)                                                                         ttGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTc                                                                   <  1:370807/149‑1 (MQ=255)                                                                           gCTGTTCCACCAAGGTTAACTTCCAGCAGCAGTTCACCGGGACGTAGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGCCCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCa                                                                  <  2:41782/148‑1 (MQ=255)                                                                            cTGTTGCACCAATGTTAACTTCCAGGAGCATTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGAACGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAg                                                                 <  8:182809/148‑1 (MQ=255)                                                                            cTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTAc                                                                              <  2:370974/135‑1 (MQ=255)                                                                            cTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTAc                                                                              >  1:370974/1‑135 (MQ=255)                                                                            cTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAgg                                                                <  3:88498/149‑1 (MQ=255)                                                                               ttGCACCAAGGTTAACTTCGATCAGCAGTTCAGCGGTTCGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAgg                                                                <  1:161439/146‑1 (MQ=255)                                                                               ttGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATg                                                             <  1:390226/149‑1 (MQ=255)                                                                                   aCCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCt                                                         >  7:278412/1‑149 (MQ=255)                                                                                   aCCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCt                                                         <  2:344765/149‑1 (MQ=255)                                                                                    ccAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTg                                                        >  8:286070/1‑149 (MQ=255)                                                                                    ccAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTg                                                        >  1:344227/1‑149 (MQ=255)                                                                                      aaGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGCCCCATTTTCAGGATGTCCTgg                                                       >  4:371359/1‑148 (MQ=255)                                                                                       aGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGaa                                                     <  8:77365/149‑1 (MQ=255)                                                                                        ggTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAAt                                                    >  6:374373/1‑149 (MQ=255)                                                                                          ttAACTTCCAGCAGCAGTTCAGCGGTACGGTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACAGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGCCCCATTTTCAGGATGTCCTGGAATTc                                                  >  3:225701/1‑149 (MQ=255)                                                                                             aCTTCCAGCATCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGa                                                <  4:257416/148‑1 (MQ=255)                                                                                             aCTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGa                                                >  4:267068/1‑148 (MQ=255)                                                                                             aCTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGa                                                <  6:298243/148‑1 (MQ=255)                                                                                              cTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACt                                                                                                 <  4:37012/98‑1 (MQ=255)                                                                                              cTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACt                                                                                                 >  3:37012/1‑98 (MQ=255)                                                                                              cTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGa                                                >  7:338892/1‑147 (MQ=255)                                                                                              cTTCCAGCAGCAGGTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACa                                              <  2:50631/149‑1 (MQ=255)                                                                                               ttCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATTTCCTGTAATTCGACAg                                             >  4:30410/1‑149 (MQ=255)                                                                                               ttCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAg                                             >  5:293033/1‑149 (MQ=255)                                                                                               ttCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAg                                             >  6:81897/1‑149 (MQ=255)                                                                                                     catcagTTCATCGGTAAGTTTGATGTTTTTCAGTTCTTCTTTCTGCAGGATGATGTAAGCGCGGAATTCCTGTGCGAGGGTCTTCGGTACTGCGTCGTGCTGCTGGGTACGAACCATTTTCAGGATGTCCTGGAATTCGACAGGtttt                                        <  1:30710/145‑2 (MQ=255)                                                                                                     cagcagTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTAc                                       >  7:156153/1‑149 (MQ=255)                                                                                                      agcagTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACg                                      >  3:20229/1‑149 (MQ=255)                                                                                                       gcagTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGt                                     <  2:179086/149‑1 (MQ=255)                                                                                                         agTTCAGCGGTACGTTGGATGTTTTTCACATCTTCTTTAAGCAGGATGCTGAAAGCGCGGAAATCCTCACCGAGCGTCATCGGTACTGCGTCCTGCAGCTCGGTACCACCCATATTCAGGATGTCCTGGAATAAGACAGCTATACgat                                    >  4:215644/1‑146 (MQ=255)                                                                                                          gTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTc                                   >  3:36431/1‑148 (MQ=255)                                                                                                          gTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGCTCa                                  <  3:307568/149‑1 (MQ=255)                                                                                                             cAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCaaa                                <  6:96091/148‑1 (MQ=255)                                                                                                             cAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCaaa                                <  6:293033/148‑1 (MQ=255)                                                                                                              aGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGc                              <  5:374373/149‑1 (MQ=255)                                                                                                                  gTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTc                          <  6:371739/149‑1 (MQ=255)                                                                                                                  gTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTc                          <  4:20229/149‑1 (MQ=255)                                                                                                                     cGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACg                       <  8:156153/149‑1 (MQ=255)                                                                                                                         ggATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGt                   >  8:313338/1‑149 (MQ=255)                                                                                                                            tGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTgg                >  7:337281/1‑149 (MQ=255)                                                                                                                              tttttCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGt               >  4:275866/1‑148 (MQ=255)                                                                                                                               ttttCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTg                                                        <  3:84297/106‑1 (MQ=255)                                                                                                                               ttttCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTg                                                        >  4:84297/1‑106 (MQ=255)                                                                                                                               ttttCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTg                                                        >  4:89352/1‑106 (MQ=255)                                                                                                                               ttttCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTg                                                        <  3:89352/106‑1 (MQ=255)                                                                                                                                tttCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTa             >  2:322356/1‑148 (MQ=255)                                                                                                                                  tCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATc          >  6:247484/1‑149 (MQ=255)                                                                                                                                      ttcttcttTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCAt      >  7:32876/1‑149 (MQ=255)                                                                                                                                      ttcttcttTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCAt      >  7:295060/1‑149 (MQ=255)                                                                                                                                       cgcgcttTCAGCATCCTGCTGAAAGCGCGGAATTCCTGACCGAGGGTCCTCGGTACTGCGTCCTGCAGCTGGGTACGACCCCTTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCATc     >  3:266371/5‑149 (MQ=255)                                                                                                                                        cttcttTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCAt      <  8:278412/147‑1 (MQ=255)                                                                                                                                           cttTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGATGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCATCTAc  <  3:215644/148‑1 (MQ=255)                                                                                                                                           cttTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGTGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCGATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCATCTAc  <  4:171639/148‑1 (MQ=255)                                                                                                                                            tttCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCATCTAc  >  8:327550/1‑147 (MQ=255)                                                                                                                                                    |                                                                                                                                           AGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCATCTAC  >  NC_000913/4367553‑4367838

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