Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,560,362	G→A	intergenic (+496/+260)	rtcR → / ← glpG	DNA‑binding transcriptional activator RtcR/rhomboid protease GlpG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,560,362	0	G	A	94.1%	182.7 / 1.8	70	intergenic (+496/+260)	rtcR/glpG	DNA‑binding transcriptional activator RtcR/rhomboid protease GlpG
Reads supporting (aligned to +/- strand):  ref base G (4/0);  new base A (35/29);  total (41/29)
Fisher's exact test for biased strand distribution p-value = 1.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

GAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTT‑TACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAGCCAT‑TGGGGATTTGCTCCGCCACTTTGCGGGCGATGCG‑CTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTT  >  NC_000913/3560118‑3560607                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                        gAACTGGGAGAAAAAACCGAGCGTCGCTGCGCCAATGATCCCGCCATCGCGGCGGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTT‑TACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAGCCAT‑TGGGGATTTGCTCCGCCACTTTGCGGACGATgc                                                                                                                                                                                                                                                   <  1:327850/251‑1 (MQ=255)      gggAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTT‑TACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAGCCAT‑TGGGGATTTGCTCCGCCACTTTGCGGCAGATGCG‑ctct                                                                                                                                                                                                                                             >  2:661250/1‑251 (MQ=255)        gagaTAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTT‑TACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAGCCAT‑TGGGGATTTGCTCCGCCACTTTGCGGACGATGCG‑CTCttt                                                                                                                                                                                                                                           >  1:361144/1‑251 (MQ=255)        gagaTAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTT‑TACCATCAACGTATTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGGGGTGCCGCTATC‑GATAAACAGCGTCGTGACAT‑TGGGGATTTGCTCCGCCAATTTGCTGACGATACC‑CTCttt                                                                                                                                                                                                                                           >  2:524232/1‑251 (MQ=255)          gaTAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTT‑TACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAGCCAT‑TGGGGATTTGCTCCGCCACTTTGCGGACGATGCG‑CTCttttt                                                                                                                                                                                                                                         >  1:161879/1‑251 (MQ=255)            aaaaagtatgcggcGCTTCGCCAATGATCCCGCCATCGCGGCTGCTTAATTCGCCCCCGGCGAAAATGATGCGAAAATCTTCTTT‑TACCATCAACGTGTTAGCAACGTTGAGATTGTGGGTGACAATGCGCAAATTG‑CTGGGATTGAGCAGTGCGGGCGCTCCCGTTTCCGGCGTGGTGCCGATAGC‑GATAAACAGCTTCGAGCCAT‑TGGGGATTTGCTCCGCCACTTTGCGGACGATGCG‑CTCTTTttct                                                                                                                                                                                                                                       <  2:642539/238‑1 (MQ=255)              aaaTCGAGTGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTT‑TACCATCAACGTGTTAGCAACGTTGAGATTTTTGGTGACAATGCGCAAATTG‑CGGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGTGTGGGGCCGATATC‑GATAAACAGCGTCGAGCCAT‑TGGGGATTTGCTCCGCCACTTTGCGGACGATGCG‑CTCTTTttcttc                                                                                                                                                                                                                                     <  1:45309/251‑1 (MQ=255)                aTCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCTGCTGCGTAATTCGCCACCGTCGAGAATGATGCGAAAATCTTCTTT‑TACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAGCCAT‑TGGGGATTTGCTCCGCCACTTTGCGGACGATGCG‑CTCTTTTTCTTCgg                                                                                                                                                                                                                                   <  1:432963/251‑1 (MQ=255)                   gAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTT‑TACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACGGCGTCGAGCCAT‑TGGGGATTTGCTCCTCCACTTTGCGGACGATGCG‑CTCTTTTTCTTCGGTCt                                                                                                                                                                                                                                >  1:290642/1‑251 (MQ=255)                   gAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTT‑TACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAGCCAT‑TGGGGATTTGCTCCGCCACTTTGCGGACGATGCG‑CTCTTTTTCTTCGGTCt                                                                                                                                                                                                                                <  1:288166/251‑1 (MQ=255)                           cTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTT‑TACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAGCCAT‑TGGGGATTTGCTCCGCCACTTTGCGGACGATGCG‑CTCTTTTTCTTCGGTCTGGGTGGcc                                                                                                                                                                                                                        >  1:646193/1‑251 (MQ=255)                                        cccGCCATAGCGGCGGCGTAAATCGCCACCGGCGAGAATGATGCGAAAATCTGCTTG‑TACCAGCAACGTGTTAGCAACGGTGAGATTGTGGGGGACAATGCGCAAATTG‑CTGTGATTGAACAGTGCGTGCGCTCCCGCTTCGGGCGTGGTGCCGATATG‑GATAAACAGAGTCGAGC‑ATGTGGGGATTTGCTCCGCCACGTTGCGGACGATGCG‑CTCTTTTTCTGCGGTCTGGGTGGCCTTGCGATCGTGcc                                                                                                                                                                                                           <  2:280577/251‑1 (MQ=255)                                         ccGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTT‑TACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAGCCAT‑TGGGGATTTGCTCCGCCACTTTGCGGACGATGCG‑CTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCa                                                                                                                                                                                                          <  1:196238/251‑1 (MQ=255)                                                    gCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTT‑TACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAGCCAT‑TGGGGATTTGCTCCGCCACTTTGCGGACGATGCG‑CTCTTTTTCTTCGGTCTGGGTGGCTTTGCGATCGTGCCACGGCGTGTTa                                                                                                                                                                                                >  2:28515/1‑250 (MQ=255)                                                     cTGCGTCATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTT‑TACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAGCCAT‑TGGGGCTTTGCTCCGCCACTTTGCGGCCGATGCG‑CTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAc                                                                                                                                                                                              >  1:456548/1‑251 (MQ=255)                                                       tcttAATTCGCCCCCGGCGAGAATGATGCGAAAATCTACCTT‑TCCCATCAACGTGTTATCAACGTTGATAGTGTTGGTGACAATGCGCAAATTG‑TTGTTATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAGCCAT‑TGGGGATTGGCTCCGCCACTTTGCGGACGATGCG‑CTCTTTTTCTGCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAAcc                                                                                                                                                                                             <  1:29864/247‑1 (MQ=255)                                                        cGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTT‑TACCAGCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAGCCAT‑TGGGGATTTGCTCCGCCACTTTGCGGACGATGCG‑CTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGa                                                                                                                                                                                           <  1:253193/251‑1 (MQ=255)                                                           aaTTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTT‑TACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAGACAT‑TGGGGATTTGCTCCGCCACTTTGCGGGCGATGCG‑CTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACt                                                                                                                                                                                        >  2:320346/1‑251 (MQ=255)                                                               cGCCACCGGCGAGAATGATGCGAAAATCTTCTTT‑TACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAGCCAT‑TGGGGATTTGCTCCGCCACTTTGCGGACGATGCG‑CTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGa                                                                                                                                                                                     >  1:425431/1‑250 (MQ=255)                                                                 cctccGGCGAGAATGATGCGAAAATCTTATTT‑TACCAACAACGGGTGAGAAACGTTTAGATGGTGGGTGAAAAGGCGCAAAGTG‑CGGGGAGTGAGCAGTGGGTGCGCTACCGCTGCCGGCGTGGTGCCGATATC‑GAGAAACAGCGGCGAGCCAT‑GGGGGATTTGCTCCGCCACTTTGCGGACGATGCG‑CTCTTTTTCTTCGGTCGGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAgg                                                                                                                                                                                  <  2:425431/248‑1 (MQ=255)                                                                  cACCGGCGAGAATGATGCGAAAATCTTCTTT‑TACCATCCACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATGGAGCAGTGCGTGCGCGACCGCTTCCGGCGGGGTGCCGATATC‑GATAAACAGCGTCGCGCCAT‑TGGGGATTTGCTCCGCCACTTTGCGGACGATGCG‑CTCTTTTTCTTCGGTCTGGGTGGCCTTGAGATCGTGCCACGGCGTGTTTACCGAACTGGaaagc                                                                                                                                                                                 >  2:286926/1‑248 (MQ=255)                                                                    ccgtggagaaggtggccaaaagtttcttT‑TCCCATAAACTTGTTAGATACGTTGAGATTTGTGGGGAGAA‑GCCGAAATTGGCTGTGATGGATCAGTGCGGGCGTTTCCCCGTCCGGCGTGGTCCCGATATC‑GATAAACAGCGGCGAGCCAT‑GGGGGATTTGCTCCGCCTCTTTGCGGACGATGCG‑CTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACGGGAAGGCAg                                                                                                                                                                               <  2:352276/228‑1 (MQ=255)                                                                     cGGCGAGAATGATGCGAAAATCTTCTTT‑TACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAGCCAT‑TGGGGATTTGCTCCGCCACTTTGCGGACGATGCG‑CTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAgc                                                                                                                                                                              <  1:126949/251‑1 (MQ=255)                                                                      ggCGAGAATGATGCGAAAATCTTCTTT‑TACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAGCCAT‑TGGGGATTTGATCCGCCACTTTGCGGACGATGCG‑CTCTTTTTCTTCGGTCTGGGTGGCCTTGCGCTCGTGCCACGGCGTGTTAACCGAACTGGAGGGCcgcg                                                                                                                                                                             >  2:612547/1‑251 (MQ=255)                                                                           gaATGATGCGCAAATCTTCTTT‑TACCATCAACGTGATAGCAACGATGAGATAGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGAGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAGCCAG‑TGGGGATTTGCTCCGCCACTTTGCGGACGATGCG‑CTCTTTTTCTTCGGTATGTGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAATGCAGcggcgcc                                                                                                                                                                        >  2:247186/1‑250 (MQ=255)                                                                           aaaTGATGCGAAAATCT‑CTTTGTACCACCAACGTGTTAGCAGCGTTGAGATTGTTGGTGACATTGAGTAAATTG‑CGGGAATGGAGCAGTGCGGGCGTGACCGTTTCCGGCGTGGTGCCGAAATC‑AAAAAACAGCGTCGAGCCAT‑GGGGGATTTGCTCCGCCACTTTGCGGACGATGCG‑CTCTTTTTCTTCGGTCTGGGTGGCCGTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCa                                                                                                                                                                        <  2:125676/250‑1 (MQ=255)                                                                            aaTGATGCGAAAATCTTCTTT‑TACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAGCCAT‑TGGGGATTTGCTCCGCCACTTTGCGGACGATGCG‑CTCttt                                                                                                                                                                                                                                           >  2:575447/1‑183 (MQ=255)                                                                            aaTGATGCGAAAATCTTCTTT‑TACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAGCCAT‑TGGGGATTTGCTCCGCCACTTTGCGGACGATGCG‑CTCttt                                                                                                                                                                                                                                           <  1:575447/183‑1 (MQ=255)                                                                             atgatgAGAACATCTTCTTT‑TACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAGCCAT‑TGGGGATTTGCTCCGCCACTTTGCGGACGATGCG‑CTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGACCTGGAAGGCGGCGCCGCAcc                                                                                                                                                                      >  2:628107/1‑251 (MQ=255)                                                                              tgatgCGAAAATCTTCTTT‑TACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAGCCAT‑TGGGGATTTGCTCCGCCACTTTGCGGACGATGCG‑CTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCg                                                                                                                                                                     >  1:146888/1‑251 (MQ=255)                                                                               gatgtgtaatacttcttT‑TACAAGCAACTTTTTTGCATCGTTGAGATTGTGGGGAACAATGCTAAAATTG‑CTGTGAGTGAGCAGTGCGTGCGCTGCCGCTTGTGGCTTGGTGCCGATCGT‑GATAAACACGGTTGAGCCAT‑TGGGGTGTGGCTCCGCCACTTTGCGGACGATGCG‑CTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTTTTAACCGAACTGGAAGGCATCGCCGCActgc                                                                                                                                                                    <  2:511316/240‑4 (MQ=255)                                                                                 tgCGAAAATCTTCTTT‑TACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAGCCAT‑TGGGGATTTGCTCCGCCACTTTGCGGACGATGCG‑CTCTTTTTCTTCGGTCTGGTTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGCAGGCAGCGCCGCACCGccc                                                                                                                                                                  >  2:626734/1‑250 (MQ=255)                                                                                   cGAAAATCTTCTTT‑TACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGAAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGCTATC‑GATAAACAGCGTCGAGCCAT‑TGGGGATTTGCTCCGCCACTTTGCGGACGATGCG‑CTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGACGGCCGCGCCGCACCGCCatg                                                                                                                                                                >  1:280050/1‑251 (MQ=255)                                                                                    gAAAATCTTCTTT‑TACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAGCCAT‑TGGGGATTTGCTCCGCCACTTTGCGGACGATGCG‑CTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCatg                                                                                                                                                                >  1:616010/1‑250 (MQ=255)                                                                                        aTCTTCTTT‑TACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAGCCAT‑TGGGGATTTGCTCCGCCACTTTGCGGACGATGCG‑CTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCACCATGATGGc                                                                                                                                                           >  1:146671/1‑251 (MQ=255)                                                                                            tcttT‑TACCATCAACGTGTTAGCACCGTTGAGATTGTTGGTGACAAGGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACCGCTCCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAGCCAT‑GGGGGTTTTGCTCCGCCACTTTGCGGACGATGCG‑CTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGc                                                                                                                                                           <  1:163558/247‑1 (MQ=255)                                                                                                t‑tACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCATTGCGTGCTCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAGCCAT‑TGTGGATTTGCTCCGCCACTTTGCGGACGATGCG‑CTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGcgcc                                                                                                                                                                           <  1:633410/227‑1 (MQ=255)                                                                                                t‑tACCATCAACGAGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGAAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAGCCAT‑TGGGGATTTGATCCGCCACTTTGCGGACGATGCG‑CTCTTTTTCATCGGTCTGGGTGGCCTTGCGATAGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGcgcccgatcggaagagcggcgggtaggg                                                                                                                                                   >  2:633410/1‑227 (MQ=255)                                                                                                             tgtTAGCAACGTTGAGATTGTTGGTGCCAATGCGCCAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAGCCAT‑TGGGGATTTGCTCCGCCACTTTGCGGACGATGCG‑CTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTTCTcc                                                                                                                                       >  2:91801/1‑251 (MQ=255)                                                                                                                    aaCGTTGCGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAGCCAT‑TGGGGATTTGCTCCGCCACTTTGCGGACGATGCG‑CTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGTAGGCGGCGCCGCACCGCCCTGATGGCGCAGGATCCGGTTTTGCACCGccggcc                                                                                                                                >  2:272664/1‑247 (MQ=255)                                                                                                                    aaCGTTGAGATTTTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAGCCAT‑TGGGGATTTGCTCCGCCACTTTGCGGACGATGCG‑CTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGACCTGGAAGGCAACGCCCCACCGCCATGGTGGCGCAGGATAAGGTTTTGCTCCGCCAGCt                                                                                                                                >  1:158893/1‑251 (MQ=255)                                                                                                                          gagaTTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGCGCCAT‑TGGGGATTTGCTCCGCCACTTTGCGGACGATGCG‑CTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCCCCAGCTCATTg                                                                                                                           >  2:559403/1‑250 (MQ=255)                                                                                                                                       aCAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAGCCAT‑TGGGGATTTGCTCCGCCACTTTGCGGACGATGCG‑CTCTTTTTCATCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAAt                                                                                                             >  2:149887/1‑251 (MQ=255)                                                                                                                                        cAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAGCCAT‑TGGGGATTTGCTCCGCCACTTTGCGGACGATGCG‑CTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAAt                                                                                                             >  1:475267/1‑250 (MQ=255)                                                                                                                                              gcAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAGCCAT‑TGGGGATTTGCTCCGCCACTTTGCGGACGATGCG‑CTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGc                                                                                                      >  1:550661/1‑251 (MQ=255)                                                                                                                                                aaaTTG‑CTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAGCCAT‑TGGGGATTTGCTCCGCCACTTTGCGGACGATGCG‑CTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCgg                                                                                                    >  1:545581/1‑251 (MQ=255)                                                                                                                                                  atgtgcagttaTTGAGGATTGCGTGCGTTCCGGCTTCGGACGTGTTGCCGATTCC‑GATAAACGTCGTCGAGCCTG‑TGGGAATTTGCTCGGCCACTTTGGGGACTATGCG‑CTCTTTTTCGTCGGTCTGGGTGGCCTGGCGATCGTGCCACGGCGTGTTAACCGAACTTGAAGGCGGGGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCggg                                                                                                   <  1:58485/241‑1 (MQ=255)                                                                                                                                                              tGAGCAGTGCGTGCGCTACCGATTCCGGCGTGGTGCCGCTATC‑GATAAACAGCGTCGAGCCAT‑TGGGGATTTGCTCCGCCACTTTGCGGGCGATGCG‑CTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGACCTGGACGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGaa                                                                                        >  2:93648/1‑250 (MQ=255)                                                                                                                                                                aGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAGCCAT‑TGGGGATTTGCTCCGCCACTTTGCGGACGATGCG‑CTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGCAGGCCGCGCCTCACCGCCATGCTGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGCGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATg                                                                                     >  2:456398/1‑251 (MQ=255)                                                                                                                                                                        gtgcgCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAGCCAT‑TGGGGATTTGCTCCGCCACTTTGCGGGCGATGCG‑CTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCGGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAcc                                                                              >  1:55481/1‑250 (MQ=255)                                                                                                                                                                        gtgcgCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAGCCAT‑TGGGGATTTGCTCCGCCACTTTGCGGACGATGCG‑CTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCatgat                                                                                                                                                              >  1:407381/1‑170 (MQ=255)                                                                                                                                                                                  gctcccGGTGTGGTGCCGATTCC‑GAAAACCGGCGTCGGGCTGT‑TGGGTATTTGCTCCGCCACTTTGCGGACGATGCG‑CTCTTTTTCTTCGTTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGt                                                                   <  1:334442/247‑1 (MQ=255)                                                                                                                                                                                   cTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAGCCAT‑TGGGGATTTGCTCCGCCACTTTGCGGGCGATGCG‑CTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGCAATGCTCTACCAGCTCTTCGGt                                                                   >  2:332452/1‑250 (MQ=255)                                                                                                                                                                                     tCCGGCGTGGGGCC‑ATATCCGAGAAACGGCGGCGAGCCAT‑GGGGGATTTGTTCCGCCACTTTGCGGACGATGCG‑CTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTACCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTa                                                                  <  1:91801/249‑1 (MQ=255)                                                                                                                                                                                     tCCGGCGGGGTGCCGATATC‑GAAAAACAGGGTCGACCCTT‑TGGGGTTTTGCTCCGCCACTTTGCGGACGATGCG‑CTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTa                                                                  <  1:661250/249‑1 (MQ=255)                                                                                                                                                                                          cGTGGTGCCGATATC‑GATAAACAGCGTCGAGCCAT‑TGGGGATTTGCTCCGCCACTTTGCGGACGATGCG‑CTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACAt                                                           >  1:594904/1‑251 (MQ=255)                                                                                                                                                                                           gtggtgCCGATATC‑GATAAACAGCGTCGAGCCAT‑TGGGGATTTGCTCCGCCACTTTGCGGACGATGCG‑CTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACAt                                                           >  1:412722/1‑250 (MQ=255)                                                                                                                                                                                              gtgCCGATATC‑GATAAACAGCGTCGAGCCAT‑GGGGGATTTGCCCCGCCCCTTGGCGGACGATGCG‑CTCTTTTTCTTCGGTCTGGGGGGCCTTGCGATCGTGCCCCGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAcc                                                       <  1:534823/251‑1 (MQ=255)                                                                                                                                                                                               tgCCGATATC‑GATAAACAGCGTCGAGCCAT‑TGGGGATTTGCTCCGCCACTTTGCGGACGATGCG‑CTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTAc                                                                 >  1:215420/1‑240 (MQ=255)                                                                                                                                                                                                            aaacctggctcgctcctt‑tGGGAATTGGCTCCGCCACTTGGTGGACGA‑GCGGCTTTTTTTTCTCGGTTTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTACCCGACCTGGAAGGCAGCCCCCAACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAcc                                          <  1:247186/234‑1 (MQ=255)                                                                                                                                                                                                                 aGCGTCGAGCCAT‑TGGGGATTTGCTCCGCCACTTTGCGGACGATGCG‑CTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGGCGGAGAAATGCTCTACCAGCTCTTCGGTACTGTCATAACCCTGCTGTTTAACCAGTTc                                     >  2:464120/1‑251 (MQ=255)                                                                                                                                                                                                                        cgCCAT‑TGGGGTTTTGCTCCGCCACTTTGCGGACGATGCG‑CGCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGataata                              <  1:28515/250‑1 (MQ=255)                                                                                                                                                                                                                         gCCAT‑TGGGTGTTTGCCCCGCCCCTTTGTGGACGATGCG‑CTCTTTTTCTTCGGTCTGGGTGGTTTGGCGATCCTGCCAGGGCGTGTTACCCGAACGGGAAGGCAGCGCTGCACCGCCATGATGGCGCAGGATCAGGTTTTGCGCCGCCAGTTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATAc                             <  1:589359/251‑1 (MQ=255)                                                                                                                                                                                                                         gCCAT‑TGGGGATTTGCTCCGCCACTTTGCGGACGATGCG‑CTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCCGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATAc                             >  1:32561/1‑251 (MQ=255)                                                                                                                                                                                                                                ggggATTTGCTCCGCCACTTTGCGGACGATGCG‑CTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGACCTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGGCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTCACAGCTCTTCGGTACTGACTTAACCCTGCAGTTTAAACAGTTCGCTAATACCGAtgt                       >  2:467512/1‑251 (MQ=255)                                                                                                                                                                                                                                  ggttttcccccGCCACTTTGCGGACGATGCG‑CTCTTTTTCTTCGTGCTGGGTGCCCTTGCGACCGTGCCGCGGCGTTTTACCCTAACTGGAAGGCAGCGCCGCCCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGa                     <  2:594904/242‑1 (MQ=255)                                                                                                                                                                                                                                     tttccTCCGCCACTTTGCGGACGTGGCG‑CTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCTTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACg                   <  1:628107/246‑1 (MQ=255)                                                                                                                                                                                                                                        gCTCCGCCACTTTGCGGACGATGCG‑CTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAACctgct                                                  >  1:91883/1‑216 (MQ=255)                                                                                                                                                                                                                                            cGCCACTTTGCGGACGATGCG‑CTCTTTTTCTTCGGTCGGGTTGGCCTTGCGATCGGGCCCCGGCGTGTTACCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgttt           <  1:286926/251‑1 (MQ=255)                                                                                                                                                                                                                                                     gCGGACGATGCG‑CTCTTTTTCTCCGGTCTGGGGGGCCTTGCGACCGGGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCAttt  <  1:149887/251‑1 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                        GAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTT‑TACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTG‑CTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATC‑GATAAACAGCGTCGAGCCAT‑TGGGGATTTGCTCCGCCACTTTGCGGGCGATGCG‑CTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTTGTTTCATTT  >  NC_000913/3560118‑3560607

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑