Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,413,751	0	A	C	43.1%	84.3 / 25.4	58	N59K (AAT→AAG)	ralR	Rac prophage; restriction alleviation protein
Reads supporting (aligned to +/- strand):  ref base A (9/24);  new base C (25/0);  total (34/24)
Fisher's exact test for biased strand distribution p-value = 3.10e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGA  >  NC_000913/1413602‑1413901                                                                                                                                                       |                                                                                                                                                        gCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACA‑TGAACACCTCCa                                                                                                                                                         <  2:512789/150‑1 (MQ=255)    ctctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATA‑AGAAAACCTCCcttat                                                                                                                                                     >  2:69794/1‑151 (MQ=255)       tGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTATA‑CTAACACCACCattttttt                                                                                                                                                  >  1:351620/1‑149 (MQ=255)        ggAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACA‑TGAACACCTCCattattatt                                                                                                                                                 <  1:911878/151‑1 (MQ=255)          aaaCGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTc                                                                                                                                               <  2:286325/151‑1 (MQ=255)               ccTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTCA‑TTAAAACCTCCCTtttttttttcctt                                                                                                                                           >  2:611310/1‑145 (MQ=255)                          tttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTAt                                                                                                                                <  1:940010/150‑1 (MQ=255)                           ttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTTA‑TGAACACCCCCCCTATTTTTTTCACTGGtttggttttt                                                                                                                               >  2:957931/1‑142 (MQ=255)                           ttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAA‑TTAACACCTCCCTTATTATTTTCAGTGGGTCGtttttt                                                                                                                               >  1:654637/1‑150 (MQ=255)                           ttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTAtt                                                                                                                               <  2:222931/150‑1 (MQ=255)                             aCCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTGATTTTTTATA‑TGAAAACCTCCCCTATTTTTTTTCGTGGGTCGttttttt                                                                                                                              >  2:802922/1‑148 (MQ=255)                             aCCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTACA‑AGAACACCCCCCTTATTATTTTCAGTGGGTCGtttttttc                                                                                                                             >  1:425509/1‑148 (MQ=255)                                 gagaGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTAAA‑TGAAAACCCCCCTTATTTTTTTCAGTGGGTCGGTTTTTTCATCtt                                                                                                                        >  1:341364/1‑151 (MQ=255)                                    agCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTc                                                                                                                      <  1:462813/150‑1 (MQ=255)                                     gCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCTCCTTCCTTTTCAATAGTGGCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCgt                                                                                                                    <  2:714959/151‑2 (MQ=255)                                      cATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAg                                                                                                                   <  1:827609/151‑1 (MQ=255)                                       aTTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAA‑TTAAAACCCCCCTTATTTTTTTCAGTGGGTTGttttttttattttttttgt                                                                                                                  >  2:538736/1‑138 (MQ=255)                                         tttttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAA‑TTAAAACCCCCCTTTTTTTTTTCTGTGGTTTGTTTTTTTCATCTTTTGAGTg                                                                                                                 >  2:259769/1‑150 (MQ=255)                                            ttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGCCATTTTTTTCA‑TGAACACCTCCATTATTTTTTTCCGTGGGTCGTTTTTTTCATCTTTTTAGTGcttc                                                                                                             >  1:84109/1‑151 (MQ=255)                                             tCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGcttc                                                                                                             <  1:467497/150‑1 (MQ=255)                                                cATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATTAAAAACCCCCATT‑TTATTTTCCATGGGTTGTTTTTTTCATCTTTTTAGTGGTtttttt                                                                                                         >  1:607951/1‑151 (MQ=255)                                                   tCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATA‑TGAACAACTCCCTTTTTATTTTCAGTGGTTTGTTTTTTTTCTCTTTCTAGGGGTTCttttt                                                                                                        >  2:31077/1‑149 (MQ=255)                                                    cTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTTATTTTTTTTA‑TTAAAACCCCCCCTTTTTTTTTCATTGGGTTGTTTTTTTCATTTTTTCTGTGGttttttttt                                                                                                       >  2:599522/1‑148 (MQ=255)                                                      gATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACA‑TGAACACCCCCCTTTTTTTTTTCTATGGGTCGTTTTTTTCATCTTTTCAGTGCTTCtttttt                                                                                                       >  1:449273/1‑146 (MQ=255)                                                      gATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAc                                                                                                     <  1:775182/149‑1 (MQ=255)                                                        tttCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACA‑TTAACACCCCCCTTATTTTTTTCAGTGGGTCGTTTTTTTCATCTTTTTAGGGCttttttttttctt                                                                                                   >  1:335742/1‑144 (MQ=255)                                                          tCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATA‑TAAAAAACTCCCTTATTATTTTCATTGGGTTTTTTTTTTCATCTTTTTCGTGCTTCtttttttctttccc                                                                                               >  1:753480/1‑142 (MQ=255)                                                          tCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTcca                                                                                               <  2:503147/150‑1 (MQ=255)                                                              tAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAA‑TGAAAACCCCCCTTATTTTTTTTAGTGGGTCGGTTATTTCATTTTTTTAGTGCTTCtttttttctttc                                                                                                 >  1:746425/1‑138 (MQ=255)                                                                 tttAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATA‑TGAAAAACCCCCTTTTTTTTTTCAGTGGGTCGTTTTTTTCATCTTTTTAGTGCttttttttttctttcccc                                                                                              >  1:266075/1‑135 (MQ=255)                                                                       ttttGAATATCTTGTCCAGTTAGTAGGAGTGCCATCTTCCTTTTCAATAGTGGCGGTTATTTTTTAAA‑TTAACACCTCCCTTTTTTTTTTCAGTGGGTTGTTTTTTTTCTCTTTTTAGTGGTTTTTTTTTTCTTTCCCCCTTAACCgtt                                                                                    >  2:909941/1‑147 (MQ=255)                                                                             atatCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAA‑TTAAAAACCCCCTTTTTTTTTTCCGTGGGTCGGTTTTTTCCTCTTTTGAGGGGTTTTTTTTTTC‑TTTCCCCCTAACCGGGTTGGGGtt                                                                            >  1:302790/1‑150 (MQ=255)                                                                                   tGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCa                                                                        <  2:613908/148‑1 (MQ=255)                                                                                   tGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTAAA‑AGAACAACCCCCCTATTATTTTCAGTGGGTTGTTTTTTTCCTCTTTTGTGTGCTTTTTTTTTTC‑TTTCCCCCTAAACCGGTCgggggttgccttc                                                                     >  1:607359/1‑141 (MQ=255)                                                                                    gTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATcc                                                                     <  1:431189/150‑1 (MQ=255)                                                                                       cAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACA‑AGAAAACCCCCCCTTTTTTTTTCCGTGGGTTGTTTTTTTTCTCTTTTCTGTGCTTCTTTTTTTC‑TTCTCCCATAAACCGTTCcggggttgcctccctt                                                                  >  1:515545/1‑137 (MQ=255)                                                                                           tagtagGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACA‑TGAAAAACCCCCTTATTTTTTTCAGTGGGTTGTTTTTTTCCTCTTTTGAGTGCTTCTTTTTTTC‑TTTCCCCCTAAACGGGTTGGGGTTTGCATCCgcttca                                                               >  1:740402/1‑144 (MQ=255)                                                                                           tagtagGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTATAAA‑TGAAAACCCCCCCTTTTTTTTTTAGTGGTTTGTTTTTTTCCTTTTTTCTGTGCTTTTTTTTCAC‑TTTCCCCCAAAACGGTTCGGGATTCGCCTCCGtt                                                                  >  2:541931/1‑146 (MQ=255)                                                                                             gtagGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTc                                                           <  2:197024/151‑1 (MQ=255)                                                                                             gtagGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTc                                                           <  2:266852/151‑1 (MQ=255)                                                                                               agGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAA‑TGAAAACCCCCCTTTTTTTTTTCAGTGGGTCGTTTTTTTCTTCTTTCGAGTGCTTTTTTTTTTC‑TTCCACCATAACCCGGTTGGGGGTCGCATTCCTTaca                                                               >  1:720411/1‑145 (MQ=255)                                                                                               agGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCg                                                          <  2:73337/150‑1 (MQ=255)                                                                                                   gTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCgg                                                     <  1:469909/151‑1 (MQ=255)                                                                                                     gCCACCTTCCTTTTCAATAGTGGCGGCCATTTTTTTCA‑AGAACACCTCCCCTATTTTTTCTTGGGGGTCGGTTTTTTCATCTTTCTGGTGCTTTTTTTTTAT‑TTTCACCACAACCCGATTGGGGGTCCCATTCGTTaaac                                                              >  1:363811/1‑140 (MQ=255)                                                                                                         ccttccttTTCAATAGTGGCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAata                                                <  2:654566/150‑1 (MQ=255)                                                                                                                 ttCAATAGTGGCGGTAATTTTTTAAA‑TGAAAACCTCCCTTTTTTTTTTTAGTGGTTCGGTTTTTTCTTTTTTTTGGTGCTTTTTTTTTTC‑TTTCCCCATAAACGGTTTCGGGTTTGCATCCGTTTCACTTCCCCCGGTAatattccgcaa                                        >  2:957727/1‑143 (MQ=255)                                                                                                                  tCAATAGTGGCGGTAATTTTTTAAA‑TGAAAACCCCCCCTTTTTTTTTCAGTGGGTCGTTTTTTTCATTTTTTTAGTGCTTTTTTTTTTC‑TTTCACCATAACCCGTTTGGGGTTTGCCTCCGTTTCACTTCTCTCGGTAATATTCTGaaa                                        >  2:400767/1‑149 (MQ=255)                                                                                                                   cAATAGTGGCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATg                                      <  1:805760/150‑1 (MQ=255)                                                                                                                       aGTGGCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCtt                                  <  2:667769/150‑1 (MQ=255)                                                                                                                           gCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCAcc                             >  2:410655/1‑151 (MQ=255)                                                                                                                           gCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCAcc                             >  2:621973/1‑151 (MQ=255)                                                                                                                           gCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCAca                             <  2:495778/151‑2 (MQ=255)                                                                                                                             ggTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGt                           >  2:504674/1‑151 (MQ=255)                                                                                                                              gTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGtt                          >  1:445801/1‑151 (MQ=255)                                                                                                                              gTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGtt                          <  2:247666/151‑1 (MQ=255)                                                                                                                                aaTTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTAc                        <  1:640773/151‑1 (MQ=255)                                                                                                                                    ttATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGa                     <  2:503572/150‑1 (MQ=255)                                                                                                                                     tataCA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGAt                    <  1:370494/150‑1 (MQ=255)                                                                                                                                       taCA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATgg                  >  2:314548/1‑150 (MQ=255)                                                                                                                                               acacCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACcac          >  2:789176/1‑151 (MQ=255)                                                                                                                                                     ccATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAAt     <  2:535864/150‑1 (MQ=255)                                                                                                                                                     ccATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATg    >  2:64089/1‑151 (MQ=255)                                                                                                                                                       attattattTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATgg   >  1:443066/1‑150 (MQ=255)                                                                                                                                                       attattattTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGa  <  1:273926/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                        GCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACA‑TGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGA  >  NC_000913/1413602‑1413901

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