BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,158,202	+C	coding (961/975 nt)	ygiRQ ←	Fe‑S oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,158,201	1	.	C	99.5%	683.4 / ‑4.1	216	coding (962/975 nt)	ygiRQ	Fe‑S oxidoreductase
Reads supporting (aligned to +/- strand): ref base . (1/0); new base C (90/125); total (91/125)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 238 total aligned reads displayed.

ACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAG‑CTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCA  >  NC_000913/3158109‑3158285

 cGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAG‑c                                                                                     <  1:3856665/93‑1 (MQ=255)

   cTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAG‑cct                                                                                   <  1:3931040/93‑3 (MQ=255)

    tGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAG‑cctg                                                                                  <  1:1098584/93‑4 (MQ=255)

     gTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAG‑cctgg                                                                                 <  2:5334607/93‑5 (MQ=255)

       ccATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGggt                                                                                <  2:3585599/93‑1 (MQ=255)

        cATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGggtt                                                                               <  2:5521582/93‑1 (MQ=255)

        cATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGggtt                                                                               >  2:3358404/1‑93 (MQ=255)

         aTCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGggttg                                                                              <  1:279551/93‑1 (MQ=255)

           ctctTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGggttggt                                                                            >  2:1804839/1‑93 (MQ=255)

             ctTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGttt                                                                          >  2:413783/1‑93 (MQ=255)

             ctTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGttt                                                                          <  2:2939929/93‑1 (MQ=255)

             ctTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCACCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGttt                                                                          >  1:1307825/1‑93 (MQ=255)

              ttCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTc                                                                         >  2:4343974/1‑93 (MQ=255)

                 ggTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCgtg                                                                      >  2:1301304/1‑93 (MQ=255)

                 ggTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCgtg                                                                      >  1:3335816/1‑93 (MQ=255)

                  gTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCgggg                                                                     >  2:2703215/1‑93 (MQ=255)

                   tAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCgtggt                                                                    <  1:1948736/93‑1 (MQ=255)

                   tAGAAAGCGGAATAGCAGGCGCGGTGATCCACACATAACGGTCGCCGTGTTTTTGCATCAATTCGCGCGCACAGCCTGGGTTGGTTTCgtggg                                                                    >  1:5670181/1‑92 (MQ=255)

                    aGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCgtggtg                                                                   <  1:2703821/93‑1 (MQ=255)

                     gAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGc                                                                  >  1:4257281/1‑93 (MQ=255)

                      aaaGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCa                                                                 <  2:4055211/93‑1 (MQ=255)

                      aaaGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCa                                                                 >  1:1280470/1‑93 (MQ=255)

                      aaaGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCa                                                                 >  2:2182149/1‑93 (MQ=255)

                      aaaGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCa                                                                 <  1:350273/93‑1 (MQ=255)

                       aaGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAg                                                                <  1:2963530/93‑1 (MQ=255)

                        aGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGa                                                               >  1:5807583/1‑93 (MQ=255)

                          cGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAAt                                                             <  1:3074491/93‑1 (MQ=255)

                          cGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAAt                                                             >  1:4337722/1‑93 (MQ=255)

                          cGGAATAGCAGGCGGGTTGATCCACACAGAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAAt                                                             <  1:3890476/93‑1 (MQ=255)

                           ggAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATa                                                            >  1:1791776/1‑93 (MQ=255)

                           ggAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATa                                                            <  1:3563885/93‑1 (MQ=255)

                            gAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATAc                                                           <  1:5654071/93‑1 (MQ=255)

                             aaTAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATAcg                                                          >  2:3527060/1‑93 (MQ=255)

                                aGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGa                                                       <  1:2596577/93‑1 (MQ=255)

                                  cAGGCGGAGTCATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAg                                                     <  1:3945383/93‑1 (MQ=255)

                                   aGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGc                                                    <  1:1888708/93‑1 (MQ=255)

                                    ggCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCa                                                   <  1:2360385/93‑1 (MQ=255)

                                    ggCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGGTGGTTTCGTGGTGCAGAATACGCGAAGCa                                                   >  1:4026046/1‑93 (MQ=255)

                                    ggCGGGTTGATCCACACACAGCGGTCGCCCTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGGTGGTTTCGTTGTGCAGAATACACGAAGCa                                                   >  2:722924/1‑93 (MQ=255)

                                      cGGGTTGATCCACACATAGCGGTCGCCGTTTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATg                                                 <  1:3299240/93‑1 (MQ=255)

                                      cGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATg                                                 <  1:2187660/93‑1 (MQ=255)

                                       gggTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATgg                                                <  1:4674603/93‑1 (MQ=255)

                                       gggTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATgg                                                <  1:5775263/93‑1 (MQ=255)

                                       gggTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATgg                                                <  2:3812321/93‑1 (MQ=255)

                                        ggTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGTTTGGTTTCGTGGTGCAGAATACGCGAAGCATggg                                               <  2:4666644/93‑1 (MQ=255)

                                        ggTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATggg                                               <  1:4098087/93‑1 (MQ=255)

                                         gtggATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGc                                              <  2:1469373/90‑1 (MQ=255)

                                         gTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGc                                              >  2:1203721/1‑93 (MQ=255)

                                           tGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGt                                            <  2:4606644/93‑1 (MQ=255)

                                            gATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTa                                           >  1:6274309/1‑93 (MQ=255)

                                              tCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACa                                         <  1:1623251/93‑1 (MQ=255)

                                                 acacaTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGcc                                      >  1:1162552/1‑92 (MQ=255)

                                                   acaTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACt                                    <  2:3727328/93‑1 (MQ=255)

                                                   acaTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACt                                    >  2:4818762/1‑93 (MQ=255)

                                                    caTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACtt                                   <  1:5180130/93‑1 (MQ=255)

                                                    caTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACtt                                   <  2:5577322/93‑1 (MQ=255)

                                                     aTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACttt                                  >  2:2500255/1‑93 (MQ=255)

                                                     aTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACttt                                  <  1:1687653/93‑1 (MQ=255)

                                                     aTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACttt                                  >  1:2312129/1‑93 (MQ=255)

                                                     aTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACttt                                  >  1:5895689/1‑93 (MQ=255)

                                                     aTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACttt                                  <  1:4976290/93‑1 (MQ=255)

                                                      tAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTa                                 >  1:4746342/1‑93 (MQ=255)

                                                      tAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTa                                 >  2:5679998/1‑93 (MQ=255)

                                                       aGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTAt                                <  1:140600/93‑1 (MQ=255)

                                                       aGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTAt                                <  1:1975704/93‑1 (MQ=255)

                                                        gCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATc                               >  1:4516011/1‑93 (MQ=255)

                                                        gCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATc                               <  1:2230447/93‑1 (MQ=255)

                                                         cGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCg                              <  1:2542647/93‑1 (MQ=255)

                                                          ggTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGc                             <  2:5120769/93‑1 (MQ=255)

                                                          ggTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGc                             >  1:5596876/1‑93 (MQ=255)

                                                           gTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGcc                            >  1:1079167/1‑93 (MQ=255)

                                                               ccGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTc                        <  2:4428826/93‑1 (MQ=255)

                                                                cGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCa                       >  2:2443030/1‑93 (MQ=255)

                                                                cGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCa                       >  2:4946972/1‑93 (MQ=255)

                                                                cGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCa                       <  1:827358/93‑1 (MQ=255)

                                                                 gtgtTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCAc                      >  2:1914796/1‑93 (MQ=255)

                                                                  tgtTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACt                     <  1:987144/93‑1 (MQ=255)

                                                                  tgtTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACt                     <  2:3873300/93‑1 (MQ=255)

                                                                   gtTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACtt                    <  1:4198884/93‑1 (MQ=255)

                                                                   gtTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACtt                    >  1:1541094/1‑93 (MQ=255)

                                                                    tttttGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACttt                   >  2:2136582/1‑93 (MQ=255)

                                                                    tttttGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACttt                   >  1:4872426/1‑93 (MQ=255)

                                                                     ttttGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTtc                  <  2:2625892/93‑1 (MQ=255)

                                                                     ttttGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTtc                  >  2:6187144/1‑93 (MQ=255)

                                                                       ttGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTtctc                >  2:4163843/1‑93 (MQ=255)

                                                                        tGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCg               >  2:2215057/1‑93 (MQ=255)

                                                                        tGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCg               <  1:5283892/93‑1 (MQ=255)

                                                                        tGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCg               >  2:313490/1‑93 (MQ=255)

                                                                         gCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGa              <  1:3180993/93‑1 (MQ=255)

                                                                           aTCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGaaa            <  2:3296793/93‑1 (MQ=255)

                                                                              aaTGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAaagaa         <  2:1631776/93‑1 (MQ=255)

                                                                              aaTGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAaagaa         >  1:3524653/1‑93 (MQ=255)

                                                                              aaTGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAaagaa         >  1:1907630/1‑93 (MQ=255)

                                                                               aTGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAaagaag        <  2:5780119/93‑1 (MQ=255)

                                                                                tGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAgg       <  1:6182566/93‑1 (MQ=255)

                                                                                tGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAgg       >  2:5246139/1‑93 (MQ=255)

                                                                                tGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAgg       >  1:4943864/1‑93 (MQ=255)

                                                                                 gcgcGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGgc      <  1:2448275/93‑1 (MQ=255)

                                                                                    cgGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGgcagc   <  1:258004/93‑1 (MQ=255)

                                                                                     gggCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGgcagca  >  1:4366197/1‑93 (MQ=255)

ACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAG‑CTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCA  >  NC_000913/3158109‑3158285

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