Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,384,808	C→T	K305K (AAG→AAA)	yfbK ←	IPR002035/DUF3520 domain‑containing protein YfbK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,384,808	0	C	T	93.5%	192.9 / ‑5.3	62	K305K (AAG→AAA)	yfbK	IPR002035/DUF3520 domain‑containing protein YfbK
Reads supporting (aligned to +/- strand):  ref base C (0/4);  new base T (35/23);  total (35/27)
Fisher's exact test for biased strand distribution p-value = 3.15e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATACTTTCTGCGCTTCAGAGAGGGTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGT‑TGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACC‑AGACTCCCGCTGTTTTTTGACCATTG‑ATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGA‑TTGATGCCGCCCTTAATAAACCCCTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAA‑TCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTTCATCAGAAATCATTGAACCAGA  >  NC_000913/2384566‑2385053                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                          aaTACTTTCTGCGCTTCAGAGAGGGTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGT‑TGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACC‑AGACTCCCGCTGTTTTTTGACCATTG‑ATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGA‑TTGATGACGCCCTTAATAAACCCTTT‑C‑GTCGc                                                                                                                                                                                                                                                 >  1:568119/1‑251 (MQ=255)       ttCTGCGCTTCAGAGAGGGTATCAATGTAGCTGTAGTTGCCGTTACCAACAACGGCAATTCGCACCATCATTGCCTCGT‑TGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACC‑AGACTCCCGCTGTTTTTTGACCATTG‑ATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGA‑TTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCctgct                                                                                                                                                                                                                                            <  2:117455/250‑1 (MQ=255)         cTGCGATTCAGACAGGATATCAATGTCGCTGTAGTTGCCGTTACAAACATCGGCCATTCGCACCATCATTGCCTCGA‑TGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACC‑AGACTCCCGCTGTTTTTTGACCATTG‑ATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGA‑TTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCTTGCTGa                                                                                                                                                                                                                                          >  2:328154/1‑250 (MQ=255)            cgcTTCAGAGAGGGTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGT‑TGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACC‑AGACTCCCGCTGTTTTTTGACCATTG‑ATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGA‑TTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAg                                                                                                                                                                                                                                      <  1:424434/251‑1 (MQ=255)                   gagagGGTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGT‑TGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACC‑AGACTCCCGCTGTTTTTTGACCATTG‑ATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGA‑TTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTTCTGATAAGCC‑AGTTc                                                                                                                                                                                                                              <  1:542699/251‑1 (MQ=255)                         gTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGT‑TGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACC‑AGACTCCCGCTGTTTTTTGACCATTG‑ATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGA‑TTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTCCAGccc                                                                                                                                                                                                                        <  1:598996/251‑1 (MQ=255)                           aTCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGT‑TGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACC‑AGACTCCCGCTGTTTTTTGACCATTG‑ATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGA‑TTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTCCAGCCCgg                                                                                                                                                                                                                      >  1:623867/1‑251 (MQ=255)                               aTGTAGCTGTAGTGGCCGTTACCAACATCGTCATTTCGCACCATCATTGCCTCGT‑TGTAATTGCTATTCCCCACGCCAAACGTCAACAGAGTAACACC‑AGACTCCCGCTGTTTTTTGACCATTG‑ATTCAATCGTTTTGGGATCGTCAATGCCAACGTTAAAGTCACCGCCAGTGGCTAATAAAATGCGA‑TTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTCCAGCCCGGCAcc                                                                                                                                                                                                                  <  1:495546/251‑1 (MQ=255)                                  tAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGT‑TGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACC‑AGACTCCCGCTGTTTTTTGACCATTG‑ATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGA‑TTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTCCAGCCCGGCACCGcc                                                                                                                                                                                                               >  1:322845/1‑251 (MQ=255)                                                                aTTCGCACCATCATTGCCTCGT‑TGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACC‑AGACTCCCGCTGTTTTTTGACCATTG‑ATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGA‑TTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAA‑Tc                                                                                                                                                                                >  1:251866/1‑251 (MQ=255)                                                                       ccATCATTGCCTCGT‑TGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACC‑AGACTCCCGCTGTTTTTTGACCATTG‑ATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGA‑TTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTCCAGCCCGGCACCGCCATTGGTACTTCCTTCGGCATCCAGCGAA‑TCAATTGCg                                                                                                                                                                         >  2:689400/1‑251 (MQ=255)                                                                        aatcatTGCCT‑GTGTGTAATTTCTATCCCCCACCCCAAACTTCGACAGAGTAACCCC‑AGACTCCCGCTGTTTTTTGACCATTG‑ATTCAAGCGATTTGGGAGCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAAAAAAAGGCGA‑TGGAGGCCGCCCTAAATAACCCCTTT‑C‑GCCGCCTGCTAATAACCA‑AGTCCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAA‑TCAATTGCgg                                                                                                                                                                        <  2:322845/250‑1 (MQ=255)                                                                            tgtGCCTTGT‑TTTAATTGCCTTCCCCCACGCCAAACGTCAACAGAG‑AACACCCAGCCCCCCGCTTTTTTTTG‑CCCTTGTATTCATTGGATTTTGGGTCGGCAATGCCAACTTTAAATTCTCTGGCAGTGGCTAATAAAATTCGA‑TGGAGCCCGCCCTTAATAA‑CCCCTTGC‑TTGGCGTGTTTATAGGCC‑AGTTCCAGCCCGGCGCCTCCATTGGTACTGCCTTCGGCATCCAGCGAA‑TCAATTGCGGCAt                                                                                                                                                                     <  2:484033/248‑1 (MQ=255)                                                                                 cTCGT‑TGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACC‑AGACTCCCGCTGTTTTTTGACCATTG‑ATTCAATCGATTTTGGATCGTCAATGCCCACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGA‑TTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAA‑TCAATTGCGGCATTAAttt                                                                                                                                                               >  1:409968/1‑251 (MQ=255)                                                                                    gT‑TGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACC‑AGACTCCCGCTGTTTTTTGACCATTG‑ATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGA‑TTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAA‑TCAATTGCGGCATTAATTTCCg                                                                                                                                                            <  1:689400/251‑1 (MQ=255)                                                                                         tAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACC‑AGACTCCCGCTGTTTTTTGACCATTG‑ATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGA‑TTGATGCCGCCCTTAATAAACCATTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAA‑TCAATTGCGGCATTAATTTCCGCCtt                                                                                                                                                        >  1:491967/1‑251 (MQ=255)                                                                                          aaTTGCTATTCCCCACGCCAAACGTCGCCAGAGTAACACC‑AGACTCCCGCTGTTTTTTGACCATTG‑ATTCAATCGATTTTGGATCGTCAATGCCACCGTTAACGTCACCGTCAGTGGCTAATAAAATGCGA‑TTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAA‑TCAATTGCGGCATTAATTTCCGCCtt                                                                                                                                                        >  1:78498/1‑250 (MQ=255)                                                                                              gCTATTCCCCCCGCCAAACGTCGACAGAGTACCACC‑AGTCCCCCGTTGTTTTTTGACCATTG‑ATTAACTCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCATTGTCTAATAAAATGCGA‑TTAATGCCGCCCTTAATAACCCCTTT‑C‑TTCGCCTGCTGATAACCC‑AGTCCCATCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAA‑TCAATTGCGGCATTAATTTCCGCCTTATGa                                                                                                                                                    <  2:255551/250‑1 (MQ=255)                                                                                                         cGCCAAACGTCGACAGAGTAACACC‑AGACTCCCGCTGTTTTTTGACCATTG‑ATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGA‑TTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTCCAGCCCGGCAc                                                                                                                                                                                                                   <  2:235578/177‑1 (MQ=255)                                                                                                         cGCCAAACGTCGACAGAGTAACACC‑AGACTCCCGCTGTTTTTTGACCATTG‑ATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGA‑TTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTCCAGCCCGGCAc                                                                                                                                                                                                                   >  1:235578/1‑177 (MQ=255)                                                                                                           ccAAACGTCGACAGAGTACCACC‑AGACTCCCGCTGTTTTTTGACCATTG‑ATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGA‑TTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAA‑TCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGAt                                                                                                                                         >  1:224414/1‑248 (MQ=255)                                                                                                                cGTCGACAGAGTAACACC‑CGACTCCCGCTGTTTTTTGACCATTG‑ATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGC‑TTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAA‑TCAATTGCGGCCTTAATTTCCGCCTTATGACTCCCGGAGATAGAAGgc                                                                                                                                  >  2:330899/1‑250 (MQ=255)                                                                                                                     aCAGAGTAACACC‑AGACTCCCGCTGTTTTTTGACCATTG‑ATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGA‑TTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGAGAA‑TCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATCGAAGgcaatgc                                                                                                                             >  1:654125/1‑250 (MQ=255)                                                                                                                      cAGAGTAACACC‑AGACTCCCGCTGTTTTTTGACCATTG‑ATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGA‑TTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCC‑Ag                                                                                                                                                                                                                                 >  1:7016/1‑150 (MQ=255)                                                                                                                      cAGAGTAACACC‑AGACTCCCGCTGTTTTTTGACCATTG‑ATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGA‑TTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCC‑Ag                                                                                                                                                                                                                                 <  2:7016/150‑1 (MQ=255)                                                                                                                      cAGAGTAACACC‑AGACTCCCGCTGTTTTTTGACCATTG‑ATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGACTAATAAAATGCGA‑TTGATGCCGCCCTTAATAAACCCTTT‑C‑GACGCCTGCTGATAAGCC‑AGTTCCAGACCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAA‑TCAATTGCGGCATTCATTTCCGCCTTATGACTCCCGGAGATAGAAGgcaatgct                                                                                                                            >  2:126261/1‑249 (MQ=255)                                                                                                                       agagGAACACC‑AGACTCCCGCTGTTTTTTGACCATTG‑ATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGA‑TTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTCCAGCCCGGCCCCGCCATTGGTACTGCCTTCGGCATCCAGCGAA‑TCAATTGCGGCCTTAATTTCCGCCTTATGACTCCCGGAGATCGAAGgcaatgcaat                                                                                                                          >  2:231770/1‑251 (MQ=255)                                                                                                                        gagTAACCCC‑AGACTCCCGCTGTTTTTTGACCATTG‑ATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGA‑TTGATGCCGCCCTTAATAACCCCTTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAA‑TCAATTGCGGCATTa                                                                                                                                                                   <  1:523901/209‑1 (MQ=255)                                                                                                                        gagTAACCCC‑AGACTCCCGCTGTTTTTTGACCATTG‑ATTCAATCGATTTTGGATCGTCAATGCAAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGA‑TTGATGCCGCCCTTAATAACCCCTTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAA‑TCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATa                                                                                                                         <  1:545873/251‑1 (MQ=255)                                                                                                                        gagTAACACC‑AGACTCCCGCTGTTTTTTGACCATTG‑ATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGA‑TTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAA‑TCAATTGCGGCATTa                                                                                                                                                                   >  2:523901/1‑209 (MQ=255)                                                                                                                        aggtgaactccacaCTCCCGTTGTTTTTTGTCCATAT‑ATTAAATGGATTTGGGTGCGTCAATGCCACCTTGAAAGCCACCGTCAGGGGCAATTAAAATGCTA‑TTGCGGCCGACCTGATTAAACCCTTT‑G‑GTCACGTCCTGATAGGCG‑AGTTCCATCCCGGCTCCGCCTTTGGTACTGCCGTCGGCATCCA‑CGAAATCAAATGCGGCATTAATGTCCGCCTTATGACTCCCGGAGATAGAAGgcgatgcaat                                                                                                                          <  1:420966/238‑1 (MQ=255)                                                                                                                             acccc‑AGACTCCCGCTGTTTTTTGACCATTG‑ATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGA‑TTGAGGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAA‑TCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGa                                                                                                                    <  1:99958/248‑1 (MQ=255)                                                                                                                                    aaCTCCCTCTTTTTTTTGACCATTG‑ATTTAACCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCTTCAGTGGCTAATAAAAGGCGA‑TTGATGCCGACCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAA‑TCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGcc                                                                                                              <  1:633105/250‑1 (MQ=255)                                                                                                                                         ccGCTGTTTTTTGACCATTG‑ATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGA‑TTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAA‑TCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCg                                                                                                          >  2:271222/1‑250 (MQ=255)                                                                                                                                                tttttGACCATTG‑ATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGA‑TTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTCCAGCCCGGCACCGCCATTGGCACTGCCTTCGGCATCCAGCGAA‑TCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCAc                                                                                                  >  2:121773/1‑251 (MQ=255)                                                                                                                                                      cacctttgatcaaaTCGATTTTGGATCGACAATCCCAACGTCAAAGTCACCTTCATTGGCTAAAAAAAGTCTA‑TTGAGCCCCCCCTTAATAA‑CCCCTT‑CCGTCCCCTGCTGATAAGCC‑AGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAA‑TCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATg                                                                                              <  2:623867/238‑1 (MQ=255)                                                                                                                                                      aCCATTG‑ATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGA‑TTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAA‑TCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCCCGATGGc                                                                                            >  2:330340/1‑251 (MQ=255)                                                                                                                                                       ccATTG‑ATTCAATCTAGTTTGGTTCTTCAATCCCAACGTTAAATTCACCGTCGGTGCTTAATAAAAGGCGA‑TTGATGCCGCCCTTAATAACCCCTTT‑C‑GTCGCTTGCTGA‑AAGCCCAGTTCCAGCCCGGCTCCGCCATTGGTACTGCCTTCGGCATCCAGCGAA‑TCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCa                                                                                           <  2:688668/251‑1 (MQ=255)                                                                                                                                                         aTTG‑ATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGA‑TTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAA‑TCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAAt                                                                                         >  1:563469/1‑251 (MQ=255)                                                                                                                                                          ttG‑ATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGA‑TTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAA‑TCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATg                                                                                        >  1:698655/1‑251 (MQ=255)                                                                                                                                                               ttCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGA‑TTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATc                                                                                                                                                                                          >  2:512733/1‑150 (MQ=255)                                                                                                                                                               ttCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGA‑TTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATc                                                                                                                                                                                          <  1:512733/150‑1 (MQ=255)                                                                                                                                                                  aagcGATTTTGTATCGTCAATCCCAACGTAAATGCCTCCGTAGGTGGTTATAAAAATG‑GAGTGGATGCCTCCCTTAA‑AAACCCCTTGC‑GTCGCCTGCTGATAAGCC‑GGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAA‑TCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTAGCCAGCGTAGGTCACGATGGCAATGTTATCCt                                                                                 <  2:289199/248‑1 (MQ=255)                                                                                                                                                                     cGATTTTGGATCGTCAATGCCAACGTTAAAGCCACCTTCAGTGGCTAATAAAATGCGA‑TTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAA‑TCAATTGCGGCATTAATGTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGc                                                                               <  1:121773/250‑1 (MQ=255)                                                                                                                                                                         tttGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGA‑TTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAA‑TCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATCCGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTAACg                                                                          >  2:439632/1‑251 (MQ=255)                                                                                                                                                                          ttGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGA‑TTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAA‑TCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACg                                                                          >  1:42995/1‑250 (MQ=255)                                                                                                                                                                          ttGGATCGTCAATGCCAACGTTAAAGTCAACGTCAGTGGCTAATAAAATGCGA‑TTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAA‑TCAATTGAGGCATTAATTTCCGCCTTATGACTCCCGGAGCTAGAAGTCAATGCAATAAGGGAGTCGCAAGCGTAGGACACGATGGCAATGTTATCCTGccaacgg                                                                         >  2:213430/1‑245 (MQ=255)                                                                                                                                                                          ttGGATCGTCAATGCAAACGTTAAAGCCACCGTCAGTGGCTAATAAAATGCGA‑TTGATGCCGCCCTTAA‑AAACCCCTTTC‑GTCGCCTGCTGATAAGCC‑AGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAA‑TCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACg                                                                          <  1:363342/250‑1 (MQ=255)                                                                                                                                                                            ggATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGA‑TTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAA‑TCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGa                                                                         >  1:111098/1‑249 (MQ=255)                                                                                                                                                                                    aaggCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGA‑TTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAA‑TCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGa                                                               <  1:35608/248‑1 (MQ=255)                                                                                                                                                                                     aTGCCAACGTTAAAGTCACCGTCAGTGGATAATAAAATGCGA‑TTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAA‑TCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGAc                                                              <  1:277380/251‑1 (MQ=255)                                                                                                                                                                                            cGTTAAAGTCACCGTCAGTGGCTAATAAACTGCGA‑TTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAA‑TCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCaataat                                                       >  1:297555/1‑251 (MQ=255)                                                                                                                                                                                                 aaGTCACCGAGAGTTGCTAATAAAGTCCGA‑TTGATGCCGCCCTTGCGACACACTTT‑C‑GACGCCTGCTGATCAGCC‑CCTTCCAGACCTGCACCGCCCTTGGTACGGCCTTCGGAATCCAGCGAA‑TCACTTGCGGCATTAATTTCCGCCTTACGCCGCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGACTACGTAAGACACGATATCAGTGTTATCCTGCTCCCACAGTTCTTTGCCCAATACTTTCa                                                   >  2:206866/1‑250 (MQ=255)                                                                                                                                                                                                      aCCGTCAGTGGCTAATAAAATGCGA‑TTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAA‑TCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGa                                               >  1:265772/1‑249 (MQ=255)                                                                                                                                                                                                                 cTAATAAAATGCGA‑TTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAA‑TCAATTGCGGCATTAATTTCCGCCTTATGACACCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGAGCACGATGGCAATGATATCCTGCTCACGAAGTTCTTTGACCAATAATTTAAACGAAGAAGGGATAAgg                                  >  1:387512/1‑250 (MQ=255)                                                                                                                                                                                                                  taataaAATGCGA‑TTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAA‑TCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTg                                 >  1:118652/1‑251 (MQ=255)                                                                                                                                                                                                                            cGA‑TTGATGCCGCCCTTAATAAACCCTTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTCCAGCCCGGCACCGCCATTGGTAATGCCTTCGGCATCCAGCGAA‑TCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTtc                       >  1:433951/1‑251 (MQ=255)                                                                                                                                                                                                                                         cccTTAATACACCCTTT‑C‑GTCGCCTTCTTATAAGCC‑TTTCCCAGCC‑GGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAA‑TCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTTCATCAGAAATCAtt          <  1:94941/251‑1 (MQ=255)                                                                                                                                                                                                                                             taataaACCCTTG‑C‑GCCGCCTGCTGATAAGCC‑AGTTCCAGCCCGGCCCCGCCATTGGTACTGCCTTCGGCATCCAGCGAA‑TCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTTCATCAGAAATCATTg         <  2:265772/249‑1 (MQ=255)                                                                                                                                                                                                                                                 aaaCCCTTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTACAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCATCGAA‑TCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGCGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGAAATGTTATCCTGCTCACGAAGTTCTTTAACCAATAATTTCAACGAAGCCTGGATAAGTGGCGAACGTTCATAAGAAATCATTGAAca     >  2:472138/1‑248 (MQ=255)                                                                                                                                                                                                                                                  caCCCTTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAA‑TCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGCGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGCCTGGATAAGTGGCAAACGTTCATCAGCAAGCATTGACCCaa   >  2:580677/2‑249 (MQ=255)                                                                                                                                                                                                                                                  aaCCCTTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAA‑TCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTTCATCAGAAATCATTGAACCAGa  <  1:408977/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                          AATACTTTCTGCGCTTCAGAGAGGGTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGT‑TGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACC‑AGACTCCCGCTGTTTTTTGACCATTG‑ATTCAATCGATTTTGGATCGTCAATGCCAACGTTAAAGTCACCGTCAGTGGCTAATAAAATGCGA‑TTGATGCCGCCCTTAATAAACCCCTT‑C‑GTCGCCTGCTGATAAGCC‑AGTTCCAGCCCGGCACCGCCATTGGTACTGCCTTCGGCATCCAGCGAA‑TCAATTGCGGCATTAATTTCCGCCTTATGACTCCCGGAGATAGAAGGCAATGCAATACGGGAGTCGCCAGCGTAGGTCACGATGGCAATGTTATCCTGCTCACGAAGTTCTTTGACCAATAATTTCAACGAAGACTGGATAAGTGGCAAACGTTCATCAGAAATCATTGAACCAGA  >  NC_000913/2384566‑2385053

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