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	99.6%	899.4 / NA	258	Y111H (TAT→CAT)	metC	cystathionine beta‑lyase
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (140/117); total (140/118)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 264 total aligned reads displayed.

GTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATC  >  NC_000913/3150491‑3150670

gTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCAt                                                                                         >  1:842061/1‑93 (MQ=255)

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

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

    tATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATTGAACAGGGCGATCATGTGTTGATTACCAACACCGACCATGCAc                                                                                     >  1:5420486/1‑93 (MQ=255)

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

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

      tttCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCg                                                                                   <  1:5269757/93‑1 (MQ=255)

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

       ttCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGa                                                                                  >  1:5963925/1‑93 (MQ=255)

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

        tCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAg                                                                                 >  2:4741238/1‑93 (MQ=255)

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

        tCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCCAg                                                                                 >  1:4272696/1‑93 (MQ=255)

         cccTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGt                                                                                >  2:3915817/1‑93 (MQ=255)

          ccTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTc                                                                               <  2:2409034/93‑1 (MQ=255)

          ccTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTc                                                                               <  2:2934880/93‑1 (MQ=255)

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

              cGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGa                                                                           <  2:1113500/93‑1 (MQ=255)

               ggggCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGAt                                                                          <  2:4493125/93‑1 (MQ=255)

               ggggCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGAt                                                                          <  1:3088596/93‑1 (MQ=255)

               gggcCGGCAGCGGTTGCTATTTCCATTCTTGCTTTTATCGAACAGGGCTATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGAt                                                                          <  1:976522/93‑1 (MQ=255)

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

                 ggcggcAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGAttt                                                                        <  2:1563247/93‑1 (MQ=255)

                 ggcggcAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGAttt                                                                        >  1:2770302/1‑93 (MQ=255)

                  gcggcAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTc                                                                       >  1:2873281/1‑93 (MQ=255)

                   cggcAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCt                                                                      >  1:5361187/1‑93 (MQ=255)

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

                     gcAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGt                                                                    <  1:1139531/93‑1 (MQ=255)

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

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

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

                        gCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACGCCGCCCATGAACCGAGTCAGGATTTCTGTAGc                                                                 <  2:3992525/93‑1 (MQ=255)

                        gCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGc                                                                 >  2:2743151/1‑93 (MQ=255)

                        gCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGc                                                                 >  2:3712913/1‑93 (MQ=255)

                        gCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGc                                                                 >  2:5929547/1‑93 (MQ=255)

                         cGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCa                                                                <  2:2478223/93‑1 (MQ=255)

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

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

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

                           gTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCaaa                                                              <  2:776463/93‑1 (MQ=255)

                             tGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAAt                                                            >  1:2257913/1‑93 (MQ=255)

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

                               cTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCGACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATcc                                                          <  2:3442528/93‑1 (MQ=255)

                               cTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATcc                                                          <  1:2163008/93‑1 (MQ=255)

                               cTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATcc                                                          >  2:2062071/1‑93 (MQ=255)

                                tAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCt                                                         >  2:3031938/1‑93 (MQ=255)

                                tAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCt                                                         >  2:957490/1‑93 (MQ=255)

                                 aaTTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTc                                                        >  2:1045070/1‑93 (MQ=255)

                                 aaTTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTc                                                        <  2:60369/93‑1 (MQ=255)

                                 aaTTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTc                                                        <  1:2592504/93‑1 (MQ=255)

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

                                     ccATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCa                                                    <  1:406702/93‑1 (MQ=255)

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

                                        ttCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAAc                                                 >  2:5980851/1‑93 (MQ=255)

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

                                          cTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTg                                               >  1:1025968/1‑93 (MQ=255)

                                           ttGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTgg                                              <  1:3882910/93‑1 (MQ=255)

                                            tGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTggg                                             >  1:4383575/1‑93 (MQ=255)

                                             gCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGc                                            >  2:5351219/1‑93 (MQ=255)

                                             gCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGc                                            <  2:2692684/93‑1 (MQ=255)

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

                                               ttttATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGt                                          >  2:3621633/1‑93 (MQ=255)

                                                tttATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTa                                         <  1:6397861/93‑1 (MQ=255)

                                                tttATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTa                                         <  1:6065622/93‑1 (MQ=255)

                                                  tATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAAc                                       <  2:3347765/93‑1 (MQ=255)

                                                      gAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGAca                                   >  1:2852304/1‑93 (MQ=255)

                                                      gAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGGCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGAca                                   >  1:2666119/1‑93 (MQ=255)

                                                        aCAGGGCGATCATGTGTTGATGACCAATGCCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGAcatc                                 <  2:305156/93‑1 (MQ=255)

                                                        aCAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGAcatc                                 >  1:1832889/1‑93 (MQ=255)

                                                        aCAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGAcatc                                 >  1:1989278/1‑93 (MQ=255)

                                                        aCAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGAcatc                                 >  2:4412229/1‑93 (MQ=255)

                                                         cAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGAcatca                                >  1:5218267/1‑93 (MQ=255)

                                                         cAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGAcatca                                <  2:4914967/93‑1 (MQ=255)

                                                          aGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGAcatcat                               >  2:5461757/1‑93 (MQ=255)

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

                                                            ggCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATgg                             >  1:5835885/1‑93 (MQ=255)

                                                             gCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGt                            >  1:3258749/1‑93 (MQ=255)

                                                              cGATCATGTGTTGATGACCCACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGtt                           >  1:6317790/1‑93 (MQ=255)

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

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

                                                                   aTGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATc                      <  2:2960502/93‑1 (MQ=255)

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

                                                                    tgtgTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATcc                     <  1:2353535/93‑1 (MQ=255)

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

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

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

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

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

                                                                       gtTGATGACCAACTCCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCt                  <  1:1605257/93‑1 (MQ=255)

                                                                       gtTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCt                  <  2:695529/93‑1 (MQ=255)

                                                                          gatgaCCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGAt               >  1:3521199/1‑93 (MQ=255)

                                                                          gatgaCCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGAt               <  1:4656340/93‑1 (MQ=255)

                                                                           atgaCCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGAtt              <  2:5022414/93‑1 (MQ=255)

                                                                              aCCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGt           >  1:4746610/1‑93 (MQ=255)

                                                                               ccAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTg          <  2:2577676/93‑1 (MQ=255)

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

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

                                                                                 aaCACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGcc        >  1:1850590/1‑93 (MQ=255)

                                                                                       gCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGatatg  >  2:3964185/1‑92 (MQ=255)

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

GTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATC  >  NC_000913/3150491‑3150670

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