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.2%	906.4 / ‑1.3	265	coding (962/975 nt)	ygiRQ	Fe‑S oxidoreductase
Reads supporting (aligned to +/- strand): ref base . (0/2); new base C (106/156); total (107/158)
Fisher's exact test for biased strand distribution p-value = 5.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.50e-01

Only 100 of 294 total aligned reads displayed.

ACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAG‑C‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAACA  >  NC_000913/3158109‑3158288

aCGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATACACACATAGCGGTTGCCGTGTTTTTGCATCAACGGGCGGGCACAg                                                                                           >  1:5758348/1‑93 (MQ=255)

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

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

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

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

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

    tGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGTGCGGGCACAG‑cctg                                                                                      >  1:831842/1‑90 (MQ=255)

       ccATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGggt                                                                                   <  2:3563914/93‑1 (MQ=255)

       ccATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGggt                                                                                   <  1:5317434/93‑1 (MQ=255)

       ccATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGggt                                                                                   >  1:5371349/1‑93 (MQ=255)

       ccATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGGTGATCCACACATAGCGGACGCCGTGTGTTTGCATCAATGCGCGGGCACAGCC‑TGggt                                                                                   >  2:2697575/1‑93 (MQ=255)

        cATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGggtt                                                                                  >  1:3173776/1‑93 (MQ=255)

        cATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGggtt                                                                                  >  1:6253844/1‑93 (MQ=255)

         aTCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGggttg                                                                                 >  2:3022199/1‑93 (MQ=255)

          tctcTTCGGTAGAAAGCGGAATAGCAGGCGGGCTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGggttgg                                                                                >  2:5981176/1‑93 (MQ=255)

           ctctTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGggttggt                                                                               <  1:5391665/93‑1 (MQ=255)

           ctctTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGggttggt                                                                               <  2:4957543/93‑1 (MQ=255)

           ctctTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGggttggt                                                                               <  1:4879938/93‑1 (MQ=255)

           ctctTCGGTAGAAAGCGGAATAGCAGGCGGGGTGATCCACACATAGCGGGCTCCGTTTTGTTGCATCAATGCGCGGGCACAGCC‑TGggttggt                                                                               <  1:1979152/93‑1 (MQ=255)

            tctTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGggttggtt                                                                              >  2:5808798/1‑93 (MQ=255)

            tctTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGCGTTTTTGCATCACTGCGCGGGCACAGCC‑TGggttggat                                                                              >  1:2624075/1‑91 (MQ=255)

            tctTCGGTAGAAAGCGGAATAGCAGGCGGGGTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGggttggtt                                                                              >  2:6266723/1‑93 (MQ=255)

             ctTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGttt                                                                             <  1:5593954/93‑1 (MQ=255)

             ctTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGttt                                                                             >  2:873368/1‑93 (MQ=255)

             ctTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGttt                                                                             <  2:4416619/93‑1 (MQ=255)

             ctTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGttt                                                                             >  2:1025892/1‑93 (MQ=255)

              ttCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTc                                                                            >  2:5619145/1‑93 (MQ=255)

               tCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCg                                                                           >  2:5283307/1‑93 (MQ=255)

                 ggTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCgtg                                                                         <  2:4991636/93‑1 (MQ=255)

                 ggTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCgtg                                                                         <  2:5429934/93‑1 (MQ=255)

                  gTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCgtgg                                                                        <  1:1893733/93‑1 (MQ=255)

                    aGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCATAGCC‑TGGGTTGGTTTCgtggtg                                                                      >  1:4998453/1‑93 (MQ=255)

                      aaaGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCa                                                                    >  1:4326243/1‑93 (MQ=255)

                      aaaGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCa                                                                    >  2:3904605/1‑93 (MQ=255)

                      aaaGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCa                                                                    <  1:1122538/93‑1 (MQ=255)

                       aaGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAg                                                                   <  2:1700856/93‑1 (MQ=255)

                        aGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGa                                                                  >  1:3391728/1‑93 (MQ=255)

                         gCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGaa                                                                 <  2:4296462/93‑1 (MQ=255)

                         gCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGaa                                                                 <  2:4692139/93‑1 (MQ=255)

                          cGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAAt                                                                <  2:4877427/93‑1 (MQ=255)

                           ggAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATa                                                               >  2:1495237/1‑93 (MQ=255)

                           ggAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATa                                                               <  1:2744585/93‑1 (MQ=255)

                            gAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATAc                                                              <  1:1751625/93‑1 (MQ=255)

                             aaTAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATAcg                                                             >  2:3224913/1‑93 (MQ=255)

                              aTAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATAcgc                                                            <  2:408388/93‑1 (MQ=255)

                              aTAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATAcgc                                                            <  2:2081765/93‑1 (MQ=255)

                                  cAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAg                                                        <  1:3496288/93‑1 (MQ=255)

                                   aGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGc                                                       >  1:5990570/1‑93 (MQ=255)

                                   aGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGc                                                       >  2:2618268/1‑93 (MQ=255)

                                    ggCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCa                                                      >  2:3240636/1‑93 (MQ=255)

                                    ggCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCa                                                      <  1:2285922/93‑1 (MQ=255)

                                     gCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCAt                                                     <  2:1270228/93‑1 (MQ=255)

                                     gCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCAt                                                     <  1:4823439/93‑1 (MQ=255)

                                       gggTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATgg                                                   >  1:6331850/1‑93 (MQ=255)

                                        ggTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATggg                                                  >  2:1346508/1‑93 (MQ=255)

                                        ggTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATggg                                                  <  2:5650481/93‑1 (MQ=255)

                                         gtgtaTCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGc                                                 <  1:4709169/89‑1 (MQ=255)

                                         gTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATgggg                                                 >  1:1327462/1‑92 (MQ=255)

                                           tGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGt                                               <  2:4336394/93‑1 (MQ=255)

                                            gATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTa                                              <  1:2814699/93‑1 (MQ=255)

                                            gATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTa                                              <  1:2971329/93‑1 (MQ=255)

                                            gATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTa                                              >  2:4603732/1‑93 (MQ=255)

                                            aaTCCACACATAGCGGTCTCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTa                                              <  2:5573667/92‑1 (MQ=255)

                                             aTCCACACGTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTAc                                             <  2:5034911/93‑1 (MQ=255)

                                             aTCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTAc                                             >  1:1518785/1‑93 (MQ=255)

                                                  cacaTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCAc                                        <  1:1707164/93‑1 (MQ=255)

                                                   acaTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACt                                       <  1:6077597/93‑1 (MQ=255)

                                                    caTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACtt                                      >  2:5610317/1‑93 (MQ=255)

                                                     aTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACttt                                     >  2:6105086/1‑93 (MQ=255)

                                                      tAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTa                                    <  1:732943/93‑1 (MQ=255)

                                                       aGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTAt                                   <  1:709531/93‑1 (MQ=255)

                                                         cGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCg                                 <  2:2755758/93‑1 (MQ=255)

                                                            tCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGccc                              <  1:3885851/93‑1 (MQ=255)

                                                            tCGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGccc                              <  2:5110320/93‑1 (MQ=255)

                                                             cGCCGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCt                             <  1:4097641/93‑1 (MQ=255)

                                                               ccGTGTTTTTGCATCAATGCGCGGTCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTc                           <  2:1992066/93‑1 (MQ=255)

                                                               ccGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGCTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTc                           >  1:2477302/1‑93 (MQ=255)

                                                                cGTGTTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCa                          >  1:2582534/1‑93 (MQ=255)

                                                                  tgtTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACt                        <  1:3023161/93‑1 (MQ=255)

                                                                   gtTTTTGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACtt                       <  2:2343592/93‑1 (MQ=255)

                                                                     ttttGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTtc                     <  1:1411802/93‑1 (MQ=255)

                                                                      tttGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTtct                    <  1:3535193/93‑1 (MQ=255)

                                                                       ttGCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTtctc                   <  1:5457092/93‑1 (MQ=255)

                                                                         gCATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGa                 <  2:3348158/93‑1 (MQ=255)

                                                                          cATCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGaa                >  2:3776650/1‑93 (MQ=255)

                                                                           aTCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGaaa               >  1:6452283/1‑93 (MQ=255)

                                                                            tCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAaag              >  1:3092627/1‑93 (MQ=255)

                                                                            tCAATGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAaag              >  2:2428855/1‑93 (MQ=255)

                                                                              aaTGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAaagaa            >  2:1835505/1‑93 (MQ=255)

                                                                               aTGCGCGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAaagaag           <  2:6393007/93‑1 (MQ=255)

                                                                                 gcgcGGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAggt         >  1:3617623/1‑92 (MQ=255)

                                                                                  cgcgGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAggta        <  2:1767896/93‑3 (MQ=255)

                                                                                   gcgGGCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGgtag       <  1:5428960/93‑1 (MQ=255)

                                                                                    cgGGCACAGCC‑TGGGTTGGCTTCGTGGTGCAGAATACGGGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGgtagc      <  1:5230423/93‑1 (MQ=255)

                                                                                     gggCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGgtagca     >  1:3679295/1‑93 (MQ=255)

                                                                                     gggCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGgtagca     >  1:4487847/1‑93 (MQ=255)

                                                                                      ggCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCaa    <  2:4429859/93‑1 (MQ=255)

                                                                                      ggCACAG‑CTTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCaa    <  1:4558286/93‑1 (MQ=255)

                                                                                       gCACAGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAac   <  2:3525161/93‑1 (MQ=255)

                                                                                        cacaGCC‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGTAGCAaca  <  2:2011349/93‑1 (MQ=255)

ACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACAG‑C‑TGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAACA  >  NC_000913/3158109‑3158288

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