Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	897,968:1	+A	13.4%	intergenic (+39/‑21)	potI → / → ybjO	putrescine ABC transporter membrane subunit PotI/putative inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	897,968	1	.	A	13.4%	207.2 / 35.9	67	intergenic (+39/‑21)	potI/ybjO	putrescine ABC transporter membrane subunit PotI/putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base . (32/26);  new base A (5/4);  total (37/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

CATTGATCCTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCTGATG  >  NC_000913/897823‑898114                                                                                                                                                   |                                                                                                                                                   cATTGATCCTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATg                                                                                                                                                >  1:96044/1‑150 (MQ=255)       tCCTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCgg                                                                                                                                          <  1:256991/150‑1 (MQ=255)        ccTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGc                                                                                                                                         >  2:880509/1‑150 (MQ=255)           cGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGACATG‑CGCGGCAtt                                                                                                                                      <  2:644793/150‑1 (MQ=255)           cGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGACATG‑CGCGGCAtt                                                                                                                                      <  1:1278629/150‑1 (MQ=255)               gcgGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAATACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACTTAATGCCGCGACATG‑CGCGGCAttatgt                                                                                                                                  <  2:374334/150‑5 (MQ=255)                 ggTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGACATG‑CGCGGCAttatgtag                                                                                                                                <  2:1278175/150‑7 (MQ=255)                 ggTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTc                                                                                                                                <  2:254623/150‑1 (MQ=255)                 ggTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTc                                                                                                                                >  1:216263/1‑150 (MQ=255)                   tCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCAt                                                                                                                              <  1:629463/150‑1 (MQ=255)                   tCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCAt                                                                                                                              <  1:457588/150‑1 (MQ=255)                             cGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGaa                                                                                                                    >  2:646704/1‑150 (MQ=255)                                   tATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTg                                                                                                              >  2:1144934/1‑150 (MQ=255)                                    aTCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTgg                                                                                                             >  1:325624/1‑150 (MQ=255)                                     tCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTggg                                                                                                            >  1:681417/1‑150 (MQ=255)                                       gCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGAt                                                                                                          >  2:520353/1‑150 (MQ=255)                                        ccTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGAtt                                                                                                         >  1:251057/1‑150 (MQ=255)                                           tgTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGAttttt                                                                                                      <  2:1042824/149‑1 (MQ=255)                                            gTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGAtttttt                                                                                                     >  2:1084738/1‑150 (MQ=255)                                             tATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTa                                                                                                    >  1:1366619/1‑150 (MQ=255)                                               tCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAg                                                                                                  <  2:96044/150‑1 (MQ=255)                                                       cTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATc                                                                                          >  1:702011/1‑150 (MQ=255)                                                          gcgcAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTc                                                                                       <  1:665696/150‑1 (MQ=255)                                                             cAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCAtc                                                                                    >  1:811144/1‑150 (MQ=255)                                                             cAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGCCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCAtc                                                                                    >  2:113634/1‑150 (MQ=255)                                                                aaaaaCAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCa                                                                                 <  1:1021595/150‑1 (MQ=255)                                                                 aaaaCAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCAt                                                                                >  2:752030/1‑150 (MQ=255)                                                                 aaaaCAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCAt                                                                                >  2:468769/1‑150 (MQ=255)                                                                               gcgACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAAtgt                                                                  >  2:1233193/1‑150 (MQ=255)                                                                               gcgACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAAtgt                                                                  >  1:545086/1‑150 (MQ=255)                                                                               gcgACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAAtgt                                                                  <  1:994504/150‑1 (MQ=255)                                                                                  aCATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGcc                                                               <  2:1115949/150‑1 (MQ=255)                                                                                  aCATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTATTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGcc                                                               <  1:633648/150‑1 (MQ=255)                                                                                    aTCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTg                                                             <  1:520353/150‑1 (MQ=255)                                                                                      ccAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTgcg                                                           <  2:406076/150‑1 (MQ=255)                                                                                      ccAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTgcg                                                           >  2:1383644/1‑150 (MQ=255)                                                                                              cAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAATCGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCa                                                   <  1:853116/150‑1 (MQ=255)                                                                                              cAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCa                                                   <  2:325624/150‑1 (MQ=255)                                                                                                aGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAgg                                                 >  2:958201/1‑150 (MQ=255)                                                                                                aGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAgg                                                 <  1:638969/150‑1 (MQ=255)                                                                                                   cGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTgg                                              >  2:380356/1‑150 (MQ=255)                                                                                                     tGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCg                                                                           >  2:814209/1‑119 (MQ=255)                                                                                                     tGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCg                                                                           >  2:809138/1‑119 (MQ=255)                                                                                                     tGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCg                                                                           <  1:814209/119‑1 (MQ=255)                                                                                                     tGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCg                                                                           <  1:809138/119‑1 (MQ=255)                                                                                                      ggCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTggcgg                                           >  2:1217030/1‑150 (MQ=255)                                                                                                      ggCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTggcgg                                           <  1:468769/150‑1 (MQ=255)                                                                                                        cTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGgcggcg                                         >  1:1213462/1‑150 (MQ=255)                                                                                                           aaGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTc                                      <  1:1144934/150‑1 (MQ=255)                                                                                                                acTAAAATTTGCCAACCTGGCTACATAATGCCGCA‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCat                                 <  1:257376/150‑1 (MQ=255)                                                                                                                 cTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAgg                                                 <  1:562479/133‑1 (MQ=255)                                                                                                                 cTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAgg                                                 >  2:562479/1‑133 (MQ=255)                                                                                                                 cTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCatt                                >  1:213816/1‑150 (MQ=255)                                                                                                                   aaaaTTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCattat                              >  1:63259/1‑150 (MQ=255)                                                                                                                          gCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATc                       >  1:1279106/1‑150 (MQ=255)                                                                                                                           ccAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATcc                      >  1:1407287/1‑150 (MQ=255)                                                                                                                              aCCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTg                   <  2:749045/150‑1 (MQ=255)                                                                                                                                     tACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGa            >  2:1405568/1‑150 (MQ=255)                                                                                                                                       cATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACg          <  1:1084738/150‑1 (MQ=255)                                                                                                                                         caaTGCCGCGACATG‑CGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTg        >  2:1080835/2‑150 (MQ=255)                                                                                                                                         caaTGCCGCGACATG‑CGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTg        >  1:1150962/2‑150 (MQ=255)                                                                                                                                          aaTGCCGCGACATG‑CGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGc       >  2:719569/1‑150 (MQ=255)                                                                                                                                          aaTGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGc       <  2:681417/150‑1 (MQ=255)                                                                                                                                          aaTGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGc       <  1:1405568/150‑1 (MQ=255)                                                                                                                                           aTGCCGCGACATG‑CGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCt      >  2:840763/1‑150 (MQ=255)                                                                                                                                            tGCCGCGACATG‑CGCGGCATTGTTTTCATGGAAGACGAAACGTAGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCtg     >  2:270736/1‑150 (MQ=255)                                                                                                                                            tGCCGAG‑CAAGACGCGGAATTGTTTTCATGGAAGACCAAACGTTGGCATTTTTTAAGAAAACATCTTCATATCATGCTCCCCTGAATGTGCAAGAGATGGTGCAGGTCGCGGCGCTCGACATTATTAAGAAACTTGGCCTCGACGTGcag     >  2:1072588/1‑148 (MQ=255)                                                                                                                                               cgcg‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCtgatg  >  2:983981/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                                   CATTGATCCTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCTGATG  >  NC_000913/897823‑898114

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