Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,687,210	C→A	9.8%	D67Y (GAC→TAC)	glnB ←	nitrogen regulatory protein PII‑1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,687,210	0	C	A	9.8%	197.6 / 8.5	61	D67Y (GAC→TAC)	glnB	nitrogen regulatory protein PII‑1
Reads supporting (aligned to +/- strand):  ref base C (25/30);  new base A (4/2);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 4.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

TGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGTCATGCCGGTAATACCGACTTCGGCCAGTGCTTCGC  >  NC_000913/2687061‑2687359                                                                                                                                                      |                                                                                                                                                      tGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATgtcgtc                                                                                                                                                       >  1:1345017/1‑150 (MQ=255)  gACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTcc                                                                                                                                                      <  1:378586/150‑1 (MQ=255)  gACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTcc                                                                                                                                                      >  2:378586/1‑150 (MQ=255)    cgcgGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCgg                                                                                                                                                    >  1:1394046/1‑150 (MQ=255)     gcgGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGt                                                                                                                                                   >  2:1262965/1‑150 (MQ=255)        gTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCTACAATGTCGTCCGGTACg                                                                                                                                                <  2:951934/150‑1 (MQ=255)           tAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACa                                                                                                                                             >  2:1408790/1‑150 (MQ=255)            aaaTTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACaa                                                                                                                                            >  1:22925/1‑150 (MQ=255)            aaaTTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACaa                                                                                                                                            >  1:189790/1‑150 (MQ=255)            aaaTTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACaa                                                                                                                                            >  1:17524/1‑150 (MQ=255)              aTTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAAtc                                                                                                                                          <  1:1383215/150‑1 (MQ=255)               ttGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTACCGGACGACattgt                                                                                                                                         <  2:913814/150‑5 (MQ=255)                 gCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCa                                                                                                                                       >  2:723098/1‑150 (MQ=255)                       tcgtcCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTc                                                                                                                                 <  2:1394046/150‑1 (MQ=255)                          tcCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACt                                                                                                                              <  2:504794/150‑1 (MQ=255)                          tcCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACt                                                                                                                              <  1:784439/150‑1 (MQ=255)                             tctTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTc                                                                                                                           <  1:315432/150‑1 (MQ=255)                               ttCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCgg                                                                                                                         >  1:351643/1‑150 (MQ=255)                                tCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGc                                                                                                                        <  2:351643/150‑1 (MQ=255)                                 cACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTACCGGACGACattgtcg                                                                                                                                       <  1:1267888/134‑7 (MQ=255)                                 cACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTACCGGACGACattgtcg                                                                                                                                       >  2:1267888/1‑128 (MQ=255)                                          ggATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATcca                                                                                                              <  1:388950/150‑1 (MQ=255)                                           gATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATccac                                                                                                             <  1:1377895/150‑1 (MQ=255)                                            aTGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATccacc                                                                                                            <  1:845038/150‑1 (MQ=255)                                            aTGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATccacc                                                                                                            <  1:1408790/150‑1 (MQ=255)                                                  aTGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATAc                                                                                                      <  1:936013/150‑1 (MQ=255)                                                    gACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTc                                                                                                    >  2:129249/1‑150 (MQ=255)                                                     aCCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTcc                                                                                                   >  2:1172856/1‑150 (MQ=255)                                                               cGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGt                                                                                         <  1:723098/150‑1 (MQ=255)                                                                gTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTa                                                                                        >  2:1189187/1‑150 (MQ=255)                                                                  cAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACa                                                                                      <  2:97850/150‑1 (MQ=255)                                                                                ttACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCt                                                                        <  1:802060/150‑1 (MQ=255)                                                                                    cGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCt                                                                    <  1:775357/150‑1 (MQ=255)                                                                                       cACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGc                                                                 >  2:1363543/1‑150 (MQ=255)                                                                                        aCCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCg                                                                <  2:22925/150‑1 (MQ=255)                                                                                        aCCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCg                                                                <  2:17524/150‑1 (MQ=255)                                                                                           gATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGAcc                                                             >  2:172984/1‑150 (MQ=255)                                                                                              tttGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCaaa                                                          <  2:335126/150‑1 (MQ=255)                                                                                               ttcccGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAg                                                         >  2:810358/4‑150 (MQ=255)                                                                                                tGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGc                                                        >  1:648886/1‑150 (MQ=255)                                                                                                     gTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTc                                                   <  2:148889/150‑1 (MQ=255)                                                                                                      tttGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCa                                                  <  2:1345017/150‑1 (MQ=255)                                                                                                      tttGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCa                                                  <  1:1262965/150‑1 (MQ=255)                                                                                                          cgcCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTc                                              <  1:1326106/150‑1 (MQ=255)                                                                                                             cGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGt                                           <  1:129249/150‑1 (MQ=255)                                                                                                              gTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTc                                          >  1:382377/1‑150 (MQ=255)                                                                                                                 cGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCAcc                                       >  1:234112/1‑150 (MQ=255)                                                                                                                 cGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCAcc                                       >  1:1245890/1‑150 (MQ=255)                                                                                                                   aataatGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGt                                     >  1:1133375/1‑150 (MQ=255)                                                                                                                   aataatGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGt                                     >  2:30855/1‑150 (MQ=255)                                                                                                                     taatGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGTCa                                   <  2:648886/150‑1 (MQ=255)                                                                                                                         ggTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGTCATGcc                               >  2:272223/1‑150 (MQ=255)                                                                                                                           tATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGTCATGCCgg                             <  2:594411/150‑1 (MQ=255)                                                                                                                            aTCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGTCATGCCGGt                            >  1:804441/1‑150 (MQ=255)                                                                                                                            aTCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGTCATGCCGGt                            <  2:382377/150‑1 (MQ=255)                                                                                                                                acacAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGTCATGCCGGTAATa                        >  1:1219185/1‑150 (MQ=255)                                                                                                                                         ttgagattGTCGTACCGGACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGTCATGCCGGTAATACCGACTTCg               >  1:802345/8‑150 (MQ=255)                                                                                                                                          tgagattGTCGTACCGGACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGTCATGCCGGTAATACCGACTTCgg              >  2:1218421/7‑150 (MQ=255)                                                                                                                                          cGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGTCATGCCGGTAATACCGACTTCgg              >  2:934345/1‑150 (MQ=255)                                                                                                                                           gagattGTCGTACCGGACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGTCATGCCGGTAATACCGACTTCGGc             >  2:224257/6‑150 (MQ=255)                                                                                                                                           gACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGTCATGCCGGTAATACCGACTTCGGc             >  2:983064/1‑150 (MQ=255)                                                                                                                                            aCAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCTGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGTCATGCCGGTAATACCGACTTCGGcc            <  1:526936/150‑1 (MQ=255)                                                                                                                                            aCAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGTCATGCCGGTAATACCGACTTCGGcc            <  1:696336/150‑1 (MQ=255)                                                                                                                                               aTGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGTCATGCCGGTAATACCGACTTCGGCCAGt         <  1:1172856/150‑1 (MQ=255)                                                                                                                                                 gtcgtcCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGTCATGCCGGTAATACCGACTTCGGCCAGTGc       >  2:742111/1‑150 (MQ=255)                                                                                                                                                   cgtcCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCa                                                  <  2:331962/105‑1 (MQ=255)                                                                                                                                                   cgtcCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCa                                                  >  1:331962/1‑105 (MQ=255)                                                                                                                                                      ccGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGTCATGCCGGTAATACCGACTTCGGCCAGTGCTTcgc  >  2:1125432/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                      TGACGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGTCATGCCGGTAATACCGACTTCGGCCAGTGCTTCGC  >  NC_000913/2687061‑2687359

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: ‑