Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	298,960	0	A	C	40.3%	133.6 / 13.0	68	Y659D (TAT→GAT)	paoC	PaoABC aldehyde oxidoreductase, Moco‑containing subunit
Reads supporting (aligned to +/- strand):  ref base A (35/5);  new base C (0/27);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 4.91e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGAT‑ATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACACAGC  >  NC_000913/298814‑299110                                                                                                                                                     |                                                                                                                                                       ataCACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGAT‑ATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATAcccc                                                                                                                                                    <  1:1285623/151‑1 (MQ=255)  taCACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGAT‑ATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCg                                                                                                                                                   >  1:1366748/1‑151 (MQ=255)              gCGATAGCCGCGCTCACGCCGCACAGCCCCAGCTCACCGACCCCTTTGGCCTTCATCGGGGAGGAT‑ATGGGGTCGGTATCACCCAGAAAAATCCCCTCCTGTTTTGGG‑ATATCCGCAGGACCCGGCACCTCATCCCCCGCCATATCGTGAt                                                                                                                                       <  1:1081811/151‑1 (MQ=255)               ggATGGCCGCGCTCACGCCGCACAGGCCCACCTCCCCGACCCTTTTGGCCTTCATGGGGGAGGAT‑ATGGGGTCGTTATCACCCAGGAAAACCCCCCCCTTTTTGGGA‑ATATCCGCATGACCCGGCACCTCATACCCCGCCATATCGTGAtt                                                                                                                                      <  1:1223913/150‑1 (MQ=255)               cGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCCCCGACACCTTTGGCCTTCATCGGGGAGGAT‑ATGGGGTCGGTATCATCCAGGAAAATCCCCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATCCCCCGCCATATCGTGAtt                                                                                                                                      <  2:1135129/151‑1 (MQ=255)                   ggCCCCGCTCACGCCCCACAGCCCCAGCTCCCCGACCCCTTTGGCCTCCATCGGGGAGGAT‑ATGGGGTGGGTATCATCCAGAAAAATCCCCTCCTTTTTTGGG‑ATATCCGCATGAACCGGCCCCTCATCCCCCGCCATATCGTGATTAACg                                                                                                                                  <  1:1125972/150‑1 (MQ=255)                   cgCCGCGCTCCCGCCCCACAGCCCCGGCTCACCGCCTCCTTTGGCCTTCATCGGGGGGGAT‑AGGGGGTCGGTATCATCCAGGAAAATCCCCTCCTGTTTTGGA‑ATACCCGCATGAGCCGGCACCTCATACCCCGCCATATCGTGATTAACg                                                                                                                                  <  2:242739/150‑1 (MQ=255)                    gCCGCGCTCACGCCCCACAGCCCCAGCTCCCCGACCCCTTTGGCCTTCATCGGGGGGGAT‑ATGGGGTCGGTACCACCCAGAAAAACCCCCTCCTGTTTGGGG‑ATACCCGCATAACCCGGCACCTCATCCCCCGCCATATCGTGATTAACGa                                                                                                                                 <  1:146087/151‑1 (MQ=255)                       gcgcTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGAT‑ATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGt                                                                                                                              >  2:806157/1‑151 (MQ=255)                         gcTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGAT‑ATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGt                                                                                                                              >  1:903436/1‑149 (MQ=255)                            cACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGAT‑ATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGccc                                                                                                                         >  2:1329960/1‑151 (MQ=255)                              cgcggcACAGGCCCAGCTCACCGACACCTTTGGCCTTCATGGGGGAGGAT‑AGGGGGTCGGTATCATCCAGGAAAATCCCCTCCGGTTTTGGG‑ATATCCGAATGAACCGGCACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCaa                                                                                                                       <  2:1265538/151‑1 (MQ=255)                               gccgcACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGAT‑ATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCaa                                                                                                                       >  2:97715/1‑150 (MQ=255)                                 cgcACGGGCCCAGCTCCCCGACCCCTTTGGCCTTAATCGGGGAGGAT‑ATGGGGTCGGTATCATCCAGAAAAACCCCCTCCTTTTTTGGA‑ATACCCGAAGGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACg                                                                                                                    <  1:950191/151‑1 (MQ=255)                                 cgcACAGGCCCAGCTCCCCGACACTTTGGGCTTTCATGGGGGAGGAT‑ATGGGGTCGTTATCACCCAGAAAAACCACCTCTTTTTTTGGA‑ATACCCGCATGAACCGCCCCCTTACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACg                                                                                                                    <  1:760635/151‑1 (MQ=255)                                           cAGCTCACCGACCCCTTTGGCCTTCATCGGGGAGGAT‑ATGGGGTCGGTATCATCCAGGAAAATCCCCTCCTGTTTTGGG‑ATATCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCAcc                                                                                                          <  1:677426/151‑1 (MQ=255)                                                cACCGACACCTTTGGCCTTCATCGGGGAGGAT‑ATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGcc                                                                                                       >  2:1129852/1‑149 (MQ=255)                                                cACCGACACCTTTGGCCTTCATCGGGGAGGAT‑ATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGcc                                                                                                       >  2:1139452/1‑149 (MQ=255)                                                cACCGACACCTTTGGCCTTCATCGGGGAGGAT‑ATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGcc                                                                                                       >  2:716564/1‑149 (MQ=255)                                                 aCCGACACCTTTGGCCTTCATCGGGGAGGAT‑ATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGc                                                                                                    >  1:206920/1‑151 (MQ=255)                                                    aacccttttGCCCTTCATGGGGGAGGAT‑AGGGGGTGGTTATC‑CCCAAGAAAACCCCCCCCCTTTTTGGGAATACCCGAATGACCCGGCCCCCCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctcc                                                                                                 <  1:798208/145‑1 (MQ=255)                                                      cccttttGCCCTTCATCGGGGAGGAT‑ATGGGGCCGGTATCACCCAGAAAAACCCCCTCCTGTTTGGGG‑ATATCCGCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctcctc                                                                                               <  1:546537/147‑1 (MQ=255)                                                           ttGGCCTTCATCGGGGAGGAT‑AGGGGGTCGGTATCATCCAGAAAAACCCCCTCCTTTTTTGGG‑ATACCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATc                                                                                           <  2:764407/150‑1 (MQ=255)                                                           ttGCCCTTTATCGGGG‑GGATAATGGGGCGGTTATCATCCAGGAAAATCCCCTCTTGTTTGGGG‑ATATCCGCGTAACCCGGCACCCCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATc                                                                                           <  2:434996/150‑1 (MQ=255)                                                                cTTCATCGGGGAGGAT‑ATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGcc                                                                                     >  1:437647/1‑151 (MQ=255)                                                                  tCATCGGGGAGGAT‑ATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCcgc                                                                                   >  2:1190730/1‑151 (MQ=255)                                                                  tCATCGGGGAGGAT‑ATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCcgc                                                                                   >  2:298069/1‑151 (MQ=255)                                                                  tCATCGGGGAGGAT‑ATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCcgc                                                                                   >  2:951223/1‑151 (MQ=255)                                                                    aTCGGGGAGGAT‑ATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCgcgc                                                                                 >  2:469392/1‑151 (MQ=255)                                                                      cGGGGAGGAT‑ATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgccc                                                                               >  1:1304283/1‑151 (MQ=255)                                                                       ggggAGGAT‑ATGGGGTCGTTACCACCCAGAAAAACCCCCCCCTTTTTGGGG‑ATATCCCCATGAACCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgccca                                                                              <  2:1452504/151‑1 (MQ=255)                                                                           aggaT‑ATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgcc                                                                          >  1:326925/1‑151 (MQ=255)                                                                           aggaT‑AGGGGGTCGTTATCATCCGGGAAAACCCCCTCCTGTTTTGGG‑ATCCCCGCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgcc                                                                          <  1:730287/151‑1 (MQ=255)                                                                            ggaT‑ATGGGGTCGTTATCACCCAGAAAAACCCCCTCTTTTTTGGGG‑ATACCCACATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgccc                                                                         <  1:885501/151‑1 (MQ=255)                                                                             ggttagggggTCGGTATCACCCAGAAAAACCCCCTCCTGTTTGGGG‑ATACCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgccc                                                                         <  1:1015014/145‑1 (MQ=255)                                                                             gaT‑ATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgccc                                                                         >  2:777179/1‑150 (MQ=255)                                                                             gaT‑ATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgccc                                                                         >  2:599047/1‑150 (MQ=255)                                                                              at‑atGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCAGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgcccat                                                                       >  2:275682/1‑151 (MQ=255)                                                                                 atGGGGTCGGTATCATCCAGGAAAATCCCCTCCTGTTTTGGG‑ATTTCCGCATAACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAg                                                                     <  2:430488/151‑1 (MQ=255)                                                                                 atGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgcccat                                                                       >  1:935727/1‑149 (MQ=255)                                                                                   ggggTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCACATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTc                                                                   >  2:1240604/1‑151 (MQ=255)                                                                                   ggggTCGGTATCATCCAGAAAAATCCCCTCCTGTTTTGGG‑ATACCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTc                                                                   <  1:341878/151‑1 (MQ=255)                                                                                       tCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTg                                                               >  2:1267030/1‑151 (MQ=255)                                                                                        cGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTgc                                                              >  2:1057845/1‑151 (MQ=255)                                                                                              catcCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAAt                                                        >  1:134674/1‑151 (MQ=255)                                                                                              catcCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAAt                                                        >  2:180152/1‑151 (MQ=255)                                                                                                cccAGAAAAACCCCCTCCTGTTTTGGG‑ATACCCGCATAACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGa                                                      <  1:112661/150‑1 (MQ=255)                                                                                                  cAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGAcc                                                    >  2:238340/1‑151 (MQ=255)                                                                                                     aaaaccccccccctTTTTGGGA‑ATCTCCGCATAACCCGCCACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTgg                                                 <  1:468630/140‑1 (MQ=255)                                                                                                           cACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcgc                                           >  2:662755/1‑151 (MQ=255)                                                                                                                   gTTTGGGA‑ATACCCGAATGAACCGGCCCCCCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGtt                                     <  1:17151/149‑1 (MQ=255)                                                                                                                    tttgGGA‑ATACCCGCATGACCCGGCACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGtttt                                   <  1:985927/150‑1 (MQ=255)                                                                                                                     ttgggg‑ATACCCGCAGGAACCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTc                                  <  1:401740/147‑1 (MQ=255)                                                                                                                       tGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCgg                                >  2:261292/1‑150 (MQ=255)                                                                                                                       tGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCgg                                >  2:1172938/1‑150 (MQ=255)                                                                                                                        ggg‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGAt                              >  1:429516/1‑151 (MQ=255)                                                                                                                                ccGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGAt                       >  2:1079374/1‑151 (MQ=255)                                                                                                                                  gCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGc                     >  2:1096153/1‑151 (MQ=255)                                                                                                                                  gCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGc                     >  1:1128049/1‑151 (MQ=255)                                                                                                                                  gCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGc                     >  1:682820/1‑151 (MQ=255)                                                                                                                                    aTGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGt                   >  1:1198053/1‑151 (MQ=255)                                                                                                                                        cccGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTg               <  2:552542/150‑1 (MQ=255)                                                                                                                                         cggcccccccATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCACCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTGCGGATTCAGGATGCGTACTgc              <  1:182181/142‑1 (MQ=255)                                                                                                                                         ccGGCACCTAATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTgc              <  2:903439/151‑1 (MQ=255)                                                                                                                                            cccccTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTgcagc           <  1:455547/148‑1 (MQ=255)                                                                                                                                               ccccATGCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGAAGCGTCCTGCAGcaca        <  1:707073/148‑1 (MQ=255)                                                                                                                                                  cATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcacacac     <  2:1366744/151‑1 (MQ=255)                                                                                                                                                  cATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcacacac     <  2:538140/151‑1 (MQ=255)                                                                                                                                                     aCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACACAGc  >  2:812619/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                       ATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGAT‑ATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACACAGC  >  NC_000913/298814‑299110

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