Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,729,436	G→A	R117Q (CGG→CAG)	lhr →	putative ATP‑dependent helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,729,436	0	G	A	100.0%	315.9 / NA	77	R117Q (CGG→CAG)	lhr	putative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (43/34);  total (43/34)

TCCGCGAAGGCGGCGAAGATACCCGCGAGGCGCATAAGCGTAAAACCTCACGCATCCTCTATATTTCACCGATAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCACTCTATCTGATGCTGACCTCCCG  >  NC_000913/1729295‑1729580                                                                                                                                              |                                                                                                                                                 tCCGCGAAGGCGGCGAAGATACCCGCGAGGCGCATAAGCGTAAAACCTCACGCATCCTCTATATTTCACCGATAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGgcagcggcg                                                                                                                                            <  1:81119/148‑1 (MQ=255)  ccGCGAAGGCGGCGAAGATACCCGCGAGGCGCATAAGCGTAAAACCTCACGCATCCTCTATATTTCACCGATAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAgcggcgcggc                                                                                                                                        >  2:706362/1‑151 (MQ=255)  ccGCGAAGGCGGCGAAGATACCCGCGAGGCGCATAAGCGTAAAACCTCACGCATCCTCTATATTTCACCGATAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAgcggcgcggc                                                                                                                                        <  2:496137/151‑1 (MQ=255)   cgcgAAGGCGGCGAAGATACCCGCGAGGCGCATAAGCGTAAAACCTCACGCATCCTCTATATTTCACCGATAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGcggcgcggcg                                                                                                                                       <  2:675581/151‑1 (MQ=255)       aaGGCGGCGAAGATACCCGCGAGGCGCATAAGCGTAAAACCTCACGCATCCTCTATATTTCACCGATAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAAc                                                                                                                                   <  2:15002/151‑1 (MQ=255)          gcggcgAAGATACCCGCGAGGCGCATAAGCGTAAAACCTCACGCATCCTCTATATTTCACCGATAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGa                                                                                                                                <  1:166866/151‑1 (MQ=255)               gAAGATACCCGCGAGGCGCATAAGCGTAAAACCTCACGCATCCTCTATATTTCACCGATAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTc                                                                                                                            >  1:596334/1‑150 (MQ=255)                    tACCCGCGAGGCGCATAAGCGTAAAACCTCACGCATCCTCTATATTTCACCGATAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCtt                                                                                                                      >  2:118239/1‑151 (MQ=255)                         gcgAGGCGCATAAGCGTAAAACCTCACGCATCCTCTATATTTCACCGATAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGt                                                                                                                 >  1:3961/1‑151 (MQ=255)                         gcgAGGCGCATAAGCGTAAAACCTCACGCATCCTCTATATTTCACCGATAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGt                                                                                                                 <  1:37331/151‑1 (MQ=255)                           gAGGCGCATAAGCGTAAAGCCGCACGCATCCTCTATATTACACCGATAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTAAATCTTCGCGTa                                                                                                                <  2:216916/150‑1 (MQ=255)                               cgcATAAGCGTAAAACCTCACGCATCCTCTATATTTCACCGATAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGAt                                                                                                           >  2:728150/1‑151 (MQ=255)                                gcATAAGCGTAAAACCTCACGCATCCTCTATATTTCACCGATAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATc                                                                                                          >  2:16272/1‑151 (MQ=255)                                  aTAAGCGTAAAACCTCACGCATCCTCTATATTTCACCGATAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCg                                                                                                        >  1:227822/1‑151 (MQ=255)                                  aTAAGCGTAAAACCTCACGCATCCTCTATATTTCACCGATAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCg                                                                                                        >  1:143905/1‑151 (MQ=255)                                   tAAGCGTAAAACCTCACGCATCCTCTATATTTCACCGATAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCg                                                                                                        >  2:285466/1‑150 (MQ=255)                                     aGCGTAAAACCTCACGCATCCTCTATATTTCACCGATAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTAc                                                                                                     <  2:193286/151‑1 (MQ=255)                                      gCGTAAAACCTCACGCATCCTCTATATTTCACCGATAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACt                                                                                                    >  1:219376/1‑151 (MQ=255)                                        gTAAAACCTCACGCATCCTCTATATTTCACCGATAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTCGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTgg                                                                                                  >  1:40843/1‑151 (MQ=255)                                         tAAAACCTCACGCATCCTCTATATTTCACCGATAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGc                                                                                                 >  1:692097/1‑151 (MQ=255)                                          aaaaCCTCACGCATCCTCTATATTTCACCGATAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTgga                                                                                                 <  2:572524/150‑2 (MQ=255)                                           aaaCCTCACGCATCCTCTATATTTCACCGATAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCg                                                                                                >  2:554592/1‑150 (MQ=255)                                           aaaCCTCACGCATCCTCTATATTTCACCGATAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGa                                                                                               <  2:244336/151‑1 (MQ=255)                                            aaCCTCACGCATCCTCTATATTTCACCGATAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGAt                                                                                              <  1:660006/151‑1 (MQ=255)                                             aCCTCACGCATCCTCTATATTTCACCGATAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATa                                                                                             <  1:288961/151‑1 (MQ=255)                                             aCCTCACGCATCCTCTATATTTCACCGATAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATa                                                                                             <  1:274921/151‑1 (MQ=255)                                                tCACGCATCCTCTATATTTCACCGATAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGc                                                                                          >  2:271308/1‑151 (MQ=255)                                                  aCGCATCCTCTATATTTCACCGATAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCt                                                                                        <  2:364651/151‑1 (MQ=255)                                                   cGCATCCTCTATATTTCACCGATAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTg                                                                                       <  1:326021/151‑1 (MQ=255)                                                   cGCATCCTCTATATTTCACCGATAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTg                                                                                       <  2:66218/151‑1 (MQ=255)                                                     cATCCTCTATATTTCACCGATAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGca                                                                                     <  2:305234/151‑1 (MQ=255)                                                       tCCTCTATATTTCACCGATAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGcac                                                                                    >  2:617342/1‑150 (MQ=255)                                                        ccTCTATATTTCACCGATAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGcac                                                                                    >  1:176900/1‑149 (MQ=255)                                                             atatTTCACCGATAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAAc                                                                              <  1:364752/150‑1 (MQ=255)                                                               atTTCACCGATAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACgca                                                                           <  1:766303/151‑1 (MQ=255)                                                               atTTCACCGATAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACgc                                                                            <  1:42201/150‑1 (MQ=255)                                                                tttCACCGATAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACgc                                                                            >  1:447516/1‑149 (MQ=255)                                                                 ttCACCGATAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACgcagc                                                                         >  1:660110/1‑151 (MQ=255)                                                                 ttCACCGATAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACgcagc                                                                         >  2:281653/1‑151 (MQ=255)                                                                  tCACCGATAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACgcagc                                                                         >  1:370575/1‑150 (MQ=255)                                                                       gATAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTc                                                                   <  1:268019/151‑1 (MQ=255)                                                                        aTAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTc                                                                   >  1:225350/1‑150 (MQ=255)                                                                        aTAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTc                                                                   >  1:603185/1‑150 (MQ=255)                                                                         tAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCAc                                                                 >  1:115968/1‑151 (MQ=255)                                                                           aaaGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCAccc                                                               >  2:701421/1‑151 (MQ=255)                                                                            aaGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCg                                                              >  1:355130/1‑151 (MQ=255)                                                                              gCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGt                                                             >  2:308106/1‑150 (MQ=255)                                                                                    ggCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATccgc                                                      <  1:663509/151‑1 (MQ=255)                                                                                     gCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATccgcc                                                     >  1:734691/1‑151 (MQ=255)                                                                                        ccGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCgg                                                   >  2:244221/1‑150 (MQ=255)                                                                                         cGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGat                                                 <  1:313161/151‑1 (MQ=255)                                                                                          gACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGat                                                 >  2:741446/1‑150 (MQ=255)                                                                                          gACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGat                                                 >  1:87842/1‑150 (MQ=255)                                                                                            cGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATAtt                                              <  1:653758/151‑1 (MQ=255)                                                                                               tCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTg                                           >  2:109968/1‑151 (MQ=255)                                                                                                cAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTg                                           <  2:221867/150‑1 (MQ=255)                                                                                                cAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGa                                          <  2:370890/151‑1 (MQ=255)                                                                                                  gcgcAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCACGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATc                                        <  2:39476/151‑1 (MQ=255)                                                                                                   cgcAACTTGCAGCTCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGCCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTg                                           <  2:121136/147‑1 (MQ=255)                                                                                                      aaCTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATcacca                                    <  2:762758/151‑1 (MQ=255)                                                                                                       aCTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATcaccac                                   >  1:414999/1‑151 (MQ=255)                                                                                                         ttGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCacac                                 >  2:216236/1‑151 (MQ=255)                                                                                                         ttGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCacac                                 >  1:387586/1‑151 (MQ=255)                                                                                                          tGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACcaca                                  <  1:627035/149‑1 (MQ=255)                                                                                                            cAGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCTGCCGGATATTCTGATCACCACACCCg                              >  2:524858/1‑151 (MQ=255)                                                                                                             aGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATcc                                                        >  1:720240/1‑124 (MQ=255)                                                                                                             aGATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATcc                                                        <  2:720196/124‑1 (MQ=255)                                                                                                              gATCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGaa                            <  2:726072/151‑1 (MQ=255)                                                                                                               aTCCAGTTCAAGGGTATTGCCGATGAACGGCAGAGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAAt                           <  2:322433/151‑1 (MQ=255)                                                                                                                tCCCGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATc                          >  1:212005/1‑151 (MQ=255)                                                                                                                 cccGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATc                          >  1:248969/1‑150 (MQ=255)                                                                                                                   cGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCAct                       >  1:730016/1‑151 (MQ=255)                                                                                                                   cGTTGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCAct                       <  1:414732/151‑1 (MQ=255)                                                                                                                      tGAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCActct                     <  2:716682/150‑1 (MQ=255)                                                                                                                       gAAGGGTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCACTCTAt                   >  1:375790/1‑151 (MQ=255)                                                                                                                           ggTATTGCCGATGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCACTCTATCTg                >  1:497350/1‑150 (MQ=255)                                                                                                                                      tGAACGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCACTCTATCTGATGCTGACCTcc    >  1:141364/1‑151 (MQ=255)                                                                                                                                        aaCGGCAGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCACTCTATCTGATGCTGACCTCCcg  >  2:746021/1‑151 (MQ=255)                                                                                                                                              |                                                                                                                                                 TCCGCGAAGGCGGCGAAGATACCCGCGAGGCGCATAAGCGTAAAACCTCACGCATCCTCTATATTTCACCGATAAAAGCCCTGGGCACCGACGTTCAGCGCAACTTGCAGATCCCGTTGAAGGGTATTGCCGATGAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCACTCTATCTGATGCTGACCTCCCG  >  NC_000913/1729295‑1729580

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