BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,427,261	T→C	V130A (GTC→GCC)	b3279 →	putative transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,427,261	0	T	C	99.5%	678.5 / NA	212	V130A (GTC→GCC)	b3279	putative transferase
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (91/120); minor base A (1/0); total (92/120)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 219 total aligned reads displayed.

TTGGCAATCGAGTTTTGGTTGGGATGGGCTCAATTTTACTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGTCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCAAACAGATCCGCCCGTTAAGTGATGAAGAGAAGGCTG  >  NC_000913/3427171‑3427350

ttGGCGATGGAGTTTTGGTTGGGATGGGCTCAATTTTACTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGccc                                                                                         <  1:2800504/93‑1 (MQ=255)

  ggCAATCGAGTTTTGGTTGGGATGGGCTCAATTTTACTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCca                                                                                       >  2:5808840/1‑93 (MQ=255)

  ggCAATCGAGTTTTGGTTGGGATGGGCTCAATTTTACTTAATGGCGCAATATTAGAAGATGATGTGATGATAGGTGCGGGTAGTCTGGCCCca                                                                                       <  2:1139043/93‑1 (MQ=255)

   gCAATCGAGTTTTGGTTGGGATGGGCTCAATTTTACTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCcac                                                                                      <  2:4819507/93‑1 (MQ=255)

    cAATCGAGTTTTGGTTGGGATGGGCTCAATTTTACTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCcaca                                                                                     <  1:5016629/93‑1 (MQ=255)

     aaTCGAGTTTTGGTTGGGATGGGCTCAATTTTACTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACaa                                                                                    <  1:6149516/93‑1 (MQ=255)

        cGAGTTTTGGTTGGGATGGGCTCAATTTTACTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAAt                                                                                 <  2:214736/93‑1 (MQ=255)

         gAGTTTTGGTTGGGATGGGCTCAATTTTACTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATa                                                                                <  2:3185826/93‑1 (MQ=255)

         gAGTTTTGGTTGGGATGGGCTCAATTTTACTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATa                                                                                <  1:4230623/93‑1 (MQ=255)

         gAGTTTTGGTTGGGATGGGCTCAATTTTACTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATa                                                                                <  1:592572/93‑1 (MQ=255)

         gAGTTTTGGTTGGGATGGGCTCAATTTTACTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATa                                                                                <  2:2237037/93‑1 (MQ=255)

          aGTTTTGGTTGGGATGGGCTCAATTTTACTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATaa                                                                               >  2:2189650/1‑93 (MQ=255)

           gTTTTGGTGGGGATGGGCTCAATTTTACTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATaaa                                                                              <  2:1544615/93‑1 (MQ=255)

             tttGGTTGGGATGGGCTCAATTTTACTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACg                                                                            <  1:3875274/93‑1 (MQ=255)

              ttggttggGATGGGCTCAATTTTACTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACgg                                                                           >  1:5705008/1‑93 (MQ=255)

               tggttggGATGGGCTCAATTTTACTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGc                                                                          >  1:4428898/1‑93 (MQ=255)

               tggttggGATGGGCTCAATTTTACTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGc                                                                          <  1:3120369/93‑1 (MQ=255)

                ggttggTATGTGCTCAATTTTACTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCt                                                                         <  2:5609645/93‑1 (MQ=255)

                ggttggGATGGGCTCAATTTTACTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCt                                                                         <  1:2109434/93‑1 (MQ=255)

                  ttggGATGGGCTCAATTTTACTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTgg                                                                       <  1:192667/93‑1 (MQ=255)

                  ttggGATGGGCTCAATTTTACTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTgg                                                                       >  1:3799344/1‑93 (MQ=255)

                      gATGGGCTCAATTTTACTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGagag                                                                   <  1:1641988/93‑1 (MQ=255)

                       aTGGGCTCAATTTTACTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGc                                                                  >  1:1713942/1‑93 (MQ=255)

                        tGGGCTCAATTTTACTTGATGACGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCg                                                                 >  2:4008770/1‑93 (MQ=255)

                         gggCTCAATTTTACTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCgg                                                                <  2:2360232/93‑1 (MQ=255)

                         gggCTCAATTTTACTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCgg                                                                >  2:3707551/1‑93 (MQ=255)

                          ggCTCAATTTTACTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGa                                                               >  2:4419097/1‑93 (MQ=255)

                            cTCAATTTTACTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGata                                                             >  2:4605295/1‑93 (MQ=255)

                            cTCAATTTTACTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGata                                                             >  1:2914683/1‑93 (MQ=255)

                              cAATTTTACTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATc                                                           >  2:4680820/1‑93 (MQ=255)

                              cAATTTTACTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATc                                                           >  1:4560277/1‑93 (MQ=255)

                              cAATTTTACTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATc                                                           <  1:4624320/93‑1 (MQ=255)

                              cAATTTTACTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATc                                                           <  2:1518579/93‑1 (MQ=255)

                               aaTTTTACTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCt                                                          <  1:3421294/93‑1 (MQ=255)

                               aaTTTTACTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCt                                                          >  1:1358265/1‑93 (MQ=255)

                                aTTTTACTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTg                                                         >  1:3268443/1‑93 (MQ=255)

                                aTTTTACTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTg                                                         <  2:5424732/93‑1 (MQ=255)

                                 ttttACTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGt                                                        <  2:5730883/93‑1 (MQ=255)

                                  tttACTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTa                                                       <  1:198806/93‑1 (MQ=255)

                                    tACTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTAtc                                                     >  1:4273370/1‑93 (MQ=255)

                                    tACTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTAtc                                                     >  2:4781471/1‑93 (MQ=255)

                                     aCTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTAtct                                                    <  1:1516021/93‑1 (MQ=255)

                                     aCTTGATGGCGCAATAGTAGAAGATGATGCGATGATTGGAGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTAtcc                                                    >  1:1026268/1‑92 (MQ=255)

                                      cTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTAtctc                                                   >  2:2413959/1‑93 (MQ=255)

                                      cTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTAtctc                                                   >  1:5652543/1‑93 (MQ=255)

                                      cTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTAtctc                                                   <  2:5532663/93‑1 (MQ=255)

                                       ttGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCg                                                  >  2:1415923/1‑93 (MQ=255)

                                       ttGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCg                                                  <  1:3409349/93‑1 (MQ=255)

                                          aTGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTa                                               >  2:3534335/1‑93 (MQ=255)

                                              cgcAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAACCGGCTGGAGAGCGGATATCTGTATCTCGGTAGccc                                           <  1:4095913/93‑1 (MQ=255)

                                               gcAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGcccc                                          <  2:5710187/93‑1 (MQ=255)

                                               gcAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGcccc                                          <  2:2554825/93‑1 (MQ=255)

                                                cAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCg                                         >  1:2255146/1‑93 (MQ=255)

                                                cAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCg                                         >  2:4258718/1‑93 (MQ=255)

                                                 aaTAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGt                                        >  1:5517635/1‑93 (MQ=255)

                                                 aaTAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGt                                        <  2:2292695/93‑1 (MQ=255)

                                                   tagtagAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCa                                      >  2:5194887/1‑93 (MQ=255)

                                                   tagtagAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCa                                      >  1:3843791/1‑93 (MQ=255)

                                                      tagAAGATGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCAAAc                                   >  1:55916/1‑93 (MQ=255)

                                                       agaagaTGATGTGATGATTGCTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCAAACa                                  >  2:6065761/1‑93 (MQ=255)

                                                        gaagaTGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAACAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCAAACAg                                 >  1:1722466/1‑93 (MQ=255)

                                                         aagaTGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCAAACAGa                                >  1:5914882/1‑93 (MQ=255)

                                                         aagaTGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCAAACAGa                                <  2:6140913/93‑1 (MQ=255)

                                                         aagaTGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCAAACAGa                                <  1:2134107/93‑1 (MQ=255)

                                                         aagaTGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCAAACAGa                                <  1:3023552/93‑1 (MQ=255)

                                                          agaTGATGTGATGATTGGTGCGGGTAGTCTGGCCCCAGAAAATAAACGGCTGGAAAGCGGATATCTGTATCTCGGTAGCCCCGTCAAACAGAt                               <  2:3153608/93‑1 (MQ=255)

                                                          agaTGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCAAACAGAt                               <  2:558279/93‑1 (MQ=255)

                                                          agaTGATGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCAAACAGAt                               <  2:1074831/93‑1 (MQ=255)

                                                           gatgatGTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCAAACAGATc                              >  2:2762668/1‑93 (MQ=255)

                                                            atgatgTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCAAACAGATcc                             >  1:3644403/1‑93 (MQ=255)

                                                            atgatgTGATGATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCAAACAGATcc                             <  1:3591772/93‑1 (MQ=255)

                                                              gatgtgatgATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCAAACAGATCCGc                           >  1:116083/1‑93 (MQ=255)

                                                                 gtgatgATTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCAAACAGATCCGCCCg                        >  2:3289105/1‑93 (MQ=255)

                                                                  tgatgaTTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCAAACAGATCCGCCCGt                       >  1:79727/1‑93 (MQ=255)

                                                                     tgatTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCAAACAGATCCGCCCGTTaa                    <  2:746923/93‑1 (MQ=255)

                                                                       atTGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCAAACAGATCCGCCCGTTAAGt                  <  2:337227/93‑1 (MQ=255)

                                                                         tGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCAAACAGATCCGCCCGTTAAGtga                <  1:6178403/93‑1 (MQ=255)

                                                                         tGGTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCAAACAGATCCGCCCGTTAAGtga                <  2:453894/93‑1 (MQ=255)

                                                                          ggTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCAAACAGATCCGCCCGTTAAGtgat               <  1:5955757/93‑1 (MQ=255)

                                                                           gTGCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCAAACAGATCCGCCCGTTAAGtgatg              <  2:5346757/93‑1 (MQ=255)

                                                                             gCGGGTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCAAACAGATCCGCCCGTTAAGTGATGaa            <  2:5675707/93‑1 (MQ=255)

                                                                               gggTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCAAACAGATCCGCCCGTTAAGTGATGAaga          <  2:1709553/93‑1 (MQ=255)

                                                                               gggTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCAAACAGATCCGCCCGTTAAGTGATGAaga          <  2:2824211/93‑1 (MQ=255)

                                                                               gggTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCAAACAGATCCGCCCGTTAAGTGATGAaga          <  1:3091945/93‑1 (MQ=255)

                                                                                 gTAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCAAACAGATCCGCCCGTTAAGTGATGAAgaga        <  1:1425787/93‑1 (MQ=255)

                                                                                  tAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCAAACAGATCCGCCCGTTAAGTGATGAAGAGaa       <  2:3853519/93‑1 (MQ=255)

                                                                                  tAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCAAACAGATCCGCCCGTTAAGTGATGAAGAGaa       >  2:4990889/1‑93 (MQ=255)

                                                                                  tAGTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCAAACAGATCCGCCCGTTAAGTGATGAAGAGaa       <  2:1998600/93‑1 (MQ=255)

                                                                                    gTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCAAACAGATCCGCCCGTTAAGTGATGAAGAGAAgg     >  1:4829751/1‑93 (MQ=255)

                                                                                    gTCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCAAACAGATCCGCCCGTTAAGTGATGAAGAGAAgg     >  1:650501/1‑93 (MQ=255)

                                                                                     tCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCAAACAGATCCGCCCGTTAAGTGATGAAGAGAAGGc    <  2:2847677/93‑1 (MQ=255)

                                                                                     tCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCAAACAGATCCGCCCGTTAAGTGATGAAGAGAAGGc    <  1:1817563/93‑1 (MQ=255)

                                                                                     tCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCAAACAGATCCGCCCGTTAAGTGATGAAGAGAAGGc    <  2:4184795/93‑1 (MQ=255)

                                                                                     tCTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCAAACAGATCCGCCCGTTAAGTGATGAAGAGAAGGc    <  1:1303477/93‑1 (MQ=255)

                                                                                      cTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCAAACAGATCCGCCCGTTAAGTGATGAAGAGAAGGCt   <  2:783244/93‑1 (MQ=255)

                                                                                      cTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCAAACAGATCCGCCCGTTAAGTGATGAAGAGAAGGCt   <  1:4104695/93‑1 (MQ=255)

                                                                                      cTGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCAAACAGATCCGCCCGTTAAGTGATGAAGAGAAGGCt   <  1:209193/93‑1 (MQ=255)

                                                                                       tGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCAAACAGATCCGCCCGTTAAGTGATGAAGAGAAGGCTg  <  2:263320/93‑1 (MQ=255)

                                                                                       tGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCAAACAGATCCGCCCGTTAAGTGATGAAGAGAAGGCTg  >  1:880684/1‑93 (MQ=255)

                                                                                       tGGCCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCAAACAGATCCGCCCGTTAAGTGATGAAGAGAAGGCTg  <  1:2705600/93‑1 (MQ=255)

TTGGCAATCGAGTTTTGGTTGGGATGGGCTCAATTTTACTTGATGGCGCAATAGTAGAAGATGATGTGATGATTGGTGCGGGTAGTCTGGTCCCACAAAATAAACGGCTGGAGAGCGGATATCTGTATCTCGGTAGCCCCGTCAAACAGATCCGCCCGTTAAGTGATGAAGAGAAGGCTG  >  NC_000913/3427171‑3427350

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