Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,281,585	T→C	pseudogene (464/504 nt)	agaA →	putative truncated N‑acetylgalactosamine‑6‑phosphate deacetylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,281,585	0	T	C	100.0%	246.9 / NA	65	pseudogene (464/504 nt)	agaA	putative truncated N‑acetylgalactosamine‑6‑phosphate deacetylase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (24/41);  total (24/41)

GGCGTAACGCCTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAA  >  NC_000913/3281452‑3281718                                                                                                                                      |                                                                                                                                      ggCGTAACGCCTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGcaaca                                                                                                                                 <  1:796545/140‑1 (MQ=255)            tGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTc                                                                                                                      >  2:425550/1‑140 (MQ=255)              cGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAg                                                                                                                    <  1:632658/140‑1 (MQ=255)              cGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAg                                                                                                                    <  2:199497/140‑1 (MQ=255)                     catccatATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAAt                                                                                                             <  1:580385/140‑1 (MQ=255)                         catATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGc                                                                                                         <  2:554710/140‑1 (MQ=255)                         catATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGc                                                                                                         <  1:905808/140‑1 (MQ=255)                           tatGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCtt                                                                                                       <  2:1116947/140‑1 (MQ=255)                           tatGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCtt                                                                                                       <  1:228374/140‑1 (MQ=255)                              ggCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGt                                                                                                                 >  2:1122131/1‑127 (MQ=255)                              ggCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGt                                                                                                                 <  1:1122131/127‑1 (MQ=255)                              ggCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGt                                                                                                    >  1:1120002/1‑140 (MQ=255)                                cGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGttt                                                                                                  >  1:556523/1‑140 (MQ=255)                                  tCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTg                                                                                                >  2:833581/1‑140 (MQ=255)                                  tCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTg                                                                                                >  1:958694/1‑140 (MQ=255)                                  tCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTg                                                                                                >  1:691160/1‑140 (MQ=255)                                   cGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGa                                                                                               >  1:150586/1‑140 (MQ=255)                                   cGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGa                                                                                               >  2:559692/1‑140 (MQ=255)                                     cTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATa                                                                                             >  2:703419/1‑140 (MQ=255)                                        cATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGtt                                                                                          <  1:224719/140‑1 (MQ=255)                                            cGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCt                                                                                      <  1:813813/140‑1 (MQ=255)                                               cgcgAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCt                                                                                   <  2:556523/140‑1 (MQ=255)                                                   aaTGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTAt                                                                               <  1:833581/140‑1 (MQ=255)                                                      gCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTgg                                                                            <  1:1047855/140‑1 (MQ=255)                                                       cTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTggg                                                                           <  1:703419/140‑1 (MQ=255)                                                        tGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGc                                                                          <  1:19092/140‑1 (MQ=255)                                                           gtgtTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCtt                                                                       <  2:156664/140‑1 (MQ=255)                                                                 aTGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTc                                                                 <  1:559692/140‑1 (MQ=255)                                                                 aTGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTc                                                                 <  2:1058501/140‑1 (MQ=255)                                                                  tGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTcc                                                                <  2:691160/140‑1 (MQ=255)                                                                  tGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAACTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTcc                                                                <  2:161214/140‑1 (MQ=255)                                                                    gtgtTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTcccc                                                              >  1:245891/1‑140 (MQ=255)                                                                      gtTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCg                                                            <  1:425550/140‑1 (MQ=255)                                                                      gtTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCg                                                            >  2:92797/1‑140 (MQ=255)                                                                               tCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTCCCCCCGTAAGGCctt                                                   >  2:473469/1‑140 (MQ=255)                                                                                 gCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCctttc                                                 >  2:444496/1‑140 (MQ=255)                                                                                 gCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCctttc                                                 >  2:833397/1‑140 (MQ=255)                                                                                         cGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCtttctttttctt                                         <  2:958694/140‑1 (MQ=255)                                                                                         cGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCtttctttttcta                                         <  1:833397/140‑2 (MQ=255)                                                                                           ggCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCtttctttttctt                                         <  2:864973/138‑1 (MQ=255)                                                                                           ggCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCtttctttttctt                                         >  1:864973/1‑138 (MQ=255)                                                                                            gCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCg                                      >  1:911386/1‑140 (MQ=255)                                                                                              aaaCGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGtt                                    <  1:92797/140‑1 (MQ=255)                                                                                              aaaCGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGtt                                    <  1:347446/140‑1 (MQ=255)                                                                                                   cgcCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGAt                               <  2:911386/140‑1 (MQ=255)                                                                                                   cgcCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGAt                               <  2:1120002/140‑1 (MQ=255)                                                                                                   cgcCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGAt                               <  1:317600/140‑1 (MQ=255)                                                                                                   cgcCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGAt                               <  2:245891/140‑1 (MQ=255)                                                                                                   cgcCAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTCCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTAGAt                               <  1:473469/140‑1 (MQ=255)                                                                                                   cgcCAGCGTCGTTGCGCTGGATAGCGGACTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGAt                               <  1:558045/140‑1 (MQ=255)                                                                                                     ccAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCt                             >  1:1048741/1‑140 (MQ=255)                                                                                                     ccAGCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCt                             <  1:444496/140‑1 (MQ=255)                                                                                                        gCGTCGTTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACGTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCtgtg                          >  1:115837/1‑140 (MQ=255)                                                                                                            cgtTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGc                      <  2:1048741/140‑1 (MQ=255)                                                                                                             gtTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCg                     <  1:789254/140‑1 (MQ=255)                                                                                                             gtTGCGCTGGATAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCg                     <  2:150586/140‑1 (MQ=255)                                                                                                                       aTAGCGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATg           <  1:867119/140‑1 (MQ=255)                                                                                                                           cGGGCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGAc       >  1:96000/1‑140 (MQ=255)                                                                                                                             ggCTACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACGTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGgcgc     <  2:115837/140‑1 (MQ=255)                                                                                                                                tACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGgcgc     <  2:309135/137‑1 (MQ=255)                                                                                                                                tACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGgcgc     >  1:309135/1‑137 (MQ=255)                                                                                                                                tACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTaa  >  1:712931/1‑140 (MQ=255)                                                                                                                                tACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTaa  <  1:842115/140‑1 (MQ=255)                                                                                                                                tACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTaa  >  2:422158/1‑140 (MQ=255)                                                                                                                                tACATGCGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTaa  >  2:842115/1‑140 (MQ=255)                                                                                                                                      |                                                                                                                                      GGCGTAACGCCTGCGGAAGCCATCCATATGGCGTCGCTGCATCCGGCGCGAATGCTGGGTGTTGATGGTGTTCTGGGATCGCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAA  >  NC_000913/3281452‑3281718

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