Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,687,212	G→C	9.1%	P66R (CCG→CGG)	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,212	0	G	C	9.1%	210.6 / 10.1	66	P66R (CCG→CGG)	glnB	nitrogen regulatory protein PII‑1
Reads supporting (aligned to +/- strand):  ref base G (27/33);  new base C (4/2);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 4.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

CGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGTCATGCCGGTAATACCGACTTCGGCCAGTGCTTCGCG  >  NC_000913/2687064‑2687360                                                                                                                                                     |                                                                                                                                                     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                                                                                                                                        >  2:1267888/1‑128 (MQ=255)                              cACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTACCGGACGACattgtcg                                                                                                                                        <  1:1267888/134‑7 (MQ=255)                                       ggATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATcca                                                                                                               <  1:388950/150‑1 (MQ=255)                                        gATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATccac                                                                                                              <  1:1377895/150‑1 (MQ=255)                                         aTGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATccacc                                                                                                             <  1:1408790/150‑1 (MQ=255)                                         aTGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATccacc                                                                                                             <  1:845038/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:17524/150‑1 (MQ=255)                                                                                     aCCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCg                                                                 <  2:22925/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                                      >  2:30855/1‑150 (MQ=255)                                                                                                                aataatGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGt                                      >  1:1133375/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)                                                                                                                                                    cGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGTCATGCCGGTAATACCGACTTCGGCCAGTGCTTcgcg  <  1:742111/150‑1 (MQ=255)                                                                                                                                                    cGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGTCATGCCGGTAATACCGACTTCGGCCAGTGCTTcgcg  <  1:1334342/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     CGCGGTTTAAATTGCCGCGTCGTCCTCTTCACCGGTACGGATGCGAATGACCCGTGCCACGTCAAAGACGAAGATTTTACCGTCACCGATTTTGCCGGTTTGCGCCGTGCGAATAATGGTATCGACACAGGTATCGACAATGTCGTCCGGTACGACAATCTCAATTTTCACTTTCGGCAGAAAATCCACCATATACTCCGCGCCGCGGTACAGCTCGGTATGGCCTTTCTGGCGACCAAAGCCTTTCACTTCGGTCACCGTCATGCCGGTAATACCGACTTCGGCCAGTGCTTCGCG  >  NC_000913/2687064‑2687360

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