Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,712,518	(C)6→7	intergenic (+360/‑251)	nth → / → dtpA	DNA glycosylase and apyrimidinic (AP) lyase (endonuclease III)/dipeptide and tripeptide permease A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,712,512	1	.	C	97.6%	282.9 / ‑0.6	82	intergenic (+354/‑257)	nth/dtpA	DNA glycosylase and apyrimidinic (AP) lyase (endonuclease III)/dipeptide and tripeptide permease A
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (42/38);  total (43/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.83e-01

GAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGTCTGA  >  NC_000913/1712374‑1712645                                                                                                                                            |                                                                                                                                      gAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCa                                                                                                                                        <  1:591383/139‑1 (MQ=255) gAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCa                                                                                                                                        <  2:875675/139‑1 (MQ=255)                 aTTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTcagtgc                                                                                                                        <  1:857260/139‑1 (MQ=255)                  tttAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTcagtgca                                                                                                                       >  2:1098986/1‑139 (MQ=255)                  tttAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑ccccccctct                                                                                                                             <  1:1164046/132‑5 (MQ=255)                  tttAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑ccccccctct                                                                                                                             >  2:1164046/1‑128 (MQ=255)                                             gACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGtttttt                                                                                            >  2:382399/1‑139 (MQ=255)                                             gACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGtttttt                                                                                            >  2:650370/1‑139 (MQ=255)                                             gACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGtttttt                                                                                            >  1:30087/1‑139 (MQ=255)                                             gACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑ccccc                                                                                                                                  <  1:83955/100‑1 (MQ=255)                                             gACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑ccccc                                                                                                                                  >  2:83955/1‑100 (MQ=255)                                                cTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGaa                                                                                         <  2:30087/139‑1 (MQ=255)                                                cTACAAAACATTACACTGGCTATTTTTCAGAAACTAGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGaa                                                                                         >  2:468593/1‑139 (MQ=255)                                                             cacTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTg                                                                            >  2:396484/1‑139 (MQ=255)                                                             cacTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑ccccc                                                                                                                                  >  2:723997/1‑84 (MQ=255)                                                             cacTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATAACCCCCCCTCTCAGTGCAGTGaaa                                                                                                                 >  1:881984/1‑102 (MQ=255)                                                             cacTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAACAGCAGAACAGATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGaaa                                                                                                                 <  2:881984/102‑1 (MQ=255)                                                             cacTGGCTATTGTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑ccccc                                                                                                                                  <  1:723997/84‑1 (MQ=255)                                                                tGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGAt                                                                         >  1:1228582/1‑139 (MQ=255)                                                                tGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGAt                                                                         >  1:221813/1‑139 (MQ=255)                                                                  gCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGTAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATaa                                                                       >  1:1333909/1‑139 (MQ=255)                                                                  gCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACg                                                                               >  2:746933/1‑131 (MQ=255)                                                                  gCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACg                                                                               <  1:746933/131‑1 (MQ=255)                                                                  gCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACg                                                                               <  1:1237948/131‑1 (MQ=255)                                                                  gCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACg                                                                               >  2:1237948/1‑131 (MQ=255)                                                                  gCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATaa                                                                       >  1:379246/1‑139 (MQ=255)                                                                     aTTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑ccccc                                                                                                                                  >  1:880095/1‑76 (MQ=255)                                                                     aTTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑ccccc                                                                                                                                  <  2:880095/76‑1 (MQ=255)                                                                           cAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCg                                                              <  1:382399/139‑1 (MQ=255)                                                                            aGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGaa                                                                                         >  2:1392893/1‑111 (MQ=255)                                                                            aGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGaa                                                                                         <  1:1392893/111‑1 (MQ=255)                                                                            aGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAAc                                                                                >  2:1419286/1‑120 (MQ=255)                                                                            aGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAAc                                                                                <  1:1419286/120‑1 (MQ=255)                                                                              aaaCTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGaa                                                           <  2:1045079/139‑1 (MQ=255)                                                                                aCTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATg                                                         >  1:410278/1‑139 (MQ=255)                                                                                aCTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATg                                                         >  1:1373960/1‑139 (MQ=255)                                                                                         tctcACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATg                                                         >  2:913210/1‑130 (MQ=255)                                                                                         tctcACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATg                                                         <  1:913210/130‑1 (MQ=255)                                                                                           tcACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTccc                                              <  1:468593/139‑1 (MQ=255)                                                                                           tcACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTccc                                              >  1:667559/1‑139 (MQ=255)                                                                                            cACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑ccccc                                                                                                                                  <  2:846132/53‑1 (MQ=255)                                                                                            cACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑ccccc                                                                                                                                  >  1:846132/1‑53 (MQ=255)                                                                                              ctcgTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTAAGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGc                                           <  2:946980/136‑1 (MQ=255)                                                                                              cTAGTAAAAACGCTAAATAGTAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTGGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGc                                           <  2:1333909/139‑1 (MQ=255)                                                                                              cTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTa                                                                                  >  1:28069/1‑100 (MQ=255)                                                                                              cTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTa                                                                                  <  2:28069/100‑1 (MQ=255)                                                                                              cTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGc                                           <  1:650370/139‑1 (MQ=255)                                                                                              cTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGc                                           >  1:435683/1‑139 (MQ=255)                                                                                              cTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGc                                           >  2:1266535/1‑139 (MQ=255)                                                                                              cTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑ccccc                                                                                                                                  <  2:736260/51‑1 (MQ=255)                                                                                              cTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑ccccc                                                                                                                                  >  1:736260/1‑51 (MQ=255)                                                                                               tAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGcc                                          <  1:396484/139‑1 (MQ=255)                                                                                                aGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCt                                         <  1:1098986/139‑1 (MQ=255)                                                                                                aGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCt                                         <  2:221813/139‑1 (MQ=255)                                                                                                aGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCt                                         >  2:165033/1‑139 (MQ=255)                                                                                                aGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCt                                         <  2:1130040/139‑1 (MQ=255)                                                                                                   aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTg                                                     <  1:155152/124‑1 (MQ=255)                                                                                                   aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTg                                                     >  2:155152/1‑124 (MQ=255)                                                                                                   aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGcc                                          <  1:574872/135‑1 (MQ=255)                                                                                                   aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGcc                                          >  2:574872/1‑135 (MQ=255)                                                                                                   aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGc                                           <  1:540029/134‑1 (MQ=255)                                                                                                   aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGc                                           >  2:540029/1‑134 (MQ=255)                                                                                                   aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAg                                      >  2:1174734/1‑139 (MQ=255)                                                                                                   aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAg                                      >  2:1098305/1‑139 (MQ=255)                                                                                                    aaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCt                                         <  2:383731/135‑1 (MQ=255)                                                                                                    aaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCt                                         >  1:383731/1‑135 (MQ=255)                                                                                                     aaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGtt                                    <  1:1174734/139‑1 (MQ=255)                                                                                                       aCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTtat                                  <  2:379246/139‑1 (MQ=255)                                                                                                       aCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTtat                                  >  1:8015/1‑139 (MQ=255)                                                                                                        cGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTtata                                 <  2:281433/139‑1 (MQ=255)                                                                                                        cGCTAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTtata                                 <  1:1266535/139‑1 (MQ=255)                                                                                                           tAAATAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGt                              >  1:1355247/1‑139 (MQ=255)                                                                                                            aaaTAGCAGAACATATCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTa                             >  1:41607/1‑139 (MQ=255)                                                                                                            aaaTAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCCTCTCAGTGCAGTGAAAAAATc                                                                                                            >  2:234788/1‑59 (MQ=255)                                                                                                            aaaTAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCCTCTCAGTGCAGTGAAAAAATc                                                                                                            <  1:234788/59‑1 (MQ=255)                                                                                                                          atCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAGCGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTaa                            >  2:919036/1‑126 (MQ=255)                                                                                                                          atCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAGCGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTaa                            <  1:919036/126‑1 (MQ=255)                                                                                                                          atCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACg                                                                               <  2:380024/75‑1 (MQ=255)                                                                                                                          atCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACg                                                                               >  1:380024/1‑75 (MQ=255)                                                                                                                          atCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAg                                      >  2:614396/1‑116 (MQ=255)                                                                                                                          atCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAg                                      <  1:614396/116‑1 (MQ=255)                                                                                                                          atCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTaa                            <  2:1377617/126‑1 (MQ=255)                                                                                                                          atCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTaa                            >  1:1377617/1‑126 (MQ=255)                                                                                                                          atCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGattat                    <  2:870792/134‑1 (MQ=255)                                                                                                                          atCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGattat                    >  1:870792/1‑134 (MQ=255)                                                                                                                          atCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACaa               >  2:1345871/1‑139 (MQ=255)                                                                                                                           tCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTccc                                              <  1:27733/107‑1 (MQ=255)                                                                                                                           tCGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTccc                                              >  2:27733/1‑107 (MQ=255)                                                                                                                            cGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTccc                                              <  2:693859/106‑1 (MQ=255)                                                                                                                            cGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTccc                                              >  1:693859/1‑106 (MQ=255)                                                                                                                            cGCCGTTAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAg             >  2:1033286/1‑139 (MQ=255)                                                                                                                                 ttAAGCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACtt        >  2:187362/1‑139 (MQ=255)                                                                                                                                     gCCATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGTCt    <  1:1345871/139‑1 (MQ=255)                                                                                                                                       cATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGTCTGa  <  2:41607/139‑1 (MQ=255)                                                                                                                                       cATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGTCTGa  <  1:187362/139‑1 (MQ=255)                                                                                                                                       cATCACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGTCTGa  <  2:1228582/139‑1 (MQ=255)                                                                                                                                            |                                                                                                                                      GAATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCA‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGTCTGA  >  NC_000913/1712374‑1712645

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