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	98.9%	600.2 / ‑1.4	187	coding (962/975 nt)	ygiRQ	Fe‑S oxidoreductase
Reads supporting (aligned to +/- strand): ref base . (0/2); new base C (92/93); total (92/95)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 206 total aligned reads displayed.

ACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑G‑CTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAAC  >  NC_000913/3158109‑3158287

aCGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑g                                                                                         <  1:76527/93‑1 (MQ=255)

  gCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑G‑cc                                                                                      <  2:2280356/93‑2 (MQ=255)

  gCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑G‑cc                                                                                      <  1:3065164/93‑2 (MQ=255)

  gCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑G‑cc                                                                                      <  2:3286591/93‑2 (MQ=255)

    tGTCCATCTCTACGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑G‑cctg                                                                                    <  2:3265021/93‑4 (MQ=255)

     gTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑G‑cctgg                                                                                   <  2:2123199/93‑5 (MQ=255)

       ccATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGggt                                                                                  <  2:3852191/93‑1 (MQ=255)

            tctTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGggttggtt                                                                             <  2:2619770/93‑1 (MQ=255)

            tctTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGggttggtt                                                                             >  2:3770089/1‑93 (MQ=255)

              ttCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTc                                                                           <  1:1321595/93‑1 (MQ=255)

               tCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTct                                                                          >  1:4697263/1‑92 (MQ=255)

               tCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCg                                                                          >  1:1697197/1‑93 (MQ=255)

               tCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCg                                                                          <  1:2243693/93‑1 (MQ=255)

                 ggTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCgtg                                                                        >  2:3188679/1‑93 (MQ=255)

                 ggTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCgtg                                                                        >  1:4364663/1‑93 (MQ=255)

                  gTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCgtgg                                                                       >  1:4076330/1‑93 (MQ=255)

                   tAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCgtggt                                                                      <  1:4298964/93‑1 (MQ=255)

                     gAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGc                                                                    >  2:556667/1‑93 (MQ=255)

                     aaaaGCGGAATAGCAGGCGGGTTGAGCCACACAGAGCGGTCGCCGTCTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGc                                                                    <  1:1810366/92‑1 (MQ=255)

                      aaaGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCa                                                                   >  2:3649461/1‑93 (MQ=255)

                      aaaGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCa                                                                   <  2:3813402/93‑1 (MQ=255)

                       aaGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGGTTCGTGGTGCAg                                                                  >  1:2172787/1‑93 (MQ=255)

                         gCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGaa                                                                <  2:179066/93‑1 (MQ=255)

                          cGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAAt                                                               <  2:4046522/93‑1 (MQ=255)

                          cGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAAt                                                               <  2:1098048/93‑1 (MQ=255)

                          cGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAAt                                                               <  2:1298311/93‑1 (MQ=255)

                             aaTAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATAcg                                                            <  1:1133463/93‑1 (MQ=255)

                               tAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATAcgcg                                                          >  2:4016154/1‑93 (MQ=255)

                                 gCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGaa                                                        >  2:2808328/1‑93 (MQ=255)

                                 gCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGaa                                                        >  1:2790595/1‑93 (MQ=255)

                                 gCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGaa                                                        >  1:1719513/1‑93 (MQ=255)

                                    ggCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCa                                                     >  2:1645634/1‑93 (MQ=255)

                                     gCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCAt                                                    <  1:3538306/93‑1 (MQ=255)

                                     gCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGCGGTGCAGAATACGCGAAGCAt                                                    >  1:1440909/1‑93 (MQ=255)

                                      cGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATg                                                   <  2:2862406/93‑1 (MQ=255)

                                      cGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATg                                                   <  1:2722808/93‑1 (MQ=255)

                                       gggTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATgg                                                  <  1:4012336/93‑1 (MQ=255)

                                       gggCTGATCCACACATAGCGGTCGCCGGGTTTTTGCATCAATCCGCAGGCACC‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCactg                                                  >  2:314947/1‑90 (MQ=255)

                                        ggTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATggg                                                 <  2:1982404/93‑1 (MQ=255)

                                        ggTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATggg                                                 <  1:1387890/93‑1 (MQ=255)

                                        ggTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATggg                                                 >  2:363326/1‑93 (MQ=255)

                                         gTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGc                                                >  2:4422157/1‑93 (MQ=255)

                                         gTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGc                                                >  2:598771/1‑93 (MQ=255)

                                           tGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGt                                              <  1:1942106/93‑1 (MQ=255)

                                            gATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTa                                             >  1:4511957/1‑93 (MQ=255)

                                            gATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTa                                             <  1:4439577/93‑1 (MQ=255)

                                             aTCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTAc                                            >  1:4615347/1‑93 (MQ=255)

                                             aTCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTAc                                            >  2:2674631/1‑93 (MQ=255)

                                             aTCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTAc                                            >  1:4157178/1‑93 (MQ=255)

                                             aTCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTAc                                            <  2:1284375/93‑1 (MQ=255)

                                              tCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACa                                           <  1:1470836/93‑1 (MQ=255)

                                              tCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACa                                           >  1:3367044/1‑93 (MQ=255)

                                               ccACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAg                                          >  1:1948338/1‑93 (MQ=255)

                                               ccACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGGTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAg                                          >  2:4233770/1‑93 (MQ=255)

                                                cacacaTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGc                                         >  1:592485/1‑93 (MQ=255)

                                                 acacaTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCa                                        >  1:4197550/1‑93 (MQ=255)

                                                 acacaTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCa                                        >  1:3296178/1‑93 (MQ=255)

                                                 acacaTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGATGGTTTCGTGGTGCAGAATACGCGAAGCATTGGCGTACAGCa                                        >  2:2485415/1‑93 (MQ=255)

                                                 acacaTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGCGTTGGTTTCGTGGTGCAGAATACGCGAAGCGTGGGTGTACAGCa                                        >  2:2150869/1‑93 (MQ=255)

                                                  cacaTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCAc                                       <  2:139471/93‑1 (MQ=255)

                                                   acaTAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACt                                      <  2:4240630/93‑1 (MQ=255)

                                                    caTAGCGGTGGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACtt                                     <  1:3817832/93‑1 (MQ=255)

                                                      tAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTa                                   >  1:3631843/1‑93 (MQ=255)

                                                      tAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTa                                   >  2:2705439/1‑93 (MQ=255)

                                                      tAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTa                                   <  2:440164/93‑1 (MQ=255)

                                                      tAGCGGTCGCCGTGTGTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTa                                   <  2:2424460/93‑1 (MQ=255)

                                                       aGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTAt                                  <  1:493962/93‑1 (MQ=255)

                                                       aGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTAt                                  <  2:2195449/93‑1 (MQ=255)

                                                        gCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATc                                 <  1:2150062/93‑1 (MQ=255)

                                                          ggTCTCTGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCGCTTTATCGc                               <  2:418056/93‑1 (MQ=255)

                                                          ggTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGc                               >  1:3315496/1‑93 (MQ=255)

                                                          ggTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGc                               <  1:2678852/93‑1 (MQ=255)

                                                           gTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGcc                              >  1:2359563/1‑93 (MQ=255)

                                                            tCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGccc                             <  2:469041/93‑1 (MQ=255)

                                                             cGCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCt                            <  1:1127541/93‑1 (MQ=255)

                                                              gCCGTGTTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCtt                           >  1:3945021/1‑93 (MQ=255)

                                                                 gtgtTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCat                        >  1:1615930/1‑92 (MQ=255)

                                                                 gtgtTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCAc                        <  2:297899/93‑1 (MQ=255)

                                                                 gtgtTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCAc                        >  1:2432930/1‑93 (MQ=255)

                                                                 gtgtTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCAc                        <  1:2291072/93‑1 (MQ=255)

                                                                   gtTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACtt                      <  1:3137839/93‑1 (MQ=255)

                                                                   gtTTTTGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACtt                      >  1:1018932/1‑93 (MQ=255)

                                                                     ttttGCATCAATGCGCGGGCGCA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTtc                    >  2:3503145/1‑93 (MQ=255)

                                                                     ttttGCAACAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTtc                    <  2:1042155/93‑1 (MQ=255)

                                                                      tttGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTtct                   <  1:3950731/93‑1 (MQ=255)

                                                                       ttGCATCAATGCGCGGGCACA‑GCCTGGGTTGGTGTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTtctc                  <  1:4569910/93‑1 (MQ=255)

                                                                         gCATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGa                >  1:3335055/1‑93 (MQ=255)

                                                                          cATCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGaa               >  1:691278/1‑93 (MQ=255)

                                                                           aTCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGaaa              <  2:2440621/93‑1 (MQ=255)

                                                                           aTCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGaaa              >  1:2301720/1‑93 (MQ=255)

                                                                           aTCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGaaa              >  1:1245913/1‑93 (MQ=255)

                                                                           aTCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCACGaaa              >  2:3531972/1‑93 (MQ=255)

                                                                            tCAATGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAaag             <  2:1764790/93‑1 (MQ=255)

                                                                              aaTGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAaagaa           <  1:2605142/93‑1 (MQ=255)

                                                                              aaTGCGCGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAaagaa           >  2:4722576/1‑93 (MQ=255)

                                                                                agcgcGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAgg         <  1:3450792/92‑1 (MQ=255)

                                                                                 gcgcGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGgc        <  2:4077694/93‑1 (MQ=255)

                                                                                 gcgcGGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGgc        <  2:6627/93‑1 (MQ=255)

                                                                                  cgcgGGCACATC‑CTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGgca       <  2:1583835/93‑1 (MQ=255)

                                                                                  cgcgGGCACA‑GCCTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGgca       <  2:2341480/93‑1 (MQ=255)

ACGCTGTCCATCTCTTCGGTAGAAAGCGGAATAGCAGGCGGGTTGATCCACACATAGCGGTCGCCGTGTTTTTGCATCAATGCGCGGGCACA‑G‑CTGGGTTGGTTTCGTGGTGCAGAATACGCGAAGCATGGGCGTACAGCACTTTATCGCCCTTCACTTTCTCGAAAGAAGGCAGCAAC  >  NC_000913/3158109‑3158287

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