BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,262,944	G→A	R242R (CGC→CGT)	tdcB ←	threonine dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,262,944	0	G	A	99.4%	634.9 / NA	178	R242R (CGC→CGT)	tdcB	threonine dehydratase
Reads supporting (aligned to +/- strand): ref base G (0/0); major base A (99/78); minor base C (0/1); total (99/79)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 192 total aligned reads displayed.

GGCAATCATACTGTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGGCGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAAAGAAGCCGCCATGCCGTGAAC  >  NC_000913/3262854‑3263033

ggCAATCATACTGTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCcggacg                                                                                         >  2:4018140/1‑93 (MQ=255)

   aaTCATACTGTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGga                                                                                      <  1:4189559/93‑1 (MQ=255)

   aaTCATACTGTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGga                                                                                      <  2:4563782/93‑1 (MQ=255)

   aaTCATACTGTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGga                                                                                      <  1:2927963/93‑1 (MQ=255)

    aTCATACTGTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGgag                                                                                     <  1:944452/93‑1 (MQ=255)

     tCATACTGTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGgaga                                                                                    >  2:4411036/1‑93 (MQ=255)

     tCATACTGTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGgaga                                                                                    >  2:1921401/1‑93 (MQ=255)

     tCATACTGTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGgaga                                                                                    <  2:2833715/93‑1 (MQ=255)

      cATACTGTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGAc                                                                                   <  1:1281016/93‑1 (MQ=255)

       aTACTGTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACa                                                                                  <  2:351957/93‑1 (MQ=255)

       aTACTGTTTCTGAGTTCGTCTTCGCTGACCAGCACGATGTCATCCACTGATTCACGCACGATTTCGTAAGTTAAATTACCCGGACGGGAGACa                                                                                  >  2:2278965/1‑93 (MQ=255)

        tACTGTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACAt                                                                                 >  1:754044/1‑93 (MQ=255)

        tACTGTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACAt                                                                                 <  1:2479311/93‑1 (MQ=255)

         aCTGTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATc                                                                                >  1:3897164/1‑93 (MQ=255)

          cTGTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATca                                                                               >  1:1509570/1‑93 (MQ=255)

            gTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATcaca                                                                             >  1:1068142/1‑93 (MQ=255)

             tttCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACaa                                                                            <  1:3110374/93‑1 (MQ=255)

              ttCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAAc                                                                           >  2:1466024/1‑93 (MQ=255)

               tCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAAcc                                                                          >  2:3425971/1‑93 (MQ=255)

                cTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCa                                                                         >  2:3184421/1‑93 (MQ=255)

                cTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCa                                                                         <  2:4440347/93‑1 (MQ=255)

                cTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCa                                                                         <  2:4624744/93‑1 (MQ=255)

                  gATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATc                                                                       >  2:322642/1‑93 (MQ=255)

                  gATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATc                                                                       >  1:2734979/1‑93 (MQ=255)

                   aTTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATcc                                                                      <  2:2420314/93‑1 (MQ=255)

                   aTTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATcc                                                                      >  1:501415/1‑93 (MQ=255)

                      tCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGcc                                                                   >  2:1634709/1‑93 (MQ=255)

                      tCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACACGGACGGGAGACATCACAACCATCCGcc                                                                   >  2:3344040/1‑93 (MQ=255)

                       cGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCa                                                                  <  2:4432047/93‑1 (MQ=255)

                        gTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAg                                                                 <  1:1667951/93‑1 (MQ=255)

                        gTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAg                                                                 >  1:3046236/1‑93 (MQ=255)

                        gTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAg                                                                 >  2:2790228/1‑93 (MQ=255)

                         tCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAgg                                                                >  2:4644724/1‑93 (MQ=255)

                         tCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAgg                                                                >  2:3150026/1‑93 (MQ=255)

                         tCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAgg                                                                <  2:2317223/93‑1 (MQ=255)

                         tCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAgg                                                                >  1:3669385/1‑93 (MQ=255)

                          cTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAggg                                                               >  1:4325600/1‑93 (MQ=255)

                          cTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAggg                                                               <  2:2378040/93‑1 (MQ=255)

                           ttCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGt                                                              >  2:2014027/1‑93 (MQ=255)

                           ttCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGt                                                              <  2:711594/93‑1 (MQ=255)

                            tCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTg                                                             >  1:187881/1‑93 (MQ=255)

                              gCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGcc                                                           >  2:118696/1‑93 (MQ=255)

                              gCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGAAATCACAACCATCCGCCAGGCTGcc                                                           >  1:1375474/1‑93 (MQ=255)

                                 gcccAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTCCCCGGCCGGGAGACATCACAACCATCCGCCAGGGTGCCGGt                                                        <  2:3265209/91‑1 (MQ=255)

                                 gACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGt                                                        <  1:3443716/93‑1 (MQ=255)

                                  aCCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTa                                                       <  1:1056847/93‑1 (MQ=255)

                                   ccAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAg                                                      <  1:428580/93‑1 (MQ=255)

                                       cACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCg                                                  >  2:1567671/1‑93 (MQ=255)

                                        aCGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCgg                                                 >  2:453734/1‑93 (MQ=255)

                                           aTGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGc                                              >  2:1360414/1‑93 (MQ=255)

                                           aTGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGc                                              >  1:3130961/1‑93 (MQ=255)

                                             gTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCgg                                            >  2:2696213/1‑92 (MQ=255)

                                              tCATCGACTAATTCACGAACGATTTCGTAAGTTTAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTg                                           >  1:3020018/1‑93 (MQ=255)

                                               cATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTgg                                          <  2:776074/93‑1 (MQ=255)

                                                aTCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGt                                         <  1:305604/93‑1 (MQ=255)

                                                 tCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGtt                                        >  2:3838720/1‑93 (MQ=255)

                                                 tCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGtt                                        >  2:742711/1‑93 (MQ=255)

                                                 tCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGtt                                        >  1:329956/1‑93 (MQ=255)

                                                      tAATTCCCGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTtc                                   <  2:962088/93‑1 (MQ=255)

                                                        aTTCACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTtctc                                 <  2:1394630/93‑1 (MQ=255)

                                                           cACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCgg                              <  2:3192276/93‑1 (MQ=255)

                                                           cACGAACGATTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCgg                              >  2:3814500/1‑93 (MQ=255)

                                                            acgaacgaTTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGa                             <  2:4202679/93‑1 (MQ=255)

                                                            acgaacgaTTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGa                             <  2:4293240/93‑1 (MQ=255)

                                                              gaacgaTTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAgg                           >  2:2490903/1‑92 (MQ=255)

                                                              gaacgaTTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGt                           >  2:1643088/1‑93 (MQ=255)

                                                              gaacgaTTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGt                           <  2:1492715/93‑1 (MQ=255)

                                                              gaacgaTTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGt                           >  1:2071397/1‑93 (MQ=255)

                                                               aacgaTTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTg                          >  1:3950768/1‑93 (MQ=255)

                                                                  gaTTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGaa                       >  2:1040397/1‑93 (MQ=255)

                                                                  gaTTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGGGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGaa                       >  1:1570463/1‑93 (MQ=255)

                                                                   aTTTCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGaaa                      <  2:4315492/93‑1 (MQ=255)

                                                                    tttCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGCGGTTATTTCTCCGGAGTGGAaag                     >  1:1095500/1‑93 (MQ=255)

                                                                     ttCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAaaga                    <  2:1777052/93‑1 (MQ=255)

                                                                     ttCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAaaga                    <  1:2343997/93‑1 (MQ=255)

                                                                      tCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAaagaa                   <  2:2174780/93‑1 (MQ=255)

                                                                      tCGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAaagaa                   <  2:3800002/93‑1 (MQ=255)

                                                                       cGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAaagaag                  <  2:3918268/93‑1 (MQ=255)

                                                                       cGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAaagaag                  <  1:3572850/93‑1 (MQ=255)

                                                                       cGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAaagaag                  <  2:1723332/93‑1 (MQ=255)

                                                                       cGTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAaagaag                  <  1:4491136/93‑1 (MQ=255)

                                                                        gTAAGTTAAATTCCTCGGACTGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTGCGGTGCGTGGTTATTTCTCCGGAGTGGAAAGAAgc                 <  2:2080597/93‑1 (MQ=255)

                                                                        gTAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAAAGAAgc                 >  2:930447/1‑93 (MQ=255)

                                                                         tAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAAAGAAgcc                >  2:4641098/1‑93 (MQ=255)

                                                                         tAAGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAAAGAAgcc                <  2:549724/93‑1 (MQ=255)

                                                                           aGTTAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAAAGAAgccgc              >  1:1196427/1‑93 (MQ=255)

                                                                             ttAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAAAGAAGCCGCCa            >  1:288786/1‑93 (MQ=255)

                                                                              tAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAAAGAAGCCGCCAt           >  1:2669106/1‑93 (MQ=255)

                                                                              tAAATTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAAAGAAGCCGCCAt           >  1:479984/1‑93 (MQ=255)

                                                                                aaTTACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAAAGAAGCCGCCATGc         >  1:1744024/1‑93 (MQ=255)

                                                                                  ttACCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAAAGAAGCCGCCATGCCg       >  1:2200172/1‑93 (MQ=255)

                                                                                    aCCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAAAGAAGCCGCCATGCCGTg     >  2:2197493/1‑93 (MQ=255)

                                                                                    aCCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAAAGAAGCCGCCATGCCGTg     >  2:3643511/1‑93 (MQ=255)

                                                                                    aCCCGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAAAGAAGCCGCCATGCCGTg     <  2:1400035/93‑1 (MQ=255)

                                                                                     cccGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAAAGAAGCCGCCATGCCGTGa    >  1:1138663/1‑93 (MQ=255)

                                                                                     cccGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAAAGAAGCCGCCATGCCGTGa    >  1:659120/1‑93 (MQ=255)

                                                                                     cccGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAAAGAAGCCGCCATGCCGTGa    <  1:1839423/93‑1 (MQ=255)

                                                                                      ccGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAAAGAAGTCGCGATGCCGTGac   >  2:4091809/1‑92 (MQ=255)

                                                                                      ccGGACGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAAAGAAGCCGCCATGCCGTGaa   <  2:2422334/93‑1 (MQ=255)

                                                                                       cggacgggAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAAAGAAGCCGCCATGCCGTGAAc  <  1:610591/93‑1 (MQ=255)

GGCAATCATACTGTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGGCGGGAGACATCACAACCATCCGCCAGGGTGCCGGTAGTTCGGTGCGTGGTTATTTCTCCGGAGTGGAAAGAAGCCGCCATGCCGTGAAC  >  NC_000913/3262854‑3263033

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