Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,030,732	A→G	5.7%	intergenic (‑64/+27)	tusE ← / ← yccA	sulfur transfer protein TusE/modulator of FtsH protease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,030,732	0	A	G	5.7%	215.7 / 2.8	70	intergenic (‑64/+27)	tusE/yccA	sulfur transfer protein TusE/modulator of FtsH protease
Reads supporting (aligned to +/- strand):  ref base A (24/42);  new base G (2/2);  total (26/44)
Fisher's exact test for biased strand distribution p-value = 6.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

CAGCCAGCGGCTCACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGTCTCACCGCCATGAATGATGTTGCTGG  >  NC_000913/1030584‑1030880                                                                                                                                                     |                                                                                                                                                     cagccagcGGCTCACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAg                                                                                                                                                     <  2:829847/150‑1 (MQ=255)     cagcGGCTCACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTa                                                                                                                                                 <  2:451018/150‑1 (MQ=255)     cagcGGCTCACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTa                                                                                                                                                 <  1:245912/150‑1 (MQ=255)      agcGGCTCACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTAt                                                                                                                                                <  1:1127973/150‑1 (MQ=255)       gcGGCTCACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGGGGCTATg                                                                                                                                               <  2:173552/150‑1 (MQ=255)       gcGGCTCACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATg                                                                                                                                               <  2:1061024/150‑1 (MQ=255)       gcGGCTCACTCCACTGGCTACTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATg                                                                                                                                               <  1:755036/150‑1 (MQ=255)          gCTCACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGaaa                                                                                                                                            <  1:1348388/150‑1 (MQ=255)           cTCACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCCGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAg                                                                                                                                           <  1:213966/150‑1 (MQ=255)               cTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGtt                                                                                                                                       >  2:75548/1‑150 (MQ=255)                    cTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCA‑TGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCggg                                                                                                                                 <  1:1147142/150‑1 (MQ=255)                     tGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCggg                                                                                                                                 <  2:45405/150‑1 (MQ=255)                     tGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCggg                                                                                                                                 <  1:24969/150‑1 (MQ=255)                       gctgctTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTg                                                                                                                               <  2:981621/150‑1 (MQ=255)                       gctgctTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTg                                                                                                                               <  2:765671/150‑1 (MQ=255)                        ctgctTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTga                                                                                                                              <  2:262986/150‑1 (MQ=255)                           ctttctttGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCTCGCTTTATGCGGGGCTTTgtttttt                                                                                                                                         <  1:56446/136‑7 (MQ=255)                           ctttctttGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCTCGCTTTATGCGGGGCTTTgtttttt                                                                                                                                         >  2:56446/1‑130 (MQ=255)                           ctttctttGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGAt                                                                                                                           >  2:770605/1‑150 (MQ=255)                               ctttGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAAt                                                                                                                       <  1:639151/150‑1 (MQ=255)                                  tGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATcgc                                                                                                                    <  2:759807/150‑1 (MQ=255)                                     gaTAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGc                                                                                                                 <  2:1089789/150‑1 (MQ=255)                                                gTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGccc                                                                                                      >  1:879308/1‑150 (MQ=255)                                                     cGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATc                                                                                                                      <  1:1364757/129‑1 (MQ=255)                                                     cGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATc                                                                                                                      >  2:1364757/1‑129 (MQ=255)                                                       tttCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGc                                                                                               <  2:318389/150‑1 (MQ=255)                                                         tCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTc                                                                                             >  2:1193538/1‑150 (MQ=255)                                                         tCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTc                                                                                             >  2:1088540/1‑150 (MQ=255)                                                           tATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTcag                                                                                           <  2:795714/150‑1 (MQ=255)                                                                ctTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCGGAATGCTCAGCAGGc                                                                                      <  2:672697/150‑1 (MQ=255)                                                                  ttACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTg                                                                                    <  2:1092822/150‑1 (MQ=255)                                                                    aCCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGAc                                                                                  <  1:828353/150‑1 (MQ=255)                                                                     ccTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACg                                                                                 >  2:811571/1‑150 (MQ=255)                                                                       ttCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGaa                                                                               >  1:827816/1‑150 (MQ=255)                                                                          gAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGAt                                                                            <  2:358198/150‑1 (MQ=255)                                                                             gATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGAtgtt                                                                         >  2:644607/1‑150 (MQ=255)                                                                              aTCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGAtgttg                                                                        <  1:931184/150‑1 (MQ=255)                                                                                  gCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACa                                                                    >  2:1047691/1‑150 (MQ=255)                                                                                    aTAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGc                                                                  <  2:103955/150‑1 (MQ=255)                                                                                     tAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCg                                                                 <  2:65744/150‑1 (MQ=255)                                                                                      aaTTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGa                                                                >  2:1263893/1‑150 (MQ=255)                                                                                         tattaCTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAAc                                                             <  1:1193538/150‑1 (MQ=255)                                                                                          attaCTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACa                                                            >  1:724373/1‑150 (MQ=255)                                                                                           ttaCTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACAt                                                           <  1:770605/150‑1 (MQ=255)                                                                                             actacCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATAc                                                         >  1:1390657/1‑150 (MQ=255)                                                                                              ctacCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACa                                                        <  1:1088540/150‑1 (MQ=255)                                                                                                acCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAgg                                                      >  2:164269/1‑150 (MQ=255)                                                                                                 ccTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGc                                                     <  1:610742/150‑1 (MQ=255)                                                                                                      aTCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGAGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAAcc                                                >  1:885523/1‑150 (MQ=255)                                                                                                                cggcAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGaa                                      >  2:166557/1‑150 (MQ=255)                                                                                                                cggcAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGaa                                      <  2:879308/150‑1 (MQ=255)                                                                                                                        tAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTc                              <  2:803804/150‑1 (MQ=255)                                                                                                                         aaCAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCg                             <  2:827816/150‑1 (MQ=255)                                                                                                                         aaCAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCg                             <  2:724373/150‑1 (MQ=255)                                                                                                                           cAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGtc                           >  2:646335/1‑150 (MQ=255)                                                                                                                             aaaaaCAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGtctc                         >  1:471786/1‑150 (MQ=255)                                                                                                                                  cataGCCTCGCTTTATGCGGGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGTCTCACCGc                    >  1:1281802/4‑150 (MQ=255)                                                                                                                                  cAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGTCTCACCGc                    >  1:849697/1‑150 (MQ=255)                                                                                                                                        cccGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGTCTCACCGCCatgaa              >  1:102246/1‑150 (MQ=255)                                                                                                                                        cccGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGTCTCACCGCCatgaa              >  1:1419213/1‑150 (MQ=255)                                                                                                                                           gCATAAATCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTATCGAAGCCCATAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGTCACTAATATAGTTCGTCTCACCGCCatgaatga           <  2:650003/150‑1 (MQ=255)                                                                                                                                           gCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGTCTCACCGCCatgaatga           >  1:1116943/1‑150 (MQ=255)                                                                                                                                           gCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGTCTCACCGCCatgaatga           <  1:811571/150‑1 (MQ=255)                                                                                                                                            cATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGTCTCACCGCCATGAatgat          <  2:1116943/150‑1 (MQ=255)                                                                                                                                            cATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGTCTCACCGCCATGAatgat          <  1:1263893/150‑1 (MQ=255)                                                                                                                                              tAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGTCTCACCGCCATGAATGATGt        <  1:555512/150‑1 (MQ=255)                                                                                                                                              tAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGTCTCACCGCCATGAATGATGt        >  1:37138/1‑150 (MQ=255)                                                                                                                                                aaGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAACCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGc                                                                  <  1:1142003/90‑1 (MQ=255)                                                                                                                                                aaGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAACCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGc                                                                  >  2:1142003/1‑90 (MQ=255)                                                                                                                                                    gAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGTCTCACCGCCATGAATGATGTTGCTgg  <  1:1047691/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     CAGCCAGCGGCTCACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTATCTCTTTACCTTCGAAGATCAGCATAATTATTACTACCTTAATCAGACTGGCGGCAGTGTAACAAAAAACAAAGCCCCGCATAAAGCGAGGCTATGAAAGTGTTAGCGGGTGAGATTAATCGCGGCTAGCGAAGCCCAGAATGCTCAGCAGGCTGACGAAGATGTTGTACAGCGAAACATACAGGCTAACCGTGGCACGAATATAGTTCGTCTCACCGCCATGAATGATGTTGCTGG  >  NC_000913/1030584‑1030880

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