Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,509,728	C→T	I33I (ATC→ATT)	yfeO →	putative transport protein YfeO

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,509,728	0	C	T	96.0%	157.2 / ‑3.8	51	I33I (ATC→ATT)	yfeO	putative transport protein YfeO
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (28/20);  total (28/23)
Fisher's exact test for biased strand distribution p-value = 1.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

GGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGC‑TGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTG‑ATTGTGGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCT‑C‑TGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTC  >  NC_000913/2509481‑2509974                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                         gggggCATAACAGCACGGCCGTAAAGCCAGGCCCCCTGTAAAGTTCAATCCTGGTCACAGAATTACTTTACCGTTGCACCACGAATGACAGCTATCCTTTTTTGCA‑GCGTGGGCAAACTAAGTATTTGTCCCCGCATAGGAAATAAACCAGGCTCCACCCGCGGGCCAGTGCCATGTTGTTATTATCGCTTCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTG‑ATTGTGGTGATGAAAATTGcc                                                                                                                                                                                                                                                       <  1:306280/251‑1 (MQ=255)            aacaacGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGC‑TGGGCAAACTAAGTATTTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCAAGAACCATGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTG‑ATTGTGGTGATGAAAATTGCCTCGGTATTAc                                                                                                                                                                                                                                             >  2:272462/1‑250 (MQ=255)             acaacGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCCCCAATTGCAGCGATGCCTTTTTGCAAGC‑TGGGCAAACTAAGTATTTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTG‑ATTGTGGTGATGAAAATTGCCTCGGTATTACAg                                                                                                                                                                                                                                           <  1:443641/251‑1 (MQ=255)                             aggcTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGC‑TGGGCAAACTAAGTATTTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGGTGTTATTATCGCTCCCCGCCGTGGCAATGGGGATTGCGTCCAGTATTATTCTG‑ATTGTGGTGATGAAAATTGCCTCGGTATTACAGACTTTGCTCTGGCAAc                                                                                                                                                                                                                           >  1:117863/1‑251 (MQ=255)                             aggcTCCCTGTAAATATCGATCTGGGTAACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGC‑TGGGCAAACTAAGTATTTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTG‑ATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAAc                                                                                                                                                                                                                           >  1:346906/1‑251 (MQ=255)                                       tAAATATCGATCTGGGTCACACCATTAATTTATCGTTCAAGCACAAATTACAGCGATGCCTTTTTGCAAGC‑TGGGCAAACTCAGTATTTGACCCAGCATAAGGACTAGAACATGCTCCATCAGCGAGCCAGAACCATGTTGTTATTATCGCTCCACGCTGTCGCCATTGGGATAGCGTCCAGTCTTATTCTG‑ATTGTGGTGATGAAAATTGCGTCGGTATTCCAGAATTTGCTCTGGCAACGACTGCCggg                                                                                                                                                                                                                 >  2:440229/1‑251 (MQ=255)                                         aataccgattggggTCACACATTTACTTTACCGTTTCAGCACCAATTGCAGCTATGCCTTTTTGAAAGG‑TGGGAAAATTAATTATTTGCCCCCGGAAAAGGAATAGAACATGCTCCACCCGCGAGCCAGAACCATGTTGTTATTATCCCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTG‑ATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGaa                                                                                                                                                                                                               <  2:346906/240‑1 (MQ=255)                                                   tGGGTCACACATTTACTTTATCGTTTCAGCCCCAATTGCAGCGATGCCTTTTTGCAAGC‑TGGGCAAACTAAGTATTTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTG‑ATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGa                                                                                                                                                                                                      <  1:377124/250‑1 (MQ=255)                                                    gggTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGC‑TGGGCAAACTAAGTATTTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCAAGAACCATGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTG‑ATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGAt                                                                                                                                                                                                     >  1:263663/1‑250 (MQ=255)                                                              attattTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGC‑TGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACATGCTCCAACCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTG‑ATTGTGGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGAtt                                                                                                                                                                                          <  2:427680/251‑1 (MQ=255)                                                              aTTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGC‑TGGGCAAACTAAGTATTTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTG‑ATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGAtt                                                                                                                                                                                          >  2:325017/1‑251 (MQ=255)                                                               ttACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGC‑TGGGCAAACTAAGTATTTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTG‑ATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAAATCTGGGGATAGCCCAGGATTc                                                                                                                                                                                         >  1:295574/1‑251 (MQ=255)                                                                     gaTCGTTTCAGCCCAAATTGCAGCGATGCTTTTTTGCAAGC‑TGGGCAAATTAAGTATTTGACCCACCAAAAGGAATAAAACATGCTCCATCCGCGACCTAGAACCATGTTTTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTG‑ATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGAGTCAcccc                                                                                                                                                                                    <  1:335660/249‑1 (MQ=255)                                                                       acgtgtcaacACCAATTGCTGCGATTCCTTTTTGCAGGC‑TGGGCAAACTAAGTATTTGACCCCGCATAAGGAATAGAGCTTGCTCCACCCCAGAGCCAGAACCGTTTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTG‑ATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGGCTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCct‑c‑t                                                                                                                                                                               <  1:383688/242‑1 (MQ=255)                                                                           ttCAGCACCAATTGCAGCGATGCCTTTTTGCAAGC‑TGGGCAAACTAAGTATTTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTG‑ATTGTGGTGATGAAAATTGCCTCGGAATTACAGAATTTGCTCTGG‑AACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCT‑C‑TGGatc                                                                                                                                                                          >  2:239465/1‑251 (MQ=255)                                                                               gcacaaagtgacgcattGCCTTTTTGCAAGC‑GGGCCAAACTAAGTATTTGACCCCGAATAAGTAAAAGAACATGCTCCATCCGCGAGCCAGAACTATGTTGTTATTATCCCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCGTATTCTG‑ATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCT‑C‑TGGatcatc                                                                                                                                                                       <  2:353493/235‑1 (MQ=255)                                                                                    aaTTGCAGCGATGCCTTTTTGCAAGC‑TGGGCAAACTAAGAATTTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTG‑ATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCT‑C‑TGGATCATCGgtgt                                                                                                                                                                  >  1:234885/1‑251 (MQ=255)                                                                                    aaTTGCAGCGAGGCATTTTTGCAAGC‑TGGGCAAACTAAGTATTTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGATGTTATTATCGATCCACGACGTGGCAATTGGGATTGCGTACAGTCTTATTCTG‑ATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCT‑C‑TGGATCATCGgtgt                                                                                                                                                                  >  2:192960/1‑251 (MQ=255)                                                                                     aTTGCAGCGATGCCTTTTTGCAAGC‑TGGGCAAACTAAGTATTTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTG‑ATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCT‑C‑TGGATCATCGgtgt                                                                                                                                                                  >  1:368988/1‑250 (MQ=255)                                                                                         cagcaatcccTTTTTGCAATC‑TGGGCAAACTAAGTATTTGACCCCGCATAAGGAATAGAACATGCCCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTG‑ATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCT‑C‑TGGATCATCGGTGTATTa                                                                                                                                                              <  1:136902/242‑1 (MQ=255)                                                                                            cGATGCCTTTTTGCAAGC‑TGGGCAAACTAAGTATTTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGAGTCCAGTCTTATTCTG‑AATGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCT‑C‑TGGATCATCGGTGTATTAACGc                                                                                                                                                          >  1:50871/1‑251 (MQ=255)                                                                                                                gggCAACCTAAGTATTTGCCCCCACTTAAGGAATAGAACATCCTCAATCCGCGAGCCAGAACCGTGTTTTTATTACCCCCCCCCGCCGTGGCAAGTGGGATTGCGTCCAGTCTTATTCTG‑AATGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCT‑C‑TTGATCGTCGGTGTATTAACGCTAACGGGTATTGCGGTggg                                                                                                                                       <  2:50871/251‑1 (MQ=255)                                                                                                                  gCAAACTAAGTATTTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTG‑ATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCT‑C‑TGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGggt                                                                                                                                     >  1:96845/1‑251 (MQ=255)                                                                                                                     aacaaaGTATTTGACCCCGCATAAGAAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTTTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTG‑ATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCT‑C‑TGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGggttgg                                                                                                                                  <  1:53733/247‑1 (MQ=255)                                                                                                                     aaCTAAGTATTTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGTGCCAATTGGGATTGCGTCCAGTCTTATTCTG‑ATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGCACTCTGGGGATAGCCCAGGATTCACCCCT‑C‑TGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGggttgg                                                                                                                                  >  2:374610/1‑251 (MQ=255)                                                                                                                         aaGTATTTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTG‑ATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCT‑C‑TGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTAt                                                                                                                              >  2:390415/1‑251 (MQ=255)                                                                                                                          aGTATTTGACCCCGCATAAGGAATAAAACATGCTCAATCCGCGAGCCAGAACAATTTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTG‑ATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAATGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCT‑C‑TGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATc                                                                                                                             <  1:439850/251‑1 (MQ=255)                                                                                                                                tGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTG‑ATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCT‑C‑TGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTc                                                                                                                       >  2:203765/1‑251 (MQ=255)                                                                                                                                     ccGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTG‑ATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCT‑C‑TGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGcc                                                                                                                   >  1:23187/1‑250 (MQ=255)                                                                                                                                           aaGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTG‑ATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCT‑C‑TGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTc                                                                                                             >  1:196483/1‑250 (MQ=255)                                                                                                                                              gAATAGAACATGCTCCGTCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTG‑ATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCT‑C‑TGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGc                                                                                                         >  1:238838/1‑251 (MQ=255)                                                                                                                                                aTAGACCATCCTCCATCCGCGAGCAAGAACCATGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTG‑ATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCT‑C‑TGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCg                                                                                                       <  2:161408/251‑1 (MQ=255)                                                                                                                                                   gAACATGATCCATCCGCGAGCCAAAACCATGTTGTAATTATCGCTCCCCGCCGTAGCAATTGGGATTGCGTCCAGCCTAATTCTG‑ATTGTGGAGATGAAAATTGTCTCGGTATTACAGAATTTGCTCTGGCCACGACACCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCT‑C‑TGGATCATCGGTGTATTAAAGCTAACGGGTATTGCGGCAGGGTTGGATCTACGTTTCAGCAAGGTTCATGCCNGAc                                                                                                    >  2:406792/1‑251 (MQ=255)                                                                                                                                                          cTCCATCCGCGAGCCAGAAGCATGTTGTTATTATAGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTG‑ATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCT‑C‑TGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGAccc                                                                                             >  1:392646/1‑251 (MQ=255)                                                                                                                                                                                gtttttttttcccgcccccggccgTGGCAATTGGAATTCCGTCAAGTCTTATTCGG‑ATTGTGGTGATAAAAATTGCATGGTTATTACAGAATTTGCTGTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGAT‑CCCCCCTGC‑TGAATCAGCGGTGTATTAACGCAAACGGTTTTGGCGGTGGGGTTGGTTACCCTTTTCAGCCAGGGTCATGCCGGACCATGCCCCGCCTGTGAACCGCTGATCGgc                                                                       <  2:263663/231‑1 (MQ=255)                                                                                                                                                                                    gttATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTG‑ATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCT‑C‑TGGATCATCGgtgt                                                                                                                                                                  >  2:40973/1‑156 (MQ=255)                                                                                                                                                                                    gttATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTG‑ATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCT‑C‑TGGATCATCGgtgt                                                                                                                                                                  <  1:40973/156‑1 (MQ=255)                                                                                                                                                                                      tattacccccccccGCCGGGGCAATTGGGATTGCGTCCAGTCTTATTCTG‑ATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCGGGCAACGACTGCCGGGAACTCTGGGGATAGCCCGGGATTCCCCCCT‑C‑TGGATCATCGGTGTATTAACGCTACGGGGTATTGCGTTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGCCCAGCCCCCGCCTGTGAACCGCTGATCGGCGCAcc                                                                  <  1:306607/240‑1 (MQ=255)                                                                                                                                                                                          tcccccccccGCCGTGGCAATGGGGATTGCGTCCAGTCTTATTCTG‑ATTGTGGTGATGAAAATTGCCTCGGTATTAAAGAATTTCCTCTGGCAACGACTGCGGGGAACTCTGGGGATAGCCCGGGATTCACCCCT‑C‑TGGATCATCGGTGTATTAACGCTAAGGGGTATGGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGtt                                                              <  1:169052/244‑1 (MQ=255)                                                                                                                                                                                             gCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTG‑ATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCT‑C‑TGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGc                                                          >  1:395787/1‑251 (MQ=255)                                                                                                                                                                                             cccccccGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTG‑ATTGTGGTGAGGAAAATTGCCTCGGTATTACAGAATTTGCTCGGGCACCGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCG‑C‑TGGATCATCGTTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTTTCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGc                                                          <  1:21085/248‑1 (MQ=255)                                                                                                                                                                                                    gcgggggaaatggggtttGCGTCAAGTCTT‑TTCTGGATGGGGGTGATGAAAATGGCCTCGGTATTACAAAATTTGCTCTGGCACCGACTGCCGGGAACTCTGGGGATA‑CCCCGGAATTACCCCC‑CGTGGATCATCGGTGTTTTACCGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCAGGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTgc                                                   <  2:9193/235‑1 (MQ=255)                                                                                                                                                                                                         ggCAATTGGGATTGCGTCCATTCTTATTCTG‑ATGGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCACCGACTGCCGGGAACTCTGGGGAGAGCCCGGGATCCACCCTT‑C‑TGGATCATCGGTGTATTAACGCAAACGGGTATTGCGGTGGGGTTGGTTTCCCGTTTCACCCGGGGTCATGCCGGCCCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTAc                                              <  2:6649/251‑1 (MQ=255)                                                                                                                                                                                                          gCAATTGGGATTGCGTCCAGTCTTATTCTG‑ATTGTGGTGAGGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCT‑C‑TGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGATGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTAcc                                             >  2:168966/1‑251 (MQ=255)                                                                                                                                                                                                                       cGTCCAGTCTTATTCTG‑ATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCT‑C‑TGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGt                                 >  1:411515/1‑250 (MQ=255)                                                                                                                                                                                                                            attCTTATTCAG‑ATTGGGGTAATGAAAATTGCCTCGAAATTACAGAATTTGCTCTGGCAACGACTGCGGGGAACTCTGGGGAGACCCCAGGATCCACCCCT‑C‑GGGATCATCGGTGTATTACCTCAGACGGGTGTTGCGGGGGGGTTGTTTACCCGTTTCAGCCAGGGCCAGGCCGGCCCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCAtt                           <  2:168478/249‑1 (MQ=255)                                                                                                                                                                                                                            aGTCTTATTCTG‑ATTGTGGTGATGAAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCT‑C‑TGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGNAtt                           >  2:417255/1‑251 (MQ=255)                                                                                                                                                                                                                                 tATTCTG‑ATTGTGGTGCTGTAATTC‑CCTCGTTATTCCATAATTTGCTCTGGCACCCACTGCCGGGAACTCTGGGGATAGCCCAGTATTCCCCCCT‑C‑TGGATCATCGGTTTATTAACGCTAACGGTTTTTGCGGTGGGTTTTTTTTTCCTTTTCACCCAGGTCCCTGCCGCACCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAAtt                       <  1:900/249‑1 (MQ=255)                                                                                                                                                                                                                                                 gatgaAAATTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCT‑C‑TGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTggcg       >  1:342033/1‑251 (MQ=255)                                                                                                                                                                                                                                                   agaAAATTGCCTCGGTATTACAAAATTTGCTCGGGCAACGCCTGCCGGGAACTCTGGGGATAGCCCAGGTTTCACCCCT‑C‑GGGATCATGGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTACCCTTTTCCGCCAGGGTCATCCCGGCCCAGCCCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTggcggc     <  2:356814/250‑1 (MQ=255)                                                                                                                                                                                                                                                      aaaTTGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCT‑C‑TGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTc  >  1:160787/1‑251 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                         GGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGC‑TGGGCAAACTAAGTATCTGACCCCGCATAAGGAATAGAACATGCTCCATCCGCGAGCCAGAACCATGTTGTTATTATCGCTCCCCGCCGTGGCAATTGGGATTGCGTCCAGTCTTATTCTG‑ATTGTGGTGATGAAAATCGCCTCGGTATTACAGAATTTGCTCTGGCAACGACTGCCGGGAACTCTGGGGATAGCCCAGGATTCACCCCT‑C‑TGGATCATCGGTGTATTAACGCTAACGGGTATTGCGGTGGGGTTGGTTATCCGTTTCAGCCAGGGTCATGCCGGACCAGACCCCGCCTGTGAACCGCTGATCGGCGCACCGGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTC  >  NC_000913/2509481‑2509974

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