Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,569,172	G→A	pseudogene (175/921 nt)	yjiP →	pseudogene, transposase_31 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,569,172	0	G	A	100.0%	307.1 / NA	78	pseudogene (175/921 nt)	yjiP	pseudogene, transposase_31 family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (36/42);  total (36/42)

GCATGATGCATTATTTAAAACCTTTCTCACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGGCGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAGCGCCATA  >  NC_000913/4569024‑4569316                                                                                                                                                     |                                                                                                                                                 gcATGATGCATTATTTAAAACCTTTCTCACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGacg                                                                                                                                                >  1:842458/1‑151 (MQ=255)   atgatgCATTATTTAAAACCTTTCTCACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCt                                                                                                                                              >  1:444241/1‑151 (MQ=255)     gatgCATTATTTAAAACCTTTCTCACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTa                                                                                                                                             <  1:796896/150‑1 (MQ=255)      atgCATTATTTAAAACCTTTCTCACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTaca                                                                                                                                           <  1:92712/151‑1 (MQ=255)       tgCATTATTTAAAACCTTTCTCACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTacac                                                                                                                                          >  1:631643/1‑151 (MQ=255)        gCATTATTTAAAACCTTTCTCACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTacac                                                                                                                                          <  2:1124734/150‑1 (MQ=255)        gCATTATTTAAAACCTTTCTCACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACt                                                                                                                                         <  2:791508/151‑1 (MQ=255)              tttAAAACCTTTCTCACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGTTTTACGTGAACTGTGCGATCTCGATAGCTTAAACCTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGat                                                                                                                                    >  1:1113379/1‑150 (MQ=255)                 aaaaCCTTTCTCACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATAttt                                                                                                                                <  2:680649/151‑1 (MQ=255)                 aaaaCCTTTCTCACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATAtt                                                                                                                                 <  2:117432/150‑1 (MQ=255)                  aaaCCTTTCTCACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATAtttt                                                                                                                               <  2:806295/151‑1 (MQ=255)                  aaaCCTTTCTCACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATAtttt                                                                                                                               <  2:381084/151‑1 (MQ=255)                   aaCCTTTCTCACGCACCCTGACACTGCGCGGGATTTTATGGATATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTa                                                                                                                              <  2:684442/151‑1 (MQ=255)                   aaCCTTTCTCACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTa                                                                                                                              <  1:608989/151‑1 (MQ=255)                     ccTTTCTCACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATg                                                                                                                            <  1:597253/151‑1 (MQ=255)                     ccTTTCTCACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATg                                                                                                                            >  2:1109841/1‑151 (MQ=255)                        ttCTCACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTc                                                                                                                         >  2:220473/1‑151 (MQ=255)                         tctcACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCg                                                                                                                        >  2:858626/1‑151 (MQ=255)                         tctcACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCg                                                                                                                        >  1:35712/1‑151 (MQ=255)                          ctcACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCgg                                                                                                                       >  2:442417/1‑151 (MQ=255)                           tcACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGt                                                                                                                      >  2:239323/1‑151 (MQ=255)                            cACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCgg                                                                                                                       <  2:632238/149‑1 (MQ=255)                            cACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGt                                                                                                                      >  1:323444/1‑150 (MQ=255)                            cACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGt                                                                                                                      >  1:94509/1‑150 (MQ=255)                             aCGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTa                                                                                                                     >  2:44195/1‑150 (MQ=255)                             aCGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTAAGTGAAATGGGCGATATAGATAGATGAAAAAGGGAAGAAGCAAGATTTGGCGATGAAAAAGTGCGGACGATAAACGCAGAGATTTGATGGTAGGGaa                                                                                                                    <  2:85495/151‑1 (MQ=255)                                 aCCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGAc                                                                                                                >  1:90019/1‑151 (MQ=255)                                  cccTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGAcc                                                                                                               >  2:619593/1‑151 (MQ=255)                                       acacTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGa                                                                                                          <  2:323385/151‑1 (MQ=255)                                         acTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAg                                                                                                        <  2:482224/151‑1 (MQ=255)                                          cTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAgg                                                                                                       >  2:318108/1‑151 (MQ=255)                                            gcgcGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGc                                                                                                      <  1:543947/150‑1 (MQ=255)                                            gcgcGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCg                                                                                                     <  2:451988/151‑1 (MQ=255)                                            gcgcGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCg                                                                                                     <  2:520540/151‑1 (MQ=255)                                              gcgGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGAt                                                                                                   <  1:561827/151‑1 (MQ=255)                                                 ggATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGc                                                                                                <  1:416649/151‑1 (MQ=255)                                                        aTGGAGCTTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCt                                                                                         >  2:770653/1‑151 (MQ=255)                                                             gaTTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTaaa                                                                                                                   >  1:7959/1‑120 (MQ=255)                                                             gaTTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTaaa                                                                                                                   <  2:7958/120‑1 (MQ=255)                                                             gaTTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTAtgtg                                                                                    <  2:574203/151‑1 (MQ=255)                                                                     tgccTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACGCTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTg                                                                             <  2:196575/148‑1 (MQ=255)                                                                       ccTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAAc                                                                          >  1:666901/1‑151 (MQ=255)                                                                         tAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACAt                                                                        >  1:108483/1‑151 (MQ=255)                                                                          aaaGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATc                                                                       <  1:468510/151‑1 (MQ=255)                                                                            aGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCag                                                                     <  2:709744/151‑1 (MQ=255)                                                                             gATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATc                                                                       >  1:394900/1‑148 (MQ=255)                                                                             gATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCa                                                                      <  1:267565/149‑1 (MQ=255)                                                                               tttACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGc                                                                  >  1:1119676/1‑151 (MQ=255)                                                                                    gTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCcgcg                                                              >  2:131822/1‑150 (MQ=255)                                                                                    gTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCcgcg                                                              >  2:931717/1‑150 (MQ=255)                                                                                        aCTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGAc                                                         >  2:260605/1‑151 (MQ=255)                                                                                               gATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCat                                                   >  2:917620/1‑150 (MQ=255)                                                                                                aTCTCGATAGCTTAAACCTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACGCTCCGATATTTTCTGGTCGGTAAAGCCCCGCGACGGCGATGGCCATATCTATGTGGTGATTGAACATCATCGCCGCCAGGACATTCa                                                    <  2:281411/148‑1 (MQ=255)                                                                                                  ctcGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATgg                                               >  1:478631/1‑151 (MQ=255)                                                                                                    cGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCata                                                  <  2:444141/146‑1 (MQ=255)                                                                                                    cGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATgg                                               <  1:901684/149‑1 (MQ=255)                                                                                                    cGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATg                                                <  1:252818/148‑1 (MQ=255)                                                                                                    cGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATg                                                <  2:156761/148‑1 (MQ=255)                                                                                                    cGATAGCTTAAAACTGGAATCCGCCAGCTGTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGc                                              <  2:639603/150‑1 (MQ=255)                                                                                                      aTAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCtt                                           <  2:282064/151‑1 (MQ=255)                                                                                                        aGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTc                                         >  1:672738/1‑151 (MQ=255)                                                                                                            tAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGc                                       <  1:1013782/149‑1 (MQ=255)                                                                                                            tAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCt                                     >  1:996446/1‑151 (MQ=255)                                                                                                              aaaCTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCTTTTCGTCt                                     >  2:709535/1‑149 (MQ=255)                                                                                                              aaaCTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGa                                   <  1:173220/151‑1 (MQ=255)                                                                                                                     aaTCCGCCAGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCt                             >  2:416571/1‑150 (MQ=255)                                                                                                                             aGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATgg                    <  2:1024323/151‑1 (MQ=255)                                                                                                                             aGCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATgg                    <  2:713671/151‑1 (MQ=255)                                                                                                                              gCTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGc                   <  2:555688/151‑1 (MQ=255)                                                                                                                               cTTTGTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCg                  >  1:777358/1‑151 (MQ=255)                                                                                                                                   gTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTg               <  2:595403/150‑1 (MQ=255)                                                                                                                                   gTCGATGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGa              <  2:776879/151‑1 (MQ=255)                                                                                                                                       aTGAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGt           >  1:919194/1‑150 (MQ=255)                                                                                                                                         gAAAAATTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAgc        <  1:185696/151‑1 (MQ=255)                                                                                                                                            aaaTTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAGCGcc     >  1:391516/1‑151 (MQ=255)                                                                                                                                             aaTTGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAGCGcc     >  1:287441/1‑150 (MQ=255)                                                                                                                                               ttGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAGCGCCat   >  2:984961/1‑150 (MQ=255)                                                                                                                                               ttGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAACCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAGCGCCat   >  1:112790/1‑150 (MQ=255)                                                                                                                                                tGCGGACGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAGCGCCata  <  2:1048251/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                 GCATGATGCATTATTTAAAACCTTTCTCACGCACCCTGACACTGCGCGGGATTTTATGGAGATTCACTTACCTAAAGATTTACGTGAACTGTGCGATCTCGATAGCTTAAAACTGGAATCCGCCAGCTTTGTCGATGAAAAATTGCGGGCGCTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAGCGCCATA  >  NC_000913/4569024‑4569316

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