Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,276,716	G→T	G416C (GGT→TGT)	pdeC →	c‑di‑GMP phosphodiesterase PdeC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,276,716	0	G	T	98.0%	176.1 / ‑3.0	51	G416C (GGT→TGT)	pdeC	c‑di‑GMP phosphodiesterase PdeC
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (30/20);  total (30/21)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.71e-01

TAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCT  >  NC_000913/4276471‑4276959                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    tAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAACTTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTAc                                                                                                                                                                                                                                                >  2:406736/1‑251 (MQ=255)      gATTATGTTAGTGATAAGGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGaa                                                                                                                                                                                                                                             <  1:228377/249‑1 (MQ=255)             ttATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAAtt                                                                                                                                                                                                                                    >  1:4200/1‑251 (MQ=255)                 tGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATtgatg                                                                                                                                                                                                                                >  2:37079/1‑251 (MQ=255)                   aTAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTgatgat                                                                                                                                                                                                                              <  1:37833/251‑1 (MQ=255)                     aaTGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGAttt                                                                                                                                                                                                                            >  1:501033/1‑251 (MQ=255)                        gTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTg                                                                                                                                                                                                                          <  2:124841/250‑1 (MQ=255)                           ttCCGCGATCTGGGCGCTTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGGCACGGTTGATAGCGCGAATCAATCAGAAAACCGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGTAGTGGCAATTGATGATTTTGGTAt                                                                                                                                                                                                                      >  1:461093/1‑251 (MQ=255)                              cgcgATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTgg                                                                                                                                                                                                                   >  2:553450/1‑251 (MQ=255)                                 gATCTGGGCGATTAACTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCACTTGATGATTTTGGTATTGGCt                                                                                                                                                                                                                 >  1:427548/1‑250 (MQ=255)                                               ccTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATa                                                                                                                                                                                                  <  1:116802/251‑1 (MQ=255)                                                cTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATa                                                                                                                                                                                                  >  1:111370/1‑250 (MQ=255)                                                   gCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCt                                                                                                                                                                                              <  1:546217/251‑1 (MQ=255)                                                     aaCACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCtt                                                                                                                                                                                             >  2:527382/1‑250 (MQ=255)                                                              aGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCt                                                                                                                                                                                                                 <  2:538121/221‑1 (MQ=255)                                                              aGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCt                                                                                                                                                                                                                 >  1:538121/1‑221 (MQ=255)                                                               gATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGa                                                                                                                                                                                   >  1:477230/1‑250 (MQ=255)                                                                 tCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATg                                                                                                                                                                                >  1:505543/1‑251 (MQ=255)                                                                  cGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGCTGATTTTGGTATTGGCTACTCTAACTTGCATAACATTAAATCATTGAATGt                                                                                                                                                                               >  2:128369/1‑251 (MQ=255)                                                                            cTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTg                                                                                                                                                                     >  1:88087/1‑251 (MQ=255)                                                                             gaTTAACCTGTCGGCCGCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATGAAATTTGAAGTGACTGAGCATGCATTTCGTGATGTTTACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTGAAATCATTGAATGGCTATATTTTg                                                                                                                                                                     <  2:550593/249‑1 (MQ=255)                                                                              aTTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATAtttt                                                                                                                                                                      <  2:501033/248‑1 (MQ=255)                                                                                                gATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTttgt                                                                                                                                                 <  2:535979/251‑1 (MQ=255)                                                                                                             cACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTGTATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGacc                                                                                                                                    <  2:418444/251‑1 (MQ=255)                                                                                                                            gAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGACGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGt                                                                                                                     <  1:295386/251‑1 (MQ=255)                                                                                                                                             gagCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTTACAAAATGACTCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTTAAACGCTGACCACCCACAAAACCAGTCAGTTGATTGCGGaata                                                                                                    <  2:303076/251‑3 (MQ=255)                                                                                                                                             gagCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAAGCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAaca                                                                                                    <  2:331314/251‑1 (MQ=255)                                                                                                                                              agCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAacac                                                                                                   >  1:137523/1‑251 (MQ=255)                                                                                                                                                   tACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACAcatcat                                                                                              >  2:220154/1‑251 (MQ=255)                                                                                                                                                       cgGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCg                                                                                            >  2:446482/1‑249 (MQ=255)                                                                                                                                                                  cagcaAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGGCGCACAg                                                                               <  2:88087/251‑1 (MQ=255)                                                                                                                                                                    gcaAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGGCGCACAGcc                                                                             >  1:443299/1‑251 (MQ=255)                                                                                                                                                                            aattTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGGCGCACAGCCTGGGGTTa                                                                     <  2:485810/251‑1 (MQ=255)                                                                                                                                                                                tGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGGCGCACAGCCTGGGGTTaaaa                                                                  <  1:189726/250‑1 (MQ=255)                                                                                                                                                                                             catgcatTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAgg                                                    >  2:288114/1‑251 (MQ=255)                                                                                                                                                                                             catgcatTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAgg                                                    >  1:481018/1‑251 (MQ=255)                                                                                                                                                                                                  atTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCg                                               <  1:406736/251‑1 (MQ=255)                                                                                                                                                                                                  atTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCg                                               >  1:454462/1‑251 (MQ=255)                                                                                                                                                                                                   tttCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCg                                               >  1:95178/1‑250 (MQ=255)                                                                                                                                                                                                          aTGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACTgag                                       <  2:137523/251‑1 (MQ=255)                                                                                                                                                                                                            gTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACTgagg                                      >  1:48149/1‑250 (MQ=255)                                                                                                                                                                                                            gTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGAGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAACACAGTCATTTGATTGCGGAAAACATCATCGAGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACTgaggg                                     >  2:520100/1‑250 (MQ=255)                                                                                                                                                                                                                aCAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACTGAGGAGc                                   >  1:389624/1‑249 (MQ=255)                                                                                                                                                                                                                  aaaaGGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCAGTTGATTGCGGAACACATCATCGAGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAgtt                               <  2:340356/251‑3 (MQ=255)                                                                                                                                                                                                                   aaaTGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGtt                              >  2:408075/1‑251 (MQ=255)                                                                                                                                                                                                                     aTGACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGTTaa                            <  1:553450/251‑1 (MQ=255)                                                                                                                                                                                                                       gACGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGTTAACt                          >  1:508937/1‑251 (MQ=255)                                                                                                                                                                                                                         cGCCGATTATTCTGGCTTTCCGCCAGGCATGTTACGATGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGGAGCACAGCCTGGGGTTAAAAACGATCGCTGACGGCGTCGAAACTGAGGAGAAGGATAACTgg                        >  1:281706/1‑251 (MQ=255)                                                                                                                                                                                                                               ttattCTGGCTTTCAGCACGGCATGTTTCGAAGTGGCAATTGATGCTTTTGGTCTTGGCTACTCTAACTTGCATAACCTTAAATCATTGGATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGCAACGCTGACCACCCACCAACCCAGTCATTTGATTGCGGAACACATCATCGAGCTGACGCACAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCGAACCTGAGGAGCAGGTTAACAGGCTccgc                  >  2:533390/1‑251 (MQ=255)                                                                                                                                                                                                                                    cTGGCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAAc              >  2:409576/1‑250 (MQ=255)                                                                                                                                                                                                                                      ggCTTTCCGCCAGGCATGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAGTCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAAc              <  2:369246/248‑1 (MQ=255)                                                                                                                                                                                                                                                caggcaggTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCt  >  2:132921/1‑250 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    TAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTCGATATTTTGAAAATCGACAAATCGTTTGTTGAAACGCTGACCACCCACAAAACCAGTCATTTGATTGCGGAACACATCATCGAGCTGGCGCACAGCCTGGGGTTAAAAACGATCGCTGAAGGCGTCGAAACTGAGGAGCAGGTTAACTGGCTGCGCAAACGCGGCGTGCGCT  >  NC_000913/4276471‑4276959

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