Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,085,319	G→A	A69A (GCC→GCT)	tsuA ←	thiosulfate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,085,319	0	G	A	98.5%	234.0 / NA	65	A69A (GCC→GCT)	tsuA	thiosulfate transporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (32/32);  minor base C (1/0);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAAATTAAAAGTCTCAGCAATAGAGTTATGTTCAGTACTGTAGTGCTGCAAGGTTTGATTTAAA‑CCACTGGCATGTGGAGAACGCATCACCGCACTCATCACCATATAAGTGAAAA‑GCGCGATCCAACTGCCGATCAATCCCTCACCCGCGCGATACCAGGTCCCGGTGGCACA‑TCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAATGTTAATGAGG  >  NC_000913/2085072‑2085556                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                              aCAAATTAAAAGTCTCAGCAATAGAGTTATGTTCAGTACTGTAGTGCTGCAAGGTTTGATTTAAA‑CCACTGGCATGTGGAGAACGCATCACCGCACTCATCACCATATAAGTGAAAA‑GCGCGATCCAACTGCCGATCAATCCCTCACCCGCGCGATACCAGGTCCCGGTGGCACA‑TCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCt                                                                                                                                                                                                                                            >  1:1208869/1‑251 (MQ=255)               ccAGCAATAGAGTTATTTTCAGTACTGTAGTGCTGCAAGGTTTGATTTAAA‑CCACTGGCATGTGGAGAACGCATCACCGCACTCATCACCATATAAGTGAAAA‑GCGCGATCCAACTGCCGATCAATCCCCCACCCGCGCGATCCCAGGTCCCGGTGGCCCA‑TCCCCCCGCCAGAACAATTCCCAGCCCGAAGATATACCCCCCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAg                                                                                                                                                                                                                              <  1:833997/250‑1 (MQ=255)                aaacaaaaggtttATGTTCAGTACTGTAGTGCTGCAAGGTTTGATTTAAA‑CCCCTGGCTGGGGGAGAACCCATCCCCCCACTCACCCCCAAATAAGTGAAAA‑GCGCGATCCACCTCCCGATCAACCCCCCCCCCGCGCGATACCAGGTCCCGGTGGCACA‑TCCCCCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAgg                                                                                                                                                                                                                             <  2:582366/240‑1 (MQ=255)                aaGCAATAGAGTTAGGTTCAGTACTGTATTGCGGCAAGGTTTGATTTAAC‑CCACGGGCGTGTGGAGAACGCACCCCCGCCCTCATCCCCATATAAGTAAAAA‑GCGCGATCCACCTGCCGATAAACCCCTCCCCCGCGCGATACCAGGTCCCGGTGGCACA‑TCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAgg                                                                                                                                                                                                                             <  2:957431/250‑1 (MQ=255)                  gCAATAGAGTTATGTTCAGTACTGTAGTGCTGCAAGGTTTGATTTAAA‑CCACTGGCATGTGGAGAACGCATCACCGAACTCATCCCCATATAAGTGAAAA‑GCGCGATCCAACTGCCGATCAATCCCTCACCCGCGAGATACCAGGTCACGGAGGCACA‑TCCACCCGACAGCACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCACCGCTTCGTAAGTCAGTAGGcc                                                                                                                                                                                                                           >  1:1164792/1‑251 (MQ=255)                  gCAATAGAGTTATGTTCAGTACTGTAGTGCTGCAAGGTTTGATTAAAA‑CCACTGGCATGTGGAGAACGCATCCCCGCCCTCATCACCATATAAGGGAAAA‑GCGCGATCCAACTGCCGATCAATCCCCCACCCGCGCGATACCAGGTCCCGGTGGCACA‑TCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGcc                                                                                                                                                                                                                           <  2:335093/251‑1 (MQ=255)                   cAATAGAGTTATGTTCAGTACTGTAGTGCTGCAAGGTTTGATTTAAA‑CCACTGGCATGTGGAGAACGCATCACCGCACTCATCACCATATAAGTGAAAA‑GCGCGATCCAACTGCCGATCAATCCCTCACCCGCGCGATACCAGGTCCCGGTGGCACA‑TCCACCCGCCAGAACACTTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGAAAGTCAGTAGGccc                                                                                                                                                                                                                          >  1:48065/1‑251 (MQ=255)                    aataaatttATTTTCAGTATTATAGTGCTGAAAGGTTTGATTTAAT‑CCCCTGGCATGTGGAGAACGCATCACGGCCCTCTTCCCCATATAAGTGAAAA‑GCGCGATCCAGCTGCCGATCAATCCCTCACCCGCGCGATACCAGGCCCCGGGGGCACA‑TCCCCCCGCCAGAACAATTCCCAGCCCTAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCg                                                                                                                                                                                                                         <  1:178530/244‑1 (MQ=255)                    aaTAGAGTTATGTTCAGTACTGTAGTGCTGCAAGGTTTGATTTAAA‑CCACTGGCATGTGGAGAACGCATCACCGCACTCATCACCATATAAGTGAAAA‑GCGCGATCCAACTGCCGATCAATCCCTCACCCGCGCGATACCAGGTCCCGGTGGCCCA‑TCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCg                                                                                                                                                                                                                         >  1:325710/1‑251 (MQ=255)                       agatttATGTTCTGTACTGTAGTACTGTAAG‑TTTTATTTAAAACCCCTGCCATGTGGACACCCAATCCCCTCACTCATCCCAATATAA‑TGAAAAAGCGCGATCCAACTGCCGATCAATCCCTCACCCGCGCGATCCCAGGTCCCGTTGGCACA‑TCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCt                                                                                                                                                                                                                      <  2:1230819/247‑1 (MQ=255)                          gTTATGTTCAGTACTGTAGTGCTGCAAGGTTTGATTTAAA‑CCACTGGCATGTGGAGAACGCATCACCGCACTCATCACCATATAAGTGAAAA‑GCGCGATCCAACTGCCGATCAATCCCTCACCCGCGCGATACCAGGTCCCGGTGGCACA‑TCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGa                                                                                                                                                                                                                    >  2:23546/1‑250 (MQ=255)                            tATGTTCAGTACTGTAGTGCTGCAAGGTTTGATTTAAA‑CCACTGGCATGTGGAGAACGCATCACCGCACTCATCACCATATAAGTGAAAA‑GCGCGATCCAACTGCCGATCAATCCCTCACCCGCGCGATACCAGGTCCCGGTGGCACA‑TCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAAtt                                                                                                                                                                                                                 >  2:204161/1‑251 (MQ=255)                                    gTACTGTAGTGCTGCAAGGTTTGATTTAAA‑CCACTGGCATGTGGAGAACGCATCACCGCACTCATCACCATATAAGTGAAAA‑GCGCGATCCAACTGCCGATCAATCCCTCACCCGCGCGATACCAGGTCCCGGTGGCACA‑TCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCaa                                                                                                                                                                                                          >  1:345610/1‑250 (MQ=255)                                       cTGTAGTGCTGCAAGGTTTGATTTAAA‑CCACTGGCATGTGGAGAACGCATCACCGCACTCATCACCATATAAGTGAAAA‑GCGCGATCCAACTGCCGATCAATCCCTCACCCGCGCGATACCAGGTCCCGGTGGCACA‑TCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAAc                                                                                                                                                                                                      >  1:1065725/1‑251 (MQ=255)                                         gTAGTGCTGCAAGGTTTGATTTAAA‑CCACTGGCATGTGGAGAACGCATCACCGCACTCATCACCATATAAGCGAAAA‑GCGCGATCCAACTGCCGATCAATCCCTCACCCGCGCGATAACAGGTCCCGGCGGCACA‑TCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGATTTAACGCAAAAAccc                                                                                                                                                                                                    >  2:429652/1‑251 (MQ=255)                                                         tGATTTAAA‑CCACTGGCATGTGCAGAACACATCACCGCCCTCATCACAATAAAACCGAAAT‑GCGCGATCCAACTGCCGATCAATCCCTCACCAGCGCGATACCCGCTCCCTGTAGCACA‑TCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTCTCACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTACGTCAGTAGGCCCGCCTGACTTAACGAAAAAACCCCCACGCTTTGTAc                                                                                                                                                                                       >  2:759294/1‑248 (MQ=255)                                                             ttAAA‑CCACTGGCATGTGGAGAACGCATCACCGCACTCATCACCATATAAGTGAAAA‑GCGCGATCCAACTGCCGATCAATCCCTCACCCGCGCGATACCAGGTCCCGGTGGCACA‑TCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAAtt                                                                                                                                                                                >  1:265990/1‑251 (MQ=255)                                                              taaacccACGGGCATGTGGAGACCGCATCACCCCACTCATCACCATATAAGTAAAAA‑GCGCGATCCAACGGCCGATCAATCCCTCCCCCGCGCGATCCCAGGTCCGGGTGGCACA‑TCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAAt                                                                                                                                                                                 <  1:204161/245‑1 (MQ=255)                                                               aaa‑CCACTGGCATGTGGAGAACGCATCCCCGCACTCAGCACCATATAAGTGAAAA‑GCGCGATCCAACTGCCGATCAACCCCTCCCCCGCGCGATACCAGGTCCCGGTGGCACA‑TCCCCCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTg                                                                                                                                                                               <  1:693578/250‑1 (MQ=255)                                                               aaa‑CCACTGGCATGTGGAGAACGCATCACCGCACTCATCACCATATAAGTGAAAA‑GCGCGATCCAACTGCCGATCAATCCCTCACCCGCGCGATACCAGGTCCCGGTGGCACA‑TCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGc                                                                                                                                                                              >  1:388660/1‑251 (MQ=255)                                                                 a‑CCACTGGCATGTGGAGAACGCATCACCGCACTCATCACCATATAAGTGAAAA‑GCGCGATCCAACTGCCGATCAATCCCTCACCCGCGCGATACCAGGTCCCGGTGGCACA‑TCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAg                                                                                                                                                                                                                                 >  1:125306/1‑198 (MQ=255)                                                                 a‑CCACTGGCATGTGGAGAACGCATCACCCCACTCATCACCATATAAGTGAAAA‑GCGCGATCCAACTGCCGATCAATCCCTCACCCGCGCGATACCAGGTCCCGGTGGCACA‑TCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCa                                                                                                                                                                             >  1:911379/1‑250 (MQ=255)                                                                        ggCATGTGGAGAACGCATCACCGCACTCATCACCATATAAGTGAAAA‑GCGCGATCCAAATGCCGATCAATCCCTCACCCGCGCGATACCAGGTCCCGGTGGCACA‑TCCACCCGACAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGACACGGGAACGCCCAAGCTTCGTAAGACAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCAAGCTTTGTACCGAAATTGCAATCAGc                                                                                                                                                                       >  2:571874/1‑250 (MQ=255)                                                                        ggCATGTGGAGAACGCAGCCCCACACTCATCCCCGTATGGATGAAAA‑GCGCGACCCAACTTCCCAACAATCCCTCCCCCGCCCATTACCGGGTCCCGTTGCCACG‑TCCCCCCTCCATAACAATTCCCAACCCGAAGATATCCCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAACTTCGTAAGTCAGTAGTCCCGCCTTAATTAACGCAAAAACCCCCTCGCTTTGTACCGAAATTGCAATCAg                                                                                                                                                                        <  1:429652/249‑1 (MQ=255)                                                                          cATGTGGAGAACGCATCACCGCACTCATCACCATATAAGTGAAAA‑GCGCGATCCAACTGCCGATCAATCCCTCACCCGCGCGATACCAGGTCCCGGTGGCACA‑TCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCa                                                                                                                                                                      >  1:265132/1‑249 (MQ=255)                                                                                          cACCGCACTCATCACCATATAAGTAAAAA‑GCGCGATCCACCTGCCGATCAATCCCTCACCCGCGCGATACCAGGCCCCGGTGGCACA‑TCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGa                                                                                                                                                     <  1:878774/250‑1 (MQ=255)                                                                                                    ccccccATAAAAGGGAAAA‑GCGCGACCCAACTCCCGATCAACCCCCCCCCCGCGCGATACCGGGTCCCTGTGCCCCA‑TCCACCCGCCAGAACAATTCCCAGCCCGAAGATATCCCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTctc                                                                                                                                          <  2:196082/247‑3 (MQ=255)                                                                                                                gTGAAAA‑GCGCGATCCAACTGCCGATCAATCCCTCACCCGCGCGATACCAGGTCCCGGTGGCACA‑TCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCACGGCGTAAAACATGCGATTATTTTTCACGATATACATATc                                                                                                                              >  1:827003/1‑251 (MQ=255)                                                                                                                     aa‑GCGCGATCCAACTGCCGATCAATCCCTCACCCGCGCGATACCAGGTCCCGGTGGCACA‑TCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGAtata                                                                                                                                    >  1:540478/1‑240 (MQ=255)                                                                                                                        gcgcaccccAACGGCCAATCAATCCCTCACCCGCGCGATACCAGGTCCAGGTGCCACA‑TCCCCCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGaaaa                                                                                                                       <  2:345610/244‑1 (MQ=255)                                                                                                                        ccccgtccaacctcccGATCAATCCCTCCCCCGCTCTATCCCAGGCCCCGGTGGCCCA‑TCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGaaaa                                                                                                                       <  1:857647/237‑1 (MQ=255)                                                                                                                           cgATCAAACGGCCGATCAACCCCACCCCGGCGCGATACCAGGTCCGGGTGGCACA‑TCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAacca                                                                                                                    <  2:325710/251‑1 (MQ=255)                                                                                                                                    ccccgaacacccccccccccGGGCGATACAAGGCCCGGGTGCCAAA‑CCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAgg                                                                                                           <  2:265132/234‑1 (MQ=255)                                                                                                                                     gCCGATCAATCCCTCACCCGCGCTATACCAGCCCCCGGTGGCACA‑TCCCCCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAg                                                                                                            <  2:1208869/249‑1 (MQ=255)                                                                                                                                     ccccatcacccccccACCCCCGCGATACCAGGTCCCGGTGCCACA‑CCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGc                                                                                                          <  2:357470/237‑1 (MQ=255)                                                                                                                                        gATCAATCCCTCACCCGCGCGATACCAGGTCCCGGTGGCACA‑TCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCt                                                                                                                 >  2:21004/1‑241 (MQ=255)                                                                                                                                        gATCAATCCCTCACCCGCGCGATACCAGGTCCCGGTGGCACA‑TCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCaa                                                                                                        >  2:91470/1‑250 (MQ=255)                                                                                                                                        gATCAACCCCTCCCCCGCGCGATACCAGGTCCCGGTGGCCCA‑TCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCt                                                                                                                 <  1:21004/241‑1 (MQ=255)                                                                                                                                             aTCCCTCCCCCGCGCGATACCAGGCCCCGGGGGCACA‑TCCCCCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGc                                                                                                  <  1:1240793/251‑1 (MQ=255)                                                                                                                                                ccccccccGCGCGATACAAGGTCCCG‑TGGCCCCCTCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAAGCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCGGTCAGGCAAAAACGCCCa                                                                                               <  2:827003/246‑1 (MQ=255)                                                                                                                                                         cgcgATACCAGGTCCCGGTGGCACA‑TCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATg                                                                                      >  1:882308/1‑251 (MQ=255)                                                                                                                                                            gATACCAGGTCCCGGTGGCACA‑TCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGAtata                                                                                                                                    >  1:955042/1‑202 (MQ=255)                                                                                                                                                            gATACCACGTCCCGGTGGCACA‑TCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGACAAAAACGCCCACGCTGCATGAc                                                                                    >  2:689523/1‑250 (MQ=255)                                                                                                                                                                 cAGGTCCCGGGGGCACA‑TCCCCCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATcc                                                                              <  2:818558/251‑1 (MQ=255)                                                                                                                                                                      cccGGTGGCACA‑TCCCCCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTGCCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCcagc                                                                          <  1:763082/250‑1 (MQ=255)                                                                                                                                                                             gaaaaacccACCCGACAGAAAAAGTCCCAGCCCGAAGATATACCCACCTAAAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGcac                                                                   <  2:322362/244‑1 (MQ=255)                                                                                                                                                                               aacgcccACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCCCCTAGAACACAACCAAGCCACGGGGACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAGCCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGcacca                                                                 <  2:911379/246‑1 (MQ=255)                                                                                                                                                                                               aaCAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGat                                                                                                                                                    >  2:942222/1‑152 (MQ=255)                                                                                                                                                                                               aaCAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGat                                                                                                                                                    <  1:942222/152‑1 (MQ=255)                                                                                                                                                                                               aaCAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCtt                                                  >  1:1207830/1‑250 (MQ=255)                                                                                                                                                                                                     tCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATc                                            >  1:219889/1‑250 (MQ=255)                                                                                                                                                                                                     ccccagtacaaaaataatcccACCAATAACAGTACCAAGCCACGGTAACGCCCCAACTTCGTAATTCATTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCa                                           <  2:882308/233‑1 (MQ=255)                                                                                                                                                                                                       ccAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCAtt                                         >  1:39660/1‑251 (MQ=255)                                                                                                                                                                                                          gCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGaa                                      <  1:23546/251‑1 (MQ=255)                                                                                                                                                                                                              gAAGATATCCCCAACGATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGACAGTAGGCCCGCCTGAATTAACGCAAAAACCACCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATGTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACa                                  <  2:1108456/251‑1 (MQ=255)                                                                                                                                                                                                              gAAGATATACCCACCTATAACAGTACCAATCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTGCGCGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGGGCACCACAAATTAGCCCGCTTAATATCATTGAAAACa                                  <  1:689523/251‑1 (MQ=255)                                                                                                                                                                                                               aaGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACAt                                 >  1:478233/1‑251 (MQ=255)                                                                                                                                                                                                                              ataACAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCTCAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAAcc                    <  1:92756/249‑1 (MQ=255)                                                                                                                                                                                                                                aaCAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAAcc                    >  2:176321/1‑247 (MQ=255)                                                                                                                                                                                                                                aaCAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAAcc                    >  1:506169/1‑247 (MQ=255)                                                                                                                                                                                                                                  cAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACaa                                                              <  2:977095/203‑1 (MQ=255)                                                                                                                                                                                                                                  cAGTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACaa                                                              >  1:977095/1‑203 (MQ=255)                                                                                                                                                                                                                                    gtcccAAGCCACGTGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAAt            <  1:571874/248‑1 (MQ=255)                                                                                                                                                                                                                                    gTACCAAGCCACGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAAt            >  2:30039/1‑251 (MQ=255)                                                                                                                                                                                                                                              aCGGGAACGCCCCAGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCACACGCTTTGTACAGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACAAATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAATGTTAATGAgg  >  1:609396/1‑251 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                              ACAAATTAAAAGTCTCAGCAATAGAGTTATGTTCAGTACTGTAGTGCTGCAAGGTTTGATTTAAA‑CCACTGGCATGTGGAGAACGCATCACCGCACTCATCACCATATAAGTGAAAA‑GCGCGATCCAACTGCCGATCAATCCCTCACCCGCGCGATACCAGGTCCCGGTGGCACA‑TCCACCCGCCAGAACAATTCCCAGCCCGAAGATATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAATGTTAATGAGG  >  NC_000913/2085072‑2085556

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