Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,196,057	C→T	A15T (GCC→ACC)	thiC ←	phosphomethylpyrimidine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,196,057	0	C	T	100.0%	263.8 / NA	77	A15T (GCC→ACC)	thiC	phosphomethylpyrimidine synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (36/41);  total (36/41)

CCGCCAATTAGCGTCGGGCTAAGCTGGATCTCACGCATCGGCACGCGCACGCCGGGGTGTGTGCCAGTGATATAAATGCGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATT  >  NC_000913/4195927‑4196187                                                                                                                                   |                                                                                                                                    ccGCCAATTAGCGTCGGGCTAAGCTGGATCTCACGCATCGGCACGCGCACGCCGGGGTGTGTGCCAGTGATATAAATGCGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGcgcg                                                                                                                            >  1:644638/1‑140 (MQ=255) ccGCCAATTAGCGTCGGGCTAAGCTGGATCTCACGCATCGGCACGCGCACGCCGGGGTGTGTGCCAGTGATATAAATGCGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGcgcg                                                                                                                            <  2:1414607/140‑1 (MQ=255)    ccAATTAGCGTCGGGCTAAGCTGGATCTCACGCATCGGCACGCGCACGCCGGGGTGTGTGCCAGTGATATAAATGCGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGGGTCGATAAAATGTTGGGTCCGGGCGCgttg                                                                                                                         >  2:1029094/1‑140 (MQ=255)      aaTTAGCGTCGGGCTAAGCTGGATCTCACGCATCGGCACGCGCACGCCGGGGTGTGTGCCAGTGATATAAATGCGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCgttgtt                                                                                                                       >  1:1158499/1‑140 (MQ=255)            cGTCGGGCTAAGCTGGATCTCACGCATCGGCACGCGCACGCCGGGGTGTGTGCCAGTGATATAAATGCGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCgcggc                                                                                                                 >  1:967412/1‑140 (MQ=255)              tCGGGCTAAGCTGGATCTCACGCATCGGCACGCGCACGCCGGGGTGTGTGCCAGTGATATAAATGCGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGcggcgg                                                                                                               >  1:1594618/1‑140 (MQ=255)                 ggCTAAGCTGGATCTCACGCATCGACACGCGCACGCCGGGGTGTGTGCCAGTGATATAAATGCGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCggggc                                                                                                            >  2:128291/1‑138 (MQ=255)                       gCTGGATCTCACGCATCGGCACGCGCACGCCGGGGTGTGTGCCAGTGATATAAATGCGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGTCAGtttt                                                                                                      <  2:664267/140‑1 (MQ=255)                       gCTGGATCTCACGCATCGGCACGCGCACGCCGGGGTGTGTGCCAGTGATATAAATGCGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGTCAGtttt                                                                                                      <  1:1591715/140‑1 (MQ=255)                            aTCTCACGCATCGGCACGCGCACGCCGGGGTGTGTGCCAGTGATATAAATGCGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGc                                                                                                 <  2:682550/140‑1 (MQ=255)                            aTCTCACGCATCGGCACGCGCACGCCGGGGTGTGTGCCAGTGATATAAATGCGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGGCGGGTtttttgg                                                                                                 >  2:348633/1‑139 (MQ=255)                                cACGCATCGGCACGCGCACGCCGGGGTGTGTGCCAGTGATATAAATGCGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGTGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGGCATTTTTTTTGCAGAc                                                                                             >  1:1239251/1‑140 (MQ=255)                                 aCGCATCGGCACGCGCACGCCGGGGTGTGTGCCAGTGATATAAATGCGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGGGAGTGtttttggcgaaa                                                                                            >  2:1065201/1‑134 (MQ=255)                                  cGCATCGGCACGCGCACGCCGGGGTGTGTGCCAGTGATATAAATGCGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGGCGGGTTTTTTGccaacaa                                                                                           >  2:180276/1‑134 (MQ=255)                                       cGGCACGCGCACGCCGGGGTGTGTGCCAGTGATATAAATGCGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGGCAGGTTTTTTTCAGACATAGCTc                                                                                      >  2:1606278/1‑140 (MQ=255)                                        ggCACGCGCACGCCGGGGTGTGTGCCAGTGATATAAATGCGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGGCAGTTTTTTTTCAGACATAGCTcc                                                                                     >  1:154367/1‑139 (MQ=255)                                        ggCACGCGCACGCCGGGGTGTGTGCCAGTGATATAAATGCGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGGCAGGTTTGTTTCAGACAGAGCTcc                                                                                     >  2:1187620/1‑139 (MQ=255)                                         gcacgcgcacgcCGGGGTGTGTGCCAGTGATATAAATGCGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGTTGAGGTTTGTTAAAGACACAGGTCaa                                                                                    >  2:269785/1‑139 (MQ=255)                                             gcgcacgcCGGGGTGTGTGCCAGTGATATAAATGCGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCa                                                                               <  1:994578/140‑1 (MQ=255)                                              cgcacgcCGGGGTGTGTGCCAGTGATATAAATGCGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCaa                                                                              <  1:8793/140‑1 (MQ=255)                                              cgcacgcCGGGGTGTGTGCCAGTGATATAAATGCGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCaa                                                                              <  2:1552439/140‑1 (MQ=255)                                               gcacgcCGGGGTGTGTGCCAGTGATATAAATGCGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCaaa                                                                             <  2:967412/140‑1 (MQ=255)                                               gcacgcCGGGGTGTGTGCCAGTGATATAAATGCGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCaaa                                                                             <  2:1539055/140‑1 (MQ=255)                                               gcacgcCGGGGTGTGTGCCAGTGATATAAATGCGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCaaa                                                                             <  2:211227/140‑1 (MQ=255)                                                cacgcCGGGGTGTGTGCCAGTGATATAAATGCGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGGCAGTTTTTTTTCAGACATAGCTCA‑TTCCaaaa                                                                            >  2:82363/1‑140 (MQ=255)                                                   gcCGGGGTGTGTGCCAGTGATATAAATGCGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGGCCGGTTTTTTGCTGACATAGCCCC‑TTCCCAAAAga                                                                         >  2:354999/1‑139 (MQ=255)                                                    ccGGGGTGTGTGCCAGTGATATAAATGCGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCAAAAAGtt                                                                        <  2:1019090/140‑1 (MQ=255)                                                        ggTGTGTGCCAGTGATATAAATGCGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGGCAGTTTTTTTGCAGACATAGCTCC‑TTCCaaaaaaataaag                                                                    >  2:284610/1‑133 (MQ=255)                                                            tgtgCCAGTGATATAAATGCGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCAAAAAGTTAAGGACGt                                                                <  1:1432058/140‑1 (MQ=255)                                                             gtgCCAGTGATATAAATGCGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCAAAAAGTTAAGGACGTg                                                               <  2:1273672/140‑1 (MQ=255)                                                             gtgCCAGTGATATAAATGCGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGGCAGGTTTTTTGCAGACATAACTCA‑TTTCAAAAAAGTTAGGACggg                                                               >  1:1357247/1‑138 (MQ=255)                                                              tgCCAGTGATATAAATGCGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCAAAAAGTTAAGGACGTgg                                                              <  2:227794/140‑1 (MQ=255)                                                              tgCCAGTGATATAAATGCGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGGCAGGTTTGTTGCAGACAAAGCTCATTTCCAAAAAG‑TAAGGACGTgg                                                              >  2:1304240/1‑140 (MQ=255)                                                                ccAGTGATATAAATGCGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCAAAAAGTTAAGGACGTGGCt                                                            <  1:1103297/140‑1 (MQ=255)                                                                    tGATATAAATGCGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTTTTGCAGACATAGCTCA‑TTCCCAAAAAATAAGGACGTGGGTTGTc                                                        >  1:435891/1‑140 (MQ=255)                                                                     gATATAAATGCGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCAAAAAGTTAAGGACGTGGCTTGTCa                                                       <  2:406720/140‑1 (MQ=255)                                                                        ataAATGCGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGGCCGATTTTTTTCAAACATAGGACA‑TTTCCAAAAATTAAAGAAGGTGCTTTTCAgac                                                    >  1:453383/1‑140 (MQ=255)                                                                            aTGCGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGGCAGGTTTTTTTCAGACAAAGCTCC‑TTCCCAAAAAATAAGGGCGGGGGTTGGCAGacggcg                                                >  2:263839/1‑140 (MQ=255)                                                                              gCGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGa                                              <  1:180276/140‑1 (MQ=255)                                                                              gCGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGGCGGGTTTTTTTCAGACATAGCTCA‑TTCCAAAAAAGTAAGGGCGGGGGTTGTCAGACGACggg                                              >  2:360045/1‑139 (MQ=255)                                                                              gCGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGGCAGGTTTGTTTCAGACATAGTTTA‑TTTCCAAAAATTAAGGACGGGGGTTGTTCGAAGACggg                                              >  1:351509/1‑139 (MQ=255)                                                                               cGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGAt                                             <  2:320639/140‑1 (MQ=255)                                                                               cGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGAt                                             <  2:859698/140‑1 (MQ=255)                                                                               cGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGAt                                             <  1:1501655/140‑1 (MQ=255)                                                                               cGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGAt                                             <  1:98627/140‑1 (MQ=255)                                                                               cGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGGCAGGTTTGTTGCAGACATAGCTCA‑TTCCAAAAAAATAAGGACGGGGGTTGTCAGACGACgggg                                             >  2:596401/1‑138 (MQ=255)                                                                                gTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGGCAGTTTTTTTGCAGACATAGCCCC‑TTTCCAAAAATTAAGGACGTGGGGTGTCAGAcggcggggg                                            >  2:1164723/1‑137 (MQ=255)                                                                                 ttttGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGa                                           >  2:1491050/1‑140 (MQ=255)                                                                                 ttttGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGa                                           <  2:1158499/140‑1 (MQ=255)                                                                                 ttttGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGGCCGTTTTTTTTCAGACATAGCTCC‑TTCCCAAAAAATAAGGACGTGGCGTTTCAGACGACGGaggg                                           >  2:32424/1‑137 (MQ=255)                                                                                 ttttGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGGCAGTTTTGTTGCAGACATAGCTCA‑TTCCAAAAAGTTAAGGAAGCGGCGTGTCAGACGCCGGGTGa                                           >  2:39932/1‑140 (MQ=255)                                                                                  tttGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGGCAGTTTTTTTTCCGACATAGCTCA‑TTCCCAAAAAATAAGGACGGGGGGTGTTAGACGACGGGTGaa                                          >  2:68773/1‑140 (MQ=255)                                                                                  tttGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGGCAGGTTTTTTGCAGACATAGCTCA‑TTCCAAAAAGTTAAAGACGGGGGGTGTCAGACGACgggagaa                                          >  2:100879/1‑135 (MQ=255)                                                                                   ttGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAg                                         <  2:25254/140‑1 (MQ=255)                                                                                         tGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGTCCGTTTTGTTGCAGACATAGCTCA‑TTCCAAAAAAATAAGGACGGGGCTTGACAGACGACGGATGAAGCAagag                                   >  2:36749/1‑140 (MQ=255)                                                                                              aaGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGAt                              <  1:791537/140‑1 (MQ=255)                                                                                              aaGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGAt                              <  2:1471133/140‑1 (MQ=255)                                                                                              aaGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGAt                              <  1:496284/140‑1 (MQ=255)                                                                                                ggCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATcg                            <  1:343260/140‑1 (MQ=255)                                                                                                ggCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGTCAGGTTTGTTTTGAACACAGCTCC‑TTTTCAAAAACTAAAGACGTGGCGTGTTTGACGGCGGGAGGAGGAAGAGAccaacc                            >  2:907049/1‑135 (MQ=255)                                                                                                    gTGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATcgccgc                        <  2:569363/140‑1 (MQ=255)                                                                                                     tGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCa                       <  2:489054/140‑1 (MQ=255)                                                                                                     tGCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGGCAGTTTTTTTTCAGACATAACTCA‑TTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGAGGAAGGAAGAGAAGATcgccgcc                       >  2:28122/1‑139 (MQ=255)                                                                                                      gCCTTCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAg                      <  1:992751/140‑1 (MQ=255)                                                                                                         ttCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAgggg                   <  1:1304240/140‑1 (MQ=255)                                                                                                          tCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGgc                  <  2:668382/140‑1 (MQ=255)                                                                                                          tCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGgc                  <  1:519288/140‑1 (MQ=255)                                                                                                          tCCAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGgc                  <  1:1065201/140‑1 (MQ=255)                                                                                                            cAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGgcga                <  1:1480855/140‑1 (MQ=255)                                                                                                            cAGGGTGTCGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGgcga                <  2:1181763/140‑1 (MQ=255)                                                                                                                    cGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGGCAGTTTTTTTGCAGACATAGCTCA‑TTCCCAAAAGTTAAGGACGGGGCTTGTCAGACGGCGGGTGAAGCAAGAGGCGATCGCCGCCGGGGCGGGGCGGTAg        >  2:732866/1‑140 (MQ=255)                                                                                                                    cGATAAAATGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGGCAGGTTTTTTTTCGAAATAGATCA‑TTACCAAAAAATAAGGGAGTGGGTTGTTAGAAGACGGAAGGATCAAGAGACGCTGGCCGCCCGAGGGGTGCGATaa        >  2:548826/1‑139 (MQ=255)                                                                                                                          aaTGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGAtt  <  2:666443/140‑1 (MQ=255)                                                                                                                          aaTGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGAtt  <  1:166522/140‑1 (MQ=255)                                                                                                                          aaTGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGAtt  <  2:1357247/140‑1 (MQ=255)                                                                                                                          aaTGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGAtt  <  2:1646590/140‑1 (MQ=255)                                                                                                                          aaTGTTGGGTCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGAtt  <  1:711825/140‑1 (MQ=255)                                                                                                                                   |                                                                                                                                    CCGCCAATTAGCGTCGGGCTAAGCTGGATCTCACGCATCGGCACGCGCACGCCGGGGTGTGTGCCAGTGATATAAATGCGTTTTGAGTTGGGAAAGGCGGTGCCTTCCAGGGTGTCGATAAAATGTTGGGCCCGGGCGCGTTGTTCGCGGCGGGTCAGTTTTGTTGCAGACATAGCTCA‑TTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATT  >  NC_000913/4195927‑4196187

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