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.4%	993.5 / ‑2.0	321	coding (962/975 nt)	ygiRQ	Fe‑S oxidoreductase
Reads supporting (aligned to +/- strand): ref base . (0/2); new base C (147/171); total (148/173)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 361 total aligned reads displayed.

ACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAG‑CTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAACA  >  NC_000913/3158109‑3158288

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

 cGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAG‑c                                                                                        >  2:731960/1‑93 (MQ=255)

  gCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAG‑cc                                                                                       <  2:171511/93‑2 (MQ=255)

   cTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAG‑cct                                                                                      <  2:1427307/93‑3 (MQ=255)

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

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

         aTCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGggttg                                                                                 >  1:7657762/1‑93 (MQ=255)

          tctcTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGACTGggttgg                                                                                >  2:391922/1‑93 (MQ=255)

          tctcTTCGGTAGAAAGCGGAATACCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGggttgg                                                                                <  2:4310262/93‑1 (MQ=255)

           ctctTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGggttggt                                                                               <  2:5430567/93‑1 (MQ=255)

            tctTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGggttggtt                                                                              >  1:3431326/1‑93 (MQ=255)

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

              ttCGGTAGAAAGCGAAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTc                                                                            <  2:4560101/93‑1 (MQ=255)

               tCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCg                                                                           <  2:3485857/93‑1 (MQ=255)

                cGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCgt                                                                          <  1:4987825/93‑1 (MQ=255)

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

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

                  gTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCgtgg                                                                        >  2:3716603/1‑93 (MQ=255)

                  gTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCgtgg                                                                        <  1:7891541/93‑1 (MQ=255)

                   tAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCgtggt                                                                       >  2:3354940/1‑93 (MQ=255)

                     gAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGc                                                                     <  1:6271595/93‑1 (MQ=255)

                     gAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGc                                                                     <  1:4323307/93‑1 (MQ=255)

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

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

                         gCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGaa                                                                 <  1:3713672/93‑1 (MQ=255)

                         gCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGaa                                                                 <  2:6711797/93‑1 (MQ=255)

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

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

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

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

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

                            gAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATAc                                                              >  1:3610538/1‑93 (MQ=255)

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

                               tAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATAcgcg                                                           >  1:1040011/1‑93 (MQ=255)

                               tAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATAcgcg                                                           >  2:3282038/1‑93 (MQ=255)

                               tAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATAcgcg                                                           <  1:2225404/93‑1 (MQ=255)

                               tAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATAcgcg                                                           >  1:502347/1‑93 (MQ=255)

                                 gCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGaa                                                         <  1:7792755/93‑1 (MQ=255)

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

                                   aGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGc                                                       >  1:7729518/1‑93 (MQ=255)

                                   aGGCGGGTTGATCCACAAATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGc                                                       <  1:1284281/93‑1 (MQ=255)

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

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

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

                                    ggCGGGGTGATCCACACATAACGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGATGGTTTCGTGGGGCGGAATACGCGAAGCa                                                      >  2:5584793/1‑93 (MQ=255)

                                     gCGGGTTGATCTACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCAt                                                     <  2:7016439/93‑1 (MQ=255)

                                     gCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCAt                                                     <  1:1514906/93‑1 (MQ=255)

                                     gCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCAt                                                     >  2:3000206/1‑93 (MQ=255)

                                     gCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCAt                                                     <  1:5372788/93‑1 (MQ=255)

                                     gCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACACGAAGCAt                                                     <  2:7582577/93‑1 (MQ=255)

                                      cGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATg                                                    >  2:5929874/1‑93 (MQ=255)

                                       gggTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGGGCAGAATACGCGCAGCATgg                                                   >  1:5563615/1‑93 (MQ=255)

                                        ggTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATggg                                                  >  2:101625/1‑93 (MQ=255)

                                        ggTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATggg                                                  >  2:3634823/1‑93 (MQ=255)

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

                                          ttGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCg                                                <  1:4106055/93‑1 (MQ=255)

                                          ttGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCg                                                >  2:6919768/1‑93 (MQ=255)

                                          ttGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCg                                                >  1:2895114/1‑93 (MQ=255)

                                             aTCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGGGAAGCAGGGGCGTAc                                             >  1:2683085/1‑93 (MQ=255)

                                             aTCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTAc                                             >  2:1256797/1‑93 (MQ=255)

                                             aTCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGAATGGGCGTaa                                             >  1:5800972/1‑92 (MQ=255)

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

                                              tCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACa                                            >  1:6515428/1‑93 (MQ=255)

                                               ccACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAg                                           <  1:7341472/93‑1 (MQ=255)

                                               ccACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAg                                           >  1:6007075/1‑93 (MQ=255)

                                               ccACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGGGCAGAATACGCGAAGCATGGCCGtcaag                                           >  1:2539337/1‑89 (MQ=255)

                                                cacacaTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGc                                          >  2:773224/1‑93 (MQ=255)

                                                  cacaTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCAc                                        <  1:4647306/93‑1 (MQ=255)

                                                  cacaTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCAc                                        <  2:4152959/93‑1 (MQ=255)

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

                                                     aTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGAGTACAGCACttt                                     >  2:4065157/1‑93 (MQ=255)

                                                      tAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTa                                    <  1:6872741/93‑1 (MQ=255)

                                                      tAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTa                                    <  1:6974046/93‑1 (MQ=255)

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

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

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

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

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

                                                             cGCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCt                             <  2:2469005/93‑1 (MQ=255)

                                                              gCCGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGGTTCGTGGTGCAGAATACGCGAAGCATGGTCGTACAGCACTTTATCGCCCtt                            >  2:849572/1‑93 (MQ=255)

                                                               ccGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTc                           >  1:1774477/1‑93 (MQ=255)

                                                               ccGTGTTTTTGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTc                           >  2:1492459/1‑93 (MQ=255)

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

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

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

                                                                    tttttGCATCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTGCAGCACTTTATCGCCCTTCACttt                      >  2:3207395/1‑93 (MQ=255)

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

                                                                          cATCAATGCGCGGGCACAGCCTGGTTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGaa                <  2:5612484/93‑1 (MQ=255)

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

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

                                                                            tCAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAaag              >  2:2373789/1‑93 (MQ=255)

                                                                             cAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAaaga             <  1:1606902/93‑1 (MQ=255)

                                                                             cAATGCGCGGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAaaga             <  2:644406/93‑1 (MQ=255)

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

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

                                                                                  cgcgGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGgca        >  2:5704920/1‑93 (MQ=255)

                                                                                   gcgGGCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGgcag       >  2:6248984/1‑93 (MQ=255)

                                                                                      ggCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCaa    >  1:2301149/1‑93 (MQ=255)

                                                                                      ggCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCaa    <  1:3626673/93‑1 (MQ=255)

                                                                                      ggCACAGCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCaa    <  1:1143085/93‑1 (MQ=255)

ACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAG‑CTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAACA  >  NC_000913/3158109‑3158288

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