Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,581,948	C→T	5.3%	Q367Q (CAG→CAA)	hsdM ←	DNA methyltransferase M

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,581,948	0	C	T	5.3%	234.7 / 2.6	75	Q367Q (CAG→CAA)	hsdM	DNA methyltransferase M
Reads supporting (aligned to +/- strand):  ref base C (33/37);  new base T (2/2);  total (35/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.70e-01

TGCGCGGGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCAGATGACACTTATCCATCAGGTCACGACGAATGTCGGTGCCTTTGCCGC  >  NC_000913/4581803‑4582096                                                                                                                                                  |                                                                                                                                                     tGCGCGGGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGa                                                                                                                                                   >  6:351029/1‑149 (MQ=255)      gggCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCgg                                                                                                                                              <  5:343907/149‑1 (MQ=255)      gggCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCgg                                                                                                                                              <  7:101574/149‑1 (MQ=255)       ggCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCTTGATTCggg                                                                                                                                             >  1:402634/1‑149 (MQ=255)          ttAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGAAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGGACACTTCGTATCCTGATTCGGGTTc                                                                                                                                          >  8:370141/1‑149 (MQ=255)          gtacaCCGTGCGGGTCTTCGCCATACACGCGCTCCAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCCGCATATTGGTACGCAGGTCAGACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTc                                                                                                                                          >  2:300103/5‑149 (MQ=255)            aaaCCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGTTAAACGGTGTGCGCTTTCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGc                                                                                                                                        >  6:42796/1‑149 (MQ=255)             aaCCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAATCTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGc                                                                                                                                        <  2:126513/148‑1 (MQ=255)             aaCCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCTGTACAGTTCTTATCCTGATTCGGGTTCGc                                                                                                                                        <  6:26414/148‑1 (MQ=255)                   gCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTc                                                                                                                                 >  1:62474/1‑149 (MQ=255)                       gTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCtt                                                                                                                             >  3:54885/1‑149 (MQ=255)                          ttCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTgg                                                                                                                          >  6:355770/1‑149 (MQ=255)                              ccATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTaaa                                                                                                                      <  2:30692/149‑1 (MQ=255)                                 tACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGaa                                                                                                                   >  1:4813/1‑149 (MQ=255)                                    acGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAg                                                                                                                >  6:178464/1‑149 (MQ=255)                                         cTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGt                                                                                                           >  7:156684/1‑149 (MQ=255)                                          tCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACCCCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGtt                                                                                                          >  7:323674/1‑149 (MQ=255)                                            aaaCGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTg                                                                                                         <  1:128093/148‑1 (MQ=255)                                            aaaCGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTg                                                                                                         <  4:145565/148‑1 (MQ=255)                                             aaCGGATGCAAATGCTCGTCGGTAAACGGTGTGCTCTTGCCGAAACTCGGCATATTGGTACGCATGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTAGGGTTCGCCACCGTACCTTTGGTAAAGAACAGCAAGTTgt                                                                                                        >  6:241078/1‑147 (MQ=255)                                              aCGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTc                                                                                                      <  8:27241/149‑1 (MQ=255)                                                ggCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCtt                                                                                                    >  2:378372/1‑149 (MQ=255)                                                  cTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCa                                                                                                  <  5:351029/149‑1 (MQ=255)                                                  cTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTACCTTTGGTAAAGAACAGCACGTTGGTCTTCa                                                                                                  >  3:198135/1‑149 (MQ=255)                                                     cAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGc                                                                                               >  5:187370/1‑149 (MQ=255)                                                        aTGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTgg                                                                                                                          >  6:316131/1‑119 (MQ=255)                                                        aTGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTgg                                                                                                                          <  5:316131/119‑1 (MQ=255)                                                         tGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCt                                                                                            >  6:339140/1‑148 (MQ=255)                                                              gtcgGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGTTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGc                                                                                        <  3:137210/147‑1 (MQ=255)                                                                     aaCGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTaaaaaa                                                                                <  8:267127/148‑1 (MQ=255)                                                                      aCGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATa                                                                              <  4:54885/149‑1 (MQ=255)                                                                              cgcTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCgg                                                                      <  7:140937/149‑1 (MQ=255)                                                                               gcTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCCCCGTCCCTTTGGTAAAGAAAAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCgg                                                                      >  1:204387/1‑148 (MQ=255)                                                                                    ccGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACg                                                                <  5:355770/149‑1 (MQ=255)                                                                                     cGAAACTCGGCATATTTGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCATGCCCTGAGCGTAAAAAATACCGGTCGGCAGACg                                                                >  6:347211/1‑148 (MQ=255)                                                                                      gAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACg                                                                >  1:400917/1‑147 (MQ=255)                                                                                         aCTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGaa                                                           <  5:204847/149‑1 (MQ=255)                                                                                         aCTCGGCATATTGGTACGAAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTTCTTCGGGTTCGCCTCCGTCACTTTGGTAAAGAACCGCACGTTGGTCTTGACGCAGTGAGCGTAATAAATACCGGTCGGCAGACGCAGaa                                                           <  5:241078/149‑1 (MQ=255)                                                                                             ggCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGgt                                                       <  1:378372/149‑1 (MQ=255)                                                                                             ggCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTAGCGACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGgt                                                       <  5:42796/149‑1 (MQ=255)                                                                                               cATATTTGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCATGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGgtgt                                                     <  5:347211/149‑1 (MQ=255)                                                                                                 tatTGGTACGCAGGACATACACCCCCACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGCACAGCACGTTGGTCTTCACGCCCTGAGCGTTAAAAATACCGGTCGGCAGACGCAGAATGGtgtg                                                    >  2:385338/1‑148 (MQ=255)                                                                                                      gTACGCAGGTCATACACCCACACATCATCGTTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCATACGCAGAATGGTGTGCAGAt                                               <  2:5428/148‑1 (MQ=255)                                                                                                      gTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCAGAt                                               >  7:271696/1‑148 (MQ=255)                                                                                                       tACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCTTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGGCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCAGATGa                                             <  2:402634/149‑1 (MQ=255)                                                                                                        aCGCAGGCCAAACACCAACACATCATCGGAAAAGTTCTTATCATGATTCGGTTTCGCCACCGTCCC‑TTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCAGATGaca                                           <  2:62474/149‑1 (MQ=255)                                                                                                         cGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCAGATGaca                                           <  2:400917/149‑1 (MQ=255)                                                                                                            aGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCAGATGACACtt                                        <  5:339140/149‑1 (MQ=255)                                                                                                             tgTCATACACCCACACATCATCGGTACAGGTCTTATCCTGATTCGGGTTCGCCACCGGCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCAGATGACACTTa                                       <  2:247023/148‑1 (MQ=255)                                                                                                              gTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCAGATGACACTTAt                                      >  1:189901/1‑149 (MQ=255)                                                                                                                 aTACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTCGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGAGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCAGATGACACTTATCCa                                   >  1:213317/1‑149 (MQ=255)                                                                                                                   acacCCACACATCATCGGTACAGTTCTTATCTTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCa                                                  >  5:197725/1‑132 (MQ=255)                                                                                                                   acacCCACACATCATCGGTACAGTTCTTATCTTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCa                                                  <  6:197725/132‑1 (MQ=255)                                                                                                                    cacCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGATTGGTGTGCAGATGACACTTATCCATCa                                <  8:357506/149‑1 (MQ=255)                                                                                                                     acCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCAGATGACACTTATCCATCAg                               >  2:221024/1‑149 (MQ=255)                                                                                                                       ccACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCAGATGACACTTATCCATCAGGt                             >  4:270101/1‑149 (MQ=255)                                                                                                                         acacaTCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCAGATGACACTTATCCATCAGGTCa                           >  2:259786/1‑149 (MQ=255)                                                                                                                            catcatCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCAGATGACACTTATCCATCAGGTCacga                        <  1:107922/149‑1 (MQ=255)                                                                                                                                atcGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCAGATGACACTTATCCATCAGGTCACGACGaa                    <  6:187370/149‑1 (MQ=255)                                                                                                                                atcGGTACAGTTCTTATCCTGATGCTGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCAGATGACACTTATCCATCAGGTCACGACGaa                    <  5:119440/149‑1 (MQ=255)                                                                                                                                 tcGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCAGATGACACTTATCCATCAGGTCACGACGAAt                   <  2:4813/149‑1 (MQ=255)                                                                                                                                    gTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCAGATGACACTTATCCATCAGGTCACGACGAATGTc                <  4:198135/149‑1 (MQ=255)                                                                                                                                     tACAGTTCTTATCCTGATTCGGGTTCGCCACCGTACCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCAGATGACACTTATCCATCAGGTCACGACGAATGTCg               >  3:394625/1‑149 (MQ=255)                                                                                                                                       cAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCAGATGACACTTATCCATCAGGTCACGACGAATGTCGGt             <  7:370141/149‑1 (MQ=255)                                                                                                                                       cAGTTCTTATCCTGATTCGGGTTCGCCACCGTACCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCAGATGACACTTATCCATCAGGTCACGACGAATGTCGGt             <  4:394625/149‑1 (MQ=255)                                                                                                                                         gTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCAGCTGACACTTATCCATCAGGTCACGACGAATGTCGGTGc           >  5:209647/1‑149 (MQ=255)                                                                                                                                         attCTTATCCGGATTCGGGGGCGCCACCGGCCCTTTGGTAAAGAACAGCACGTTGGGCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCCGACGCAGAATGGGGTGCAGATGACACTTATCCATCAGGTCACGACGAATGTCGGTGc           <  8:65059/148‑1 (MQ=255)                                                                                                                                          ttCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCAGATGACACTTATCCATCAGGTCACGACGAAt                   <  8:218363/140‑1 (MQ=255)                                                                                                                                          ttCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCAGATGACACTTATCCATCAGGTCACGACGAAt                   >  7:218363/1‑140 (MQ=255)                                                                                                                                          ttCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCAGATGACACTTATCCATCAGGTCACGACGAATGTCGGTGcc          >  2:29974/1‑149 (MQ=255)                                                                                                                                          ttCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCAGATGACACTTATCCATCAGGTCACGACGAATGTCGGTGcc          >  4:233846/1‑149 (MQ=255)                                                                                                                                           tCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCa                                                  >  7:21850/1‑108 (MQ=255)                                                                                                                                           tCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCa                                                  <  8:21850/108‑1 (MQ=255)                                                                                                                                            cTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCAGATGACACTTATCCATCAGGTCACGACGAATGTCGGTGCCt         <  8:437614/148‑1 (MQ=255)                                                                                                                                                tCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCAGATGACACTTATCCATCAGGTCACGACGAATGTCGGTGCCTTTgcc    <  1:221024/149‑1 (MQ=255)                                                                                                                                                 ccTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCAGATGACACTTATCCATCAGGTCACGACGAATGTc                >  7:64591/1‑136 (MQ=255)                                                                                                                                                 ccTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCAGATGACACTTATCCATCAGGTCACGACGAATGTc                <  8:64591/136‑1 (MQ=255)                                                                                                                                                  cTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGtgtg                                                    >  1:136271/1‑99 (MQ=255)                                                                                                                                                  cTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGtgtg                                                    <  2:136271/99‑1 (MQ=255)                                                                                                                                                  cTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCAGATGACACTTATCCATCAGGTCACGACGAATGTCGGTGc           <  5:7989/140‑1 (MQ=255)                                                                                                                                                  cTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCAGATGACACTTATCCATCAGGTCACGACGAATGTCGGTGc           >  6:7989/1‑140 (MQ=255)                                                                                                                                                  cTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCAGATGACACTTATCCATCAGGTCACGACGAATGTCGGTGCCTTTgccgc  >  3:359595/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                     TGCGCGGGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCTTGCCGAAACTCGGCATATTGGTACGCAGGTCATACACCCACACATCATCGGTACAGTTCTTATCCTGATTCGGGTTCGCCACCGTCCCTTTGGTAAAGAACAGCACGTTGGTCTTCACGCCCTGAGCGTAAAAAATACCGGTCGGCAGACGCAGAATGGTGTGCAGATGACACTTATCCATCAGGTCACGACGAATGTCGGTGCCTTTGCCGC  >  NC_000913/4581803‑4582096

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: ‑