BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,150,581	T→C	Y111H (TAT→CAT)	metC →	cystathionine beta‑lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,150,581	0	T	C	84.0%	804.6 / 147.3	356	Y111H (TAT→CAT)	metC	cystathionine beta‑lyase
Reads supporting (aligned to +/- strand): ref base T (30/27); new base C (157/142); total (187/169)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.35e-01

Only 100 of 357 total aligned reads displayed.

GTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTT  >  NC_000913/3150491‑3150673

 tGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATg                                                                                           >  2:5574225/1‑93 (MQ=255)

 tGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATg                                                                                           <  1:4239754/93‑1 (MQ=255)

   cTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGaa                                                                                         >  1:3601103/1‑93 (MQ=255)

    tATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAAc                                                                                        >  1:1555806/1‑93 (MQ=255)

     aTTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAAcc                                                                                       >  2:1950225/1‑93 (MQ=255)

       ttCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGa                                                                                     <  1:6574525/93‑1 (MQ=255)

       ttCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGa                                                                                     <  1:2195897/93‑1 (MQ=255)

        tCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAg                                                                                    <  1:4558495/93‑1 (MQ=255)

          ccTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTc                                                                                  >  2:2545841/1‑93 (MQ=255)

          ccTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTc                                                                                  >  1:3201056/1‑93 (MQ=255)

           cTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCa                                                                                 <  2:114421/93‑1 (MQ=255)

            tGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAg                                                                                <  1:1380154/93‑1 (MQ=255)

             gCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAgg                                                                               <  1:3512926/93‑1 (MQ=255)

              cGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGa                                                                              <  2:1851339/93‑1 (MQ=255)

               ggggCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGAt                                                                             <  1:3701935/93‑1 (MQ=255)

               ggggCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGAt                                                                             >  1:2954487/1‑93 (MQ=255)

                gggCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGAtt                                                                            >  1:4954089/1‑93 (MQ=255)

                  gcggcAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTc                                                                          <  2:6605099/93‑1 (MQ=255)

                    ggcAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTg                                                                        >  1:3378343/1‑93 (MQ=255)

                     gcAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGt                                                                       >  2:142858/1‑93 (MQ=255)

                     gcAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGt                                                                       >  2:6522279/1‑93 (MQ=255)

                     gcAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGt                                                                       >  2:962168/1‑93 (MQ=255)

                     gcAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGt                                                                       >  1:1745635/1‑93 (MQ=255)

                     gcAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGt                                                                       >  1:2598007/1‑93 (MQ=255)

                      cAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTa                                                                      <  2:6256107/93‑1 (MQ=255)

                      cAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTa                                                                      <  1:2691188/93‑1 (MQ=255)

                      cAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTa                                                                      <  2:5532671/93‑1 (MQ=255)

                       aGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAg                                                                     >  1:885636/1‑93 (MQ=255)

                         cGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCa                                                                   >  2:1705765/1‑93 (MQ=255)

                          ggTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCaa                                                                  >  2:7142829/1‑93 (MQ=255)

                           gTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCaaa                                                                 >  1:6677139/1‑93 (MQ=255)

                            ttGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCaaaa                                                                >  2:4438568/1‑93 (MQ=255)

                              gCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATc                                                              >  1:7113248/1‑93 (MQ=255)

                              gCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATc                                                              >  1:4718668/1‑93 (MQ=255)

                              gCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATc                                                              >  2:6793733/1‑93 (MQ=255)

                                  aTTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCa                                                          >  2:4116855/1‑93 (MQ=255)

                                   ttCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAg                                                         <  2:7243868/93‑1 (MQ=255)

                                    tCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGc                                                        >  1:2368199/1‑93 (MQ=255)

                                    tCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGc                                                        >  2:7463027/1‑93 (MQ=255)

                                      cATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCaa                                                      <  1:2557750/93‑1 (MQ=255)

                                       aTTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCaaa                                                     <  2:48252/93‑1 (MQ=255)

                                        ttCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAAc                                                    >  2:5967096/1‑93 (MQ=255)

                                        ttCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAAc                                                    <  1:3106553/93‑1 (MQ=255)

                                        ttCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAAc                                                    <  1:5000944/93‑1 (MQ=255)

                                         tCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACt                                                   >  1:2512071/1‑93 (MQ=255)

                                          cTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTg                                                  <  1:3555923/93‑1 (MQ=255)

                                          cTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTg                                                  <  1:5616436/93‑1 (MQ=255)

                                          cTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTg                                                  <  1:5574225/93‑1 (MQ=255)

                                           ttGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCGCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTgg                                                 <  2:4238851/93‑1 (MQ=255)

                                             gCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGc                                               >  2:1342084/1‑93 (MQ=255)

                                              cTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCg                                              >  1:632520/1‑93 (MQ=255)

                                               ttttATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGt                                             >  1:166465/1‑93 (MQ=255)

                                               ttttATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGt                                             <  1:1438849/93‑1 (MQ=255)

                                               ttttATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGt                                             <  2:3304793/93‑1 (MQ=255)

                                                  tATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAAc                                          <  1:4441461/93‑1 (MQ=255)

                                                    tCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGa                                        >  1:2798728/1‑93 (MQ=255)

                                                    tCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGa                                        >  2:4972602/1‑93 (MQ=255)

                                                         cAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAACCCTCAGCAAACTGGGCGTAACGAcatca                                   <  2:7047217/93‑1 (MQ=255)

                                                          aGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGAcatcat                                  <  2:2382366/93‑1 (MQ=255)

                                                          aGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGAcatcat                                  <  2:3976581/93‑1 (MQ=255)

                                                           gggCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATg                                 >  2:5342317/1‑93 (MQ=255)

                                                           gggCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATg                                 >  2:6901596/1‑93 (MQ=255)

                                                           gggCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATg                                 >  2:4365252/1‑93 (MQ=255)

                                                           gggCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATg                                 <  2:3408589/93‑1 (MQ=255)

                                                           gggCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATg                                 <  1:1878132/93‑1 (MQ=255)

                                                            ggCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATgg                                <  1:5135039/93‑1 (MQ=255)

                                                                aTCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTg                            <  1:6777817/93‑1 (MQ=255)

                                                                aTCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTg                            >  2:6829479/1‑93 (MQ=255)

                                                                 tCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGa                           >  1:7424263/1‑93 (MQ=255)

                                                                  cATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGAt                          >  1:4817848/1‑93 (MQ=255)

                                                                  cATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGAt                          <  1:2415531/93‑1 (MQ=255)

                                                                  cATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGAt                          <  2:5766412/93‑1 (MQ=255)

                                                                  cATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGAt                          <  1:1127432/93‑1 (MQ=255)

                                                                  cATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGAt                          <  2:226550/93‑1 (MQ=255)

                                                                   aTGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATc                         <  2:4530585/93‑1 (MQ=255)

                                                                   aTGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATc                         <  2:4712745/93‑1 (MQ=255)

                                                                   aTGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATc                         >  1:5775132/1‑93 (MQ=255)

                                                                    tgtgTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATcc                        >  1:6063472/1‑93 (MQ=255)

                                                                    tgtgTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATcc                        >  1:1045468/1‑93 (MQ=255)

                                                                    tgtgTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATcc                        >  1:7005466/1‑93 (MQ=255)

                                                                    tgtgTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATcc                        >  1:5433733/1‑93 (MQ=255)

                                                                    tgtgTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATcc                        >  1:3640364/1‑93 (MQ=255)

                                                                      tgtTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGc                      >  1:5846985/1‑93 (MQ=255)

                                                                      tgtTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGTCATCATGGTTTGATCCGc                      <  1:3869356/93‑1 (MQ=255)

                                                                      tgtTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGc                      <  2:1294636/93‑1 (MQ=255)

                                                                      tgtTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGc                      >  1:981858/1‑93 (MQ=255)

                                                                      tgtTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGc                      >  1:6019047/1‑93 (MQ=255)

                                                                            tgaCCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTg                >  2:1054487/1‑93 (MQ=255)

                                                                            tgaCCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTg                <  2:2078070/93‑1 (MQ=255)

                                                                             gaCCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTgg               >  2:3817543/1‑93 (MQ=255)

                                                                                cAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGc            <  1:7287011/93‑1 (MQ=255)

                                                                                cAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGc            >  2:5453423/1‑93 (MQ=255)

                                                                                   cacCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGa         <  2:2890742/93‑1 (MQ=255)

                                                                                    acCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGat        >  1:7517262/1‑93 (MQ=255)

                                                                                     ccGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGata       <  2:640577/93‑1 (MQ=255)

                                                                                      cGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGatat      <  1:400310/93‑1 (MQ=255)

                                                                                      cGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGatat      <  1:1204094/93‑1 (MQ=255)

                                                                                       gCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATc     <  1:6898430/93‑1 (MQ=255)

                                                                                       gCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATc     <  2:360541/93‑1 (MQ=255)

                                                                                          tATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTCGATCCGCTGATTGGTGCCGATATCGtt  <  2:5168652/93‑1 (MQ=255)

GTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTT  >  NC_000913/3150491‑3150673

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