Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,413,735	0	A	T	46.6%	60.7 / 40.9	58	Y2* (TAT→TAA)	ydaC	DUF1187 family protein, Rac prophage; putative double‑strand break reduction protein
Reads supporting (aligned to +/- strand):  ref base A (7/24);  new base T (27/0);  total (34/24)
Fisher's exact test for biased strand distribution p-value = 2.48e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACA‑CCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TT‑CCACCA‑TAA‑CCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTG  >  NC_000913/1413595‑1413881                                                                                                                                             |                                                                                                                                                         aCTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAac                                                                                                                                                <  1:1223155/150‑1 (MQ=255)  cTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATAT‑Taaaa                                                                                                                                               >  2:452350/1‑148 (MQ=255)    ttCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAAT‑Gaaaa                                                                                                                                               >  1:817208/1‑146 (MQ=255)       tGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGAGAGTTTATTCaacgataaccccc                                                                                                                                          >  2:483167/1‑139 (MQ=255)           ctctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACA‑CCTCCatt                                                                                                                                      <  2:595679/149‑1 (MQ=255)            tctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCCCCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACA‑CCTCCatt                                                                                                                                      <  1:662798/148‑1 (MQ=255)            tctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAAT‑GAACA‑CCCCCcttttt                                                                                                                                   >  2:771535/1‑151 (MQ=255)             ctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATAT‑GAACA‑CCTCCcttttt                                                                                                                                   >  1:1342532/1‑150 (MQ=255)              tGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATAT‑TAacaacccccttttt                                                                                                                                    >  1:148467/1‑138 (MQ=255)                gAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTTAA‑TAACA‑ACTCCCTTAtttttt                                                                                                                               >  2:383968/1‑151 (MQ=255)                   aCGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACA‑CCTCCATTATTATTTCCa                                                                                                                            <  1:964053/151‑1 (MQ=255)                      ccTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACA‑CCTCCATTATTATTTCCAGTg                                                                                                                         <  2:208579/151‑1 (MQ=255)                               tttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACA‑CCTCCATTATTATTTCCAGTGGTTCGtt                                                                                                                  <  1:880047/149‑1 (MQ=255)                               tttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACA‑CCTCCATTATTATTTCCAGTGGTTCGtt                                                                                                                  <  2:1340019/149‑1 (MQ=255)                               tttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACA‑CCTCCATTATTATTTCCAGTGGTTCGt                                                                                                                   <  2:986998/148‑1 (MQ=255)                                  ttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAAT‑GAAAA‑ACCCCCTTATTTTTTTCAAGGGTTTGttttttt                                                                                                             >  2:156838/1‑150 (MQ=255)                                        gagaGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATAT‑TAAAA‑CCTCCCTTTTTTTTTTCAGTGGTTCGTTTTTTTCATCtt                                                                                                       >  1:461684/1‑151 (MQ=255)                                        gagaGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGAAAATTTTTTTAC‑GAAAA‑ACTCCCCTTTTATTTTTTGAGGTTAGGTTATTTCATCt                                                                                                        >  2:777359/1‑150 (MQ=255)                                              aTTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACA‑CCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGt                                                                                                 <  2:530756/151‑1 (MQ=255)                                              aTTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACA‑CCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGt                                                                                                 <  2:959289/151‑1 (MQ=255)                                                   ttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATAAAT‑GAAAA‑ACCCCCTTTTTTTTTTTACTGGTTCGTTTTTTTCATCttttttgtggtt                                                                                             >  1:1027976/1‑141 (MQ=255)                                                      gCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACA‑CCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCttt                                                                                         <  1:452344/151‑1 (MQ=255)                                                       cATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAAT‑AAAAA‑CTCCCCCTATTTTTTCCACTGGGGCGCTTTTTGTATACTTCTAGTGCTtttttt                                                                                        >  1:249340/1‑151 (MQ=255)                                                        aTTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACA‑CCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCttttt                                                                                       <  1:10637/151‑1 (MQ=255)                                                             gATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACA‑CCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑tt                                                                                 <  1:362363/151‑1 (MQ=255)                                                              aTTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTTATTTTTTTCAT‑TAAAA‑CCCCCCCTTTTTTTTTCTATGGGTTGTTTTTTTTTTCTTTTGTGTGCTTCtttttttctt                                                                                  >  1:204576/1‑145 (MQ=255)                                                              aTTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGGGCCACGTTCCTTTTCAATAGTGGCGGTTATATTTTTCAA‑ACAGA‑AATTCCCAATTATTTTCAATGGTTAGCTTTTTTCCTTTTTCGCGTGCTTCttttt                                                                                       >  2:560886/1‑145 (MQ=255)                                                                  cGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTTAT‑GAACA‑CCCCCCCTTTTTTTTTCCGTGGGTTGTTTTTTTCCTCTTTTGTGTGCttttttttttctttcccc                                                                             >  1:1292335/1‑141 (MQ=255)                                                                   gTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTTAAT‑TAAAA‑CCCCCCTTTTTTTTTTCTTTGGGTTGTTTTTTTCATCTTTTTTGTGGTTTTTTTTTTC‑TT‑ccccc                                                                           >  1:111324/1‑150 (MQ=255)                                                                      aGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACA‑CCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TT‑CCACCA‑Taa                                                                      <  1:165527/151‑1 (MQ=255)                                                                                 tGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTACAT‑TAAAA‑CCCCCCTTTTTTTTTTTAGTGGGTTGTTTTTTTCCTCTTTTTAGTGGTTTTTTTTTTCTTT‑CCCCCT‑AAA‑CCGGTTggg                                                            >  2:1074305/1‑150 (MQ=255)                                                                                   aaTATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTTATTTTTTACAT‑GAAAA‑CCCCCCTTATTTTTTTCAGTGGGTTGTTTTTTTCATCTTTTGTGGGGTTCTTTTTTTC‑TT‑TCCCCC‑TAA‑cccggt                                                               >  2:5194/1‑140 (MQ=255)                                                                                     tatCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACAT‑GAAAA‑CCCCCCTTATTTTTTTCAGTGGTTCGTTTATTTCATCTTTCTAGTGCTTCTTTTTTAC‑TT‑CCACCA‑TAA‑CCGGTTCcggggt                                                        >  2:738254/1‑146 (MQ=255)                                                                                      atCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACAT‑GAACACCCCCCTTTTTTTTTTCCAGGGGTTGGTTT‑TTTCATCTTTTTAGTGCTTTTTTTTTTC‑TT‑TCACCC‑TAACCCGGTTggggtt                                                         >  2:6287/1‑146 (MQ=255)                                                                                      atCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGGAATTTTTTACAA‑TAAAA‑ACCCCCTTATTTTTTTCAGGGGGTCGGTTTTTTCCTCTTTTCCGTGGTTTTTTTTTTC‑TTCCCACCC‑AAA‑CCGGTTggggt                                                          >  2:967263/1‑146 (MQ=255)                                                                                          tGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACA‑CCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TT‑CCACCA‑TAA‑CCGGTTCGGGATTCGc                                                     <  1:824384/147‑1 (MQ=255)                                                                                           gTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAAT‑TAAAA‑CCCCCCTTTTTTTTTTCACTGGGTTGGTTTTTTTTTCTTTTTTGGGCTTTTTTTTTTT‑TT‑CCCCCA‑TAA‑CCCGTTTGGGGttt                                                       >  1:791919/1‑143 (MQ=255)                                                                                            tCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTTATTTTTTTAAT‑TAACA‑CCTCCCTTATTATTTTCCGTGGGTCGTTTTTTTCATCTTTTTAGGGCTTTTTTTTTTT‑TT‑TCCCCC‑CAA‑ACCGTTCGGGATTTGCCTcc                                                 >  1:1032818/1‑149 (MQ=255)                                                                                               aGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACA‑CCTCCCTTTTTATTTTCAGTGGGTCGTTTTTTTCCTCTTTTTAGGGCTTCTTTTTTAC‑TT‑CCACCA‑CAA‑CAGGTTTGGGGTTCGCATc                                                  >  2:866518/1‑145 (MQ=255)                                                                                                 ttAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAAT‑TAAAA‑ACCCCCTTTTTTTTTTCTGTGGTTTGTTTTTTTCATTTTTTTTGTGCTTTTTTTTTTT‑TT‑CCACCA‑TAA‑ACGGTTTGGGGTTTGCATTCGtt                                              >  1:792063/1‑147 (MQ=255)                                                                                                 ttAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATAT‑GAAAA‑CCCCCACTTTTTTTTTTTGTGGGTTGTTTTTTTCCTCTTTCGGGTGCTTCTTTTTTTC‑TT‑TCCCCC‑AAA‑CCGGCTCGGGAtttcccacccgt                                              >  2:67777/1‑137 (MQ=255)                                                                                                  tagtagGAGTGCCACCTTCCTTTTCAATAGTGGCGGGTATTTTTTTTATATAAAA‑CCCCCTTTTTTTTTTTC‑CCGGGTCCTTTTTTTTTTCTTTTCTGTGCCTCTTTTTTAC‑TT‑TCACCC‑TAA‑CCCGGTTGGGGTTTGCATCCGTTatac                                          >  1:453225/1‑150 (MQ=255)                                                                                                   agtagGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATAA‑TAACA‑ACCCCCTTTTTTTTTTCAGTGGGTTGTTTTTTTCCTCTTTTTAGTGGTTCTTTTTTTC‑TT‑TCCCCACAAA‑CCGGTTGGGGtttgcctcccgttcac                                           >  2:409703/1‑136 (MQ=255)                                                                                                    gtagGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACA‑CCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TT‑CCACCA‑TAA‑CCGGTTCGGGATTCGCATCCGTTACACTTc                                       <  2:549587/151‑1 (MQ=255)                                                                                                         aGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACA‑CCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TT‑CCACCA‑TAA‑CCGGTTCGGGATTCGCATCCGTTACACTTCGCTCg                                  <  2:1216023/151‑1 (MQ=255)                                                                                                         aGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACA‑CCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TT‑CCACCA‑TAA‑CCGGTTCGGGATTCGCATCCGTTACACTTCGCTCg                                  <  1:744347/151‑1 (MQ=255)                                                                                                           tGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTCAT‑TAAAA‑CCCCCCTTTTTTTTTTCAGTGGTGCGTTTATTTCATTTTTCCTTTGCTTCTTTTTTTC‑TT‑TCACCC‑TAA‑ACCGTTAGGGGTTTGCCTCCCTTCCACTTCTCTCgg                                 >  1:279645/1‑150 (MQ=255)                                                                                                            gCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACA‑CCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TT‑CCACCA‑TAA‑CCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGt                                <  1:1214871/150‑1 (MQ=255)                                                                                                             ccACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACA‑CCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TT‑CCACCA‑TAA‑CCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTaa                              <  1:1052697/151‑1 (MQ=255)                                                                                                              cACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACA‑CCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TT‑CCACCA‑TAA‑CCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAat                             <  1:719509/151‑1 (MQ=255)                                                                                                                    ccttTTCAATAGTGGCGGTAATTTTATACAT‑GAACA‑CCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TT‑CCACCA‑TAA‑CCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTg                       <  1:925239/151‑1 (MQ=255)                                                                                                                    ccttTTCAATAGTGGCGGTAATTTTATACAT‑GAAAA‑ACTCCCTTATTTTTTTCCGGGGTTCGTTTTTTTCCTCTTTTTCGTGGTTTTTTTTTTC‑TT‑TCCCCA‑TAA‑CCCGGTCGGGGTTTGCATCCGTTTCCCTTTTCCCGGTAAAATcc                         >  2:1099544/1‑149 (MQ=255)                                                                                                                      ttttCAATAGTGGCGGTAATTTTTTTAAT‑TAAAA‑ACCCCCTTTTTTTTTTCAGTGGGTTGTTTTTTTCATCTTTTTAGTGCTTTTTTTTTTC‑TT‑TCCCCC‑TAA‑ACCGTTCGGGGTTTGCATCCGTTTCACTTCCCTCCGTTatat                           >  1:66158/1‑145 (MQ=255)                                                                                                                                 ggCGGTAATTTTATACAT‑GAACA‑CCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TT‑CCACCA‑TAA‑CCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCa           <  2:37880/150‑1 (MQ=255)                                                                                                                                   cGGTAATTTTATACAT‑GAACA‑CCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TT‑CCACCA‑TAA‑CCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCg        <  1:459365/151‑1 (MQ=255)                                                                                                                                    ggTAATTTTATACAT‑GAACA‑CCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TT‑CCACCA‑TAA‑CCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGt       >  2:709337/1‑151 (MQ=255)                                                                                                                                       aaTTTTATACAT‑GAACA‑CCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TT‑CCACCA‑TAA‑CCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTAc    >  1:275986/1‑151 (MQ=255)                                                                                                                                          tttATACAT‑GAACA‑CCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TT‑CCACCA‑TAA‑CCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTg  >  1:976653/1‑150 (MQ=255)                                                                                                                                             |                                                                                                                                                         ACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACAT‑GAACA‑CCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TT‑CCACCA‑TAA‑CCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTG  >  NC_000913/1413595‑1413881

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