Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,475,613	T→C	S36G (AGC→GGC)	tusD ←	sulfurtransferase complex subunit TusD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,475,613	0	T	C	100.0%	224.9 / NA	62	S36G (AGC→GGC)	tusD	sulfurtransferase complex subunit TusD
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (37/25);  total (37/25)

CCGCTACGCAGATATTCAGCGCCACACCATGTTGCGCATTCAGTTGTTGCCAGGCCCGTACGAGGTCAAATTCGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCTGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTGTTGCGTACCGTATGCTGGCCCGGTCACCACGATGGCAAAACGCATTACTTATCTTGCCCCTGGAAATCACCGCTCT  >  NC_000913/3475478‑3475749                                                                                                                                        |                                                                                                                                         ccGCTACGCAGATATTCAGCGCCACACCATGTTGCGCATTCAGTTGTTGCCAGGCCCGTACGAGGTCAAATTCGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCCGCtt                                                                                                                                      >  1:47821/1‑140 (MQ=255)    cTACGCAGATATTCAGCGCCACACCATGTTGCGCATTCAGTTGTTGCCAGGCCCGTACGAGGTCAAATTCGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCCGCTTAAc                                                                                                                                   >  1:257852/1‑140 (MQ=255)              atTCAGCGCCACACCATGTTGCGCATTCAGTTGTTGCCAGGCCCGTACGAGGTCAAATTCGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACgccgct                                                                                                                                       >  1:135329/1‑126 (MQ=255)              atTCAGCGCCACACCATGTTGCGCATTCAGTTGTTGCCAGGCCCGTACGAGGTCAAATTCGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACgccgct                                                                                                                                       <  2:135329/126‑1 (MQ=255)               ttCAGCGCCACACCATGTTGCGCATTCAGTTGTTGCCAGGCCCGTACGAGGTCAAATTCGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATc                                                                                                                        <  2:938004/140‑1 (MQ=255)                    cgcCACACCATGTTGCGCATTCAGTTGTTGCCAGGCCCGTACGAGGTCAAATTCGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTa                                                                                                                   >  2:752115/1‑140 (MQ=255)                       cacaCCATGTTGCGCATTCAGTTGTTGCCAGGCCCGTACGAGGTCAAATTCGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCa                                                                                                                >  1:882206/1‑140 (MQ=255)                                 tGCGCATTCAGTTGTTGCCAGGCCCGTACGAGGTCAAATTCGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTgcgc                                                                                                      >  1:594657/1‑140 (MQ=255)                                   cgcATTCAGTTGTTGCCAGGCCCGTACGAGGTCAAATTCGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCaa                                                                                                    >  1:842651/1‑140 (MQ=255)                                    gcATTCAGTTGTTGCCAGGCCCGTACGAGGTCAAATTCGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCaaa                                                                                                   <  2:728132/140‑1 (MQ=255)                                    gcATTCAGTTGTTGCCAGGCCCGTACGAGGTCAAATTCGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCaaa                                                                                                   >  1:42042/1‑140 (MQ=255)                                       ttCAGTTGTTGCCAGGCCCGTACGAGGTCAAATTCGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTg                                                                                                <  2:842651/140‑1 (MQ=255)                                            ttgttgCCAGGCCCGTACGAGGTCAAATTCGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAg                                                                                           <  1:536972/140‑1 (MQ=255)                                                      gCCCGTACGAGGTCAAATTCGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGTTAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGc                                                                                 >  2:886726/1‑140 (MQ=255)                                                      gCCCGTACGAGGTCAAATTCGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGc                                                                                 >  2:498863/1‑140 (MQ=255)                                                      gCCCGTACGAGGTCAAATTCGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGc                                                                                 >  1:540623/1‑140 (MQ=255)                                                      gCCCGTACGAGGTCAAATTCGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGc                                                                                 >  2:161541/1‑140 (MQ=255)                                                      gCCCGTACGAGGTCAAATTCGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGc                                                                                 >  2:191448/1‑140 (MQ=255)                                                        ccGTACGAGGTCAAATTCGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTTAAAAGCACTACTCGCCt                                                                               <  1:306602/140‑1 (MQ=255)                                                         cGTACGAGGTCAAATTCGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTg                                                                              <  1:191448/140‑1 (MQ=255)                                                         cGTACGAGGTCAAATTCGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTg                                                                              <  1:161541/140‑1 (MQ=255)                                                          gTACGAGGTCAAATTCGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTGt                                                                             >  2:822636/1‑140 (MQ=255)                                                           tACGAGGTCAAATTCGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTGtt                                                                            >  1:652739/1‑140 (MQ=255)                                                           tACGAGGTCAAATTCGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTGtt                                                                            <  1:581458/140‑1 (MQ=255)                                                           tACGAGGTCAAATTCGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTGtt                                                                            >  2:581458/1‑140 (MQ=255)                                                           tACGAGGTCAAATTCGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTGtt                                                                            >  2:22852/1‑140 (MQ=255)                                                              gAGGTCAAATTCGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTGTTGCg                                                                         <  2:494432/140‑1 (MQ=255)                                                                ggTCAAATTCGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTg                                                                              >  2:3650/1‑133 (MQ=255)                                                                ggTCAAATTCGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTg                                                                              <  1:3650/133‑1 (MQ=255)                                                                ggTCAAATTCGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTGTTGCGTa                                                                       >  1:387222/1‑140 (MQ=255)                                                                 gTCAAATTCGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTGTTGCGTAc                                                                      >  2:863392/1‑140 (MQ=255)                                                                 gTCAAATTCGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTGTTGCGTAc                                                                      >  2:296460/1‑140 (MQ=255)                                                                  tCAAATTCGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTGTTGCGTAcc                                                                     <  2:257852/140‑1 (MQ=255)                                                                       ttCGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTGTTGCGTACCGTATg                                                                <  1:863392/140‑1 (MQ=255)                                                                          gTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTGTTGCGTACCGTATGCTg                                                             >  1:884213/1‑140 (MQ=255)                                                                            cACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTGTTGCGTACCGTATGCTGGc                                                           <  2:618415/140‑1 (MQ=255)                                                                                  ccGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTGTTGCGTACCGTATGCTGGCCCGGTc                                                     <  2:47821/140‑1 (MQ=255)                                                                                       gagGTCAATTGGTTAGCGTTATAGACCCCTTCCCGTTAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTGTTGCGTACCGTATGCTGGCGCGGTcaccac                                                <  1:886726/140‑1 (MQ=255)                                                                                         ggTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTGTTGCGTACCGTATGCTGGCCCGGTCACCACGa                                              >  2:527454/1‑140 (MQ=255)                                                                                          gTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTGTTGCGTACCGTATGCTGGCCCGGTCACCACGAt                                             >  1:136767/1‑140 (MQ=255)                                                                                            cAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCTATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTGTTGCGTACCGTATGCTGGCCCGGTCACCACGATgg                                           >  2:672870/1‑140 (MQ=255)                                                                                              aTTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTGTTGCGTACCGTATGCTGGCCCGGTCACCACGATGGCa                                         <  1:672870/140‑1 (MQ=255)                                                                                               ttGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTGTTGCGTACCGTATGCTGGCCCGGTCACCACGATGGCaa                                        >  2:930163/1‑140 (MQ=255)                                                                                                   ttAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACACGCCTGTTGCGTACCGTATGCTGGCCCGGTCACCACGATGGCAAAACg                                    <  2:520114/140‑1 (MQ=255)                                                                                                     aGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTGTTGCGTACCGTATGCTGGCCCGGTCACCACGATGGCAAAACGCa                                  <  2:42042/140‑1 (MQ=255)                                                                                                     aGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTGTTGCGTACCGTATGCTGGCCCGGTCACCACGATGGCAAAACGCa                                  <  2:594657/140‑1 (MQ=255)                                                                                                         ttATAGACCCCTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTGTTGCGTACCGTATGCTGGCCCGGTCACCACGATGGCAAAACGCATTAc                              <  1:822636/140‑1 (MQ=255)                                                                                                          tataGACCCCTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTGTTGCGTACCGTATGCTGGCCCGGTCACCACGATGGCAAAACGCATTACt                             <  2:387222/140‑1 (MQ=255)                                                                                                                 cccTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTGTTGCGTACCGTATGCTGGCCCGGTCACCACGATGGCAAAACGCATTACTTATCTTg                      >  1:131708/1‑140 (MQ=255)                                                                                                                 cccTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTGTTGCGTACCGTATGCTGGCCCGGTCACCACGATGGCAAAACGCATTACTTATCTTg                      >  2:895612/1‑140 (MQ=255)                                                                                                                 cccTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTGTTGCGTACCGTATGCTGGCCCGGTCACCACGATGGCAAAACGCATTACTTATCTTg                      >  1:882955/1‑140 (MQ=255)                                                                                                                   cTTCCCGATAGAGAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTGTTGCGTACCGTATGCTGGCCCGGTCACCACGATGGCAAAACGCATTACTTATCTTGcc                    >  2:3281/1‑140 (MQ=255)                                                                                                                   cTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTGTTGCGTACCGTATGCTGGCCCGGTCACCACGATGGCAAAACGCATTACTTATCTTGcc                    >  2:218708/1‑140 (MQ=255)                                                                                                                   cTTCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTGTTGCGTACCGTATGCTGGCCCGGTCACCACGATGGCAAAACGCATTACTTATCTTGcc                    >  2:961235/1‑140 (MQ=255)                                                                                                                    ttCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTGTTGCGTACCGTATGCTGGCCCGGTCACCACGATGGCAAAACGCa                                  <  2:766598/125‑1 (MQ=255)                                                                                                                    ttCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTGTTGCGTACCGTATGCTGGCCCGGTCACCACGATGGCAAAACGCa                                  >  1:766598/1‑125 (MQ=255)                                                                                                                    ttCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTGTTGCGTACCGTATGCTGGCCCGGTCACCACGATGGCAAAACGCATTa                               <  1:869475/128‑1 (MQ=255)                                                                                                                    ttCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTGTTGCGTACCGTATGCTGGCCCGGTCACCACGATGGCAAAACGCATTACTTATCt                        <  2:413617/135‑1 (MQ=255)                                                                                                                    ttCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTGTTGCGTACCGTATGCTGGCCCGGTCACCACGATGGCAAAACGCATTACTTATCt                        >  1:413617/1‑135 (MQ=255)                                                                                                                    ttCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTGTTGCGTACCGTATGCTGGCCCGGTCACCACGATGGCAAAACGCATTACTTATCTTGccc                   >  2:755197/1‑140 (MQ=255)                                                                                                                    ttCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTGTTGCGTACCGTATGCTGGCCCGGTCACCACGATGGCAAAACGCATTACTTATCTTGccc                   <  1:895612/140‑1 (MQ=255)                                                                                                                    ttCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTGTTGCGTACCGTATGCTGGCCCGGTCACCACGATGGCAAAACGCATTACTTATCTTGccc                   >  1:121544/1‑140 (MQ=255)                                                                                                                    ttCCCGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTGTTGCGTACCGTATGCTGGCCCGGTCACCACGATGGCAAAACGAATTa                               >  2:869475/1‑128 (MQ=255)                                                                                                                      cccGATAGAAAAAGACGCCGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTGTTGCGTACCGTATGCTGGCCCGGTCACCACGATGGCAAAACGCATTACTTATCTTGCCCCt                 >  1:963393/1‑140 (MQ=255)                                                                                                                                     cgccgcTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTGTTGCGTACCGTATGCTGGCCCGGTCACCACGATGGCAAAACGCATTACTTATCTTGCCCCTGGAAATCACCGctct  >  2:893426/1‑140 (MQ=255)                                                                                                                                     cgccgcTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTGTTGCGTACCGTATGCTGGCCCGGTCACCACGATGGCAAAACGCATTACTTATCTTGCCCCTGGAAATCACCGctct  >  1:500329/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                         CCGCTACGCAGATATTCAGCGCCACACCATGTTGCGCATTCAGTTGTTGCCAGGCCCGTACGAGGTCAAATTCGTCACTTGCCGGAGAGGTCAATTGGTTAGCGTTATAGACCCCTTCCCGATAGAAAAAGACGCTGCTTAACTCATGGCCATCTGCTATCAGCGCCTGCGCAAACTGAAAAGCACTACTCGCCTGTTGCGTACCGTATGCTGGCCCGGTCACCACGATGGCAAAACGCATTACTTATCTTGCCCCTGGAAATCACCGCTCT  >  NC_000913/3475478‑3475749

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑