BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,615,859	T→G	100%	T403P (ACC→CCC)	yjjI ←	DUF3029 family protein, putative glycine radical enzyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,615,859	0	T	G	100.0%	614.4 / ‑6.7	177	T403P (ACC→CCC)	yjjI	DUF3029 family protein, putative glycine radical enzyme
Reads supporting (aligned to +/- strand): ref base T (0/2); new base G (95/78); total (95/82)

Only 100 of 180 total aligned reads displayed.

GCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTT  >  NC_000913/4615762‑4615954

  cTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATc                                                                                             >  1:410319/1‑100 (MQ=255)

   tGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCg                                                                                            <  2:1619779/100‑1 (MQ=255)

   tGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCg                                                                                            <  1:1249676/100‑1 (MQ=255)

      gggTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGAt                                                                                         <  1:3665873/100‑1 (MQ=255)

       ggTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGggatcggatc                                                                                        <  2:4016196/100‑1 (MQ=255)

         ttACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTg                                                                                      >  2:665153/1‑100 (MQ=255)

           aCGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGc                                                                                    <  1:3891561/100‑1 (MQ=255)

           aCGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGc                                                                                    <  2:1317205/100‑1 (MQ=255)

           aCGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGc                                                                                    <  1:196183/100‑1 (MQ=255)

           aCGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCACGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTgtc                                                                                    <  1:2804774/100‑3 (MQ=255)

            cGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCt                                                                                   >  1:3151619/1‑100 (MQ=255)

            cGTTTGATGGTTTCGTCGAGCGACATAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGGGACAGTTTGCAGATGGGGGATCGGATCTGGc                                                                                    >  2:1780118/1‑99 (MQ=255)

             gTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTc                                                                                  >  1:3478473/1‑100 (MQ=255)

                 gATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGCTGCGGCACGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCg                                                                              >  2:1559334/1‑100 (MQ=255)

                  aTGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGc                                                                             <  2:704557/100‑1 (MQ=255)

                   tGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAACAATAAGCATGATGCGGTGCGACAGGTTGCAGATGGGGGATCGGATCTGGCTCATCGc                                                                             >  2:1191316/1‑99 (MQ=255)

                     gTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCAt                                                                          >  2:1035712/1‑100 (MQ=255)

                       ttCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATAc                                                                        <  2:3982688/100‑1 (MQ=255)

                          gtcgAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGc                                                                     >  1:3622584/1‑100 (MQ=255)

                          gtcgAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGc                                                                     <  1:1072004/100‑1 (MQ=255)

                            cgAGCGTCAGAATTTCGCTGATTCTGGAATAATAATAAGCGTGATCCGGGGCTACAGTTTGCAGATGGGGGTTCGGATCTGGCTCATCGCCATACgg                                                                      <  2:2021684/97‑1 (MQ=255)

                            cgAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCaa                                                                   >  2:3136369/1‑100 (MQ=255)

                             gAGCGTCAGAACGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCaaa                                                                  >  2:3169913/1‑100 (MQ=255)

                              aGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAAc                                                                 >  2:1490986/1‑100 (MQ=255)

                                cGTCAGAATGTCTCTGATGCAGGAAAAATAATAAGCATGATGCGGGGCGACAGTTTGCACATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAAcgc                                                               >  2:1189119/1‑100 (MQ=255)

                                 gTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAAcgcg                                                              >  2:2366403/1‑100 (MQ=255)

                                  tCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAAcgcgc                                                             >  2:2136935/1‑100 (MQ=255)

                                  tCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAAcgcgc                                                             >  2:625702/1‑100 (MQ=255)

                                  tCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAAcgcgc                                                             >  1:1912115/1‑100 (MQ=255)

                                   cAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTGGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAAcgcgcg                                                            <  2:1552130/100‑1 (MQ=255)

                                     gAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGcc                                                          >  2:791706/1‑100 (MQ=255)

                                     gAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGcc                                                          >  1:3929021/1‑100 (MQ=255)

                                     gAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGcc                                                          >  1:3357807/1‑100 (MQ=255)

                                         gTCGCTGCTGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGCTGTGGGATAGGATCTCGCTCATCGCCCTACGGCAAGCGCGCGCGCgg                                                       >  2:455478/1‑99 (MQ=255)

                                          tCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCg                                                     >  1:2017278/1‑100 (MQ=255)

                                           cGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCgt                                                    <  2:2145161/100‑1 (MQ=255)

                                            gCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCgtg                                                   >  1:2932213/1‑100 (MQ=255)

                                             cTGATGCCGGAATGATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCgtg                                                   <  2:2772991/99‑1 (MQ=255)

                                               gATGCCGGAATAATAATAAGCATGATGCGCGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCGTCGCCATACGGCAAACGCGCGCCCGGCgtggtg                                                <  2:2490601/100‑1 (MQ=255)

                                                aTGCCGGAATAATAGTAAGCATGATGCGGGGCGACAGTTTGCAGATGGGAGATCGGAGCTGGCTCTGCGCCATACGGCAAACGCGCGCCCGGCGTGGTGc                                               <  2:4051910/100‑1 (MQ=255)

                                                aTGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGc                                               >  2:50760/1‑100 (MQ=255)

                                                aTGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGc                                               >  1:933197/1‑100 (MQ=255)

                                                 tGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGcc                                              <  1:2483815/100‑1 (MQ=255)

                                                 tGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGcc                                              >  1:757750/1‑100 (MQ=255)

                                                  gCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCg                                             <  1:3313869/100‑1 (MQ=255)

                                                   ccGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGa                                            <  2:1184788/100‑1 (MQ=255)

                                                   ccGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATAGCCCTACGGCAAACGCGCGCCCGGCGTGGTGCAGa                                            >  2:173385/1‑100 (MQ=255)

                                                     ggAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCg                                             >  1:1448183/1‑97 (MQ=255)

                                                      gAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGatat                                         >  1:1312075/1‑100 (MQ=255)

                                                      gAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGatat                                         >  1:358174/1‑100 (MQ=255)

                                                       aataataatAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATc                                        >  1:816079/1‑100 (MQ=255)

                                                       aataataatAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATc                                        >  1:376922/1‑100 (MQ=255)

                                                        tctataataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCg                                       <  2:653729/97‑1 (MQ=255)

                                                        ataataataAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCg                                       >  1:2877597/1‑100 (MQ=255)

                                                        ataataataAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCg                                       <  2:2932197/100‑1 (MQ=255)

                                                         taataataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCgg                                      <  1:40842/100‑1 (MQ=255)

                                                          aataataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGa                                     <  2:478429/100‑1 (MQ=255)

                                                           ataataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCgaga                                    >  1:653735/1‑97 (MQ=255)

                                                           ataataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGaa                                    >  1:641387/1‑100 (MQ=255)

                                                           ataataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGaa                                    >  1:3240667/1‑100 (MQ=255)

                                                            taataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAAc                                   >  2:2313050/1‑100 (MQ=255)

                                                            taataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAAc                                   <  1:546501/100‑1 (MQ=255)

                                                             aataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACt                                  <  1:251602/100‑1 (MQ=255)

                                                               taaGCATGATGCTGCGCGCCAGTTTGCAGATGGGTGATCGGCTCTGGCTCATCGCAATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCTGAACTGa                                <  2:2172093/100‑1 (MQ=255)

                                                               taaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGa                                >  2:1394406/1‑100 (MQ=255)

                                                                aaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGAt                               >  1:1786385/1‑100 (MQ=255)

                                                                aaGCAGGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCATGGTGCCGATATCGGAACTGAt                               >  1:3187248/1‑100 (MQ=255)

                                                                 aGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATc                              <  1:1413770/100‑1 (MQ=255)

                                                                  gCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATcc                             >  1:3069352/1‑100 (MQ=255)

                                                                     tgatgCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCg                          <  1:1683439/100‑1 (MQ=255)

                                                                     tgatgCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCg                          >  1:2980018/1‑100 (MQ=255)

                                                                     tgatgCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCg                          <  1:2845279/100‑1 (MQ=255)

                                                                     tgatgCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCg                          >  1:1212526/1‑100 (MQ=255)

                                                                      gatgCGGGGCGACAGTTTGCAGATGGGGGATCGGATCGGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGa                         >  2:2216298/1‑100 (MQ=255)

                                                                        tgcgccgCGACAGTTTGCAGATGGCGGATCTGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCTGCGTGGTGCCGATATCGGAACTGATCCCCGACt                       <  2:2667537/94‑1 (MQ=255)

                                                                        tgCGAGGCGACAGTTTGCAGATGGGGGATAGGATCTGGCTCATCGCGATACGGCAAACGAGTGCCCAGAGTGGTGGCGATATCGGAACTGGTTCCCGAc                        >  2:3551137/1‑99 (MQ=255)

                                                                         gCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACtg                      <  1:868228/100‑1 (MQ=255)

                                                                            gggCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGc                   >  1:107509/1‑100 (MQ=255)

                                                                             ggCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCg                  <  2:4009830/100‑1 (MQ=255)

                                                                               cGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCgtg                <  1:1803370/100‑1 (MQ=255)

                                                                               cGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCgtg                <  2:3366448/100‑1 (MQ=255)

                                                                                gACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCgtgt               >  2:2006957/1‑100 (MQ=255)

                                                                                gACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCgtg                <  2:878601/99‑1 (MQ=255)

                                                                                gACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCACCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCgtgt               >  1:1088701/1‑100 (MQ=255)

                                                                                 aCAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTa              <  1:2593428/100‑1 (MQ=255)

                                                                                  cAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGATTGTGCGTGTa              <  1:3542975/99‑1 (MQ=255)

                                                                                  cAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTaa             >  1:2495128/1‑100 (MQ=255)

                                                                                  cAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTaa             <  1:1002764/100‑1 (MQ=255)

                                                                                    gTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACa           >  1:2434158/1‑100 (MQ=255)

                                                                                    gTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACa           <  2:732136/100‑1 (MQ=255)

                                                                                    gTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCATGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACa           >  2:2893306/1‑100 (MQ=255)

                                                                                     tttGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACAt          >  1:3036488/1‑100 (MQ=255)

                                                                                       tggagagGGGGGAGCGGATCTGGTCCGGCGGCATACGGCAAACGCGCGCCCGTCGTGGTGCCGATATCAGAACTGATCCCCGACTGTGCGTGTa              <  2:3690526/91‑1 (MQ=255)

                                                                                       tGCAGAGGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATgg        <  2:1116063/100‑1 (MQ=255)

                                                                                        gCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATgg        >  2:2027928/1‑99 (MQ=255)

                                                                                         cAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGgcg      <  2:933191/100‑1 (MQ=255)

                                                                                           gATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGcgcg    >  1:631679/1‑100 (MQ=255)

                                                                                           gATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGcgcg    >  1:501218/1‑100 (MQ=255)

                                                                                            aTGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGt   <  1:982929/100‑1 (MQ=255)

                                                                                             tGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGtt  <  2:3330521/100‑1 (MQ=255)

GCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTT  >  NC_000913/4615762‑4615954

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