BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,053,401	T→C	A657A (GCA→GCG)	torS ←	TorS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,053,401	0	T	C	97.1%	407.2 / 0.6	136	A657A (GCA→GCG)	torS	TorS
Reads supporting (aligned to +/- strand): ref base T (2/2); new base C (69/63); total (71/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

Only 100 of 141 total aligned reads displayed.

GGTAATTCGCTGGGTTAGCGGGTTATCTTCAATTAACAGCAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTG  >  NC_000913/1053311‑1053490

ggTAATTCGCTGGGTTAGCGGGTTATCTTCAATTAACAGCAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGgcgc                                                                                         >  1:3594236/1‑93 (MQ=255)

    aTTCGCTGGGTTAGCGGGTTATCTTCAATTAACAGCAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTg                                                                                     >  2:4683322/1‑93 (MQ=255)

    aTTCGCTGGGTTAGCGGGTTATCTTCAATTAACAGCAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTg                                                                                     >  2:2058494/1‑93 (MQ=255)

    aTTCGCTGGGTTAGCGGGTTATCTTCAATTAACAGCAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTg                                                                                     <  2:1258480/93‑1 (MQ=255)

      tCGCTGGGTTAGCGGGTTATCTTCAATTAACAGCAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGc                                                                                   <  2:222391/93‑1 (MQ=255)

         cTGGGTTAGCGGGTTATCTTCAATTAACAGCAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAac                                                                                >  2:543099/1‑93 (MQ=255)

         cTGGGTTAGCGGGTTATCTTCAATTAACAGCAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAac                                                                                <  2:3008491/93‑1 (MQ=255)

         cTGGGTTAGCGGGTTATCTTCAATTAACAGCAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAac                                                                                <  2:3008490/93‑1 (MQ=255)

          tGGGTTAGCGGGTTATCTTCAATTAACAGCAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAaca                                                                               <  2:2850172/93‑1 (MQ=255)

          tGGGTTAGCGGGTTATCTTCAATTAACAGCAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAaca                                                                               <  2:2679147/93‑1 (MQ=255)

            ggTTAGCGGGTTATCTTCAATTAACAGCAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACg                                                                             <  2:666237/93‑1 (MQ=255)

            ggTTAGCGGGTTATCTTCAATTAACAGCAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACg                                                                             <  2:3019567/93‑1 (MQ=255)

             gTTAGCGGGTTATCTTCAATTAACAGCAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGt                                                                            <  1:3069715/93‑1 (MQ=255)

              ttAGCGGGTTATCTTCAATTAACAGCAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTa                                                                           >  2:1581850/1‑93 (MQ=255)

               tAGCGGGTTATCTTCAATTAACAGCAAACGTAAACCGTCAAGGCGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTaa                                                                          <  1:3926250/93‑1 (MQ=255)

               tAGCGGGTTATCTTCAATTAACAGCAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTaa                                                                          <  1:3113532/93‑1 (MQ=255)

                aGCGGGTTATCTTCAATTAACAGCAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAAc                                                                         >  2:2673459/1‑93 (MQ=255)

                   gggTTATCTTCAATTAACAGCAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGc                                                                      >  2:785231/1‑93 (MQ=255)

                   gcgTTATCTTCAATTAACAGCAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGc                                                                      <  2:3320855/91‑1 (MQ=255)

                     gTTATCTTCAATTAACAGCAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCaa                                                                    <  2:1245390/93‑1 (MQ=255)

                     gTTATCTTCAATTAACAGCAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCaa                                                                    >  2:1390442/1‑93 (MQ=255)

                     gTTATCTTCAATTAACAGCAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCaa                                                                    >  2:1756030/1‑93 (MQ=255)

                     gTTATCTTCAATTAACAGCAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCaa                                                                    <  2:2923659/93‑1 (MQ=255)

                     gTTATCTTCAATTAACAGCAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTGGGCACGGGCGCCGTGGCAACACGTAACGGCaa                                                                    <  2:2923660/93‑1 (MQ=255)

                     gTTATCTTCAATTAACAGCAAACATAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGGGGCAACCCGTAACGGCaa                                                                    <  2:2057680/93‑1 (MQ=255)

                       tATCTTCAATTAACAGCAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGc                                                                  <  1:177614/93‑1 (MQ=255)

                        aTCTTCAATTAACAGCAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCg                                                                 >  2:3727056/1‑93 (MQ=255)

                        aTCTTCAATTAACAGCAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCg                                                                 >  2:3335812/1‑93 (MQ=255)

                            tCAATTAACAGCAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCATaaa                                                             <  1:3432203/93‑1 (MQ=255)

                              aaTTAACAGCAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACa                                                           >  1:672101/1‑93 (MQ=255)

                              aaTTAACAGCAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACa                                                           >  2:3125346/1‑93 (MQ=255)

                                 tAACAGCAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACaaaa                                                        <  1:3328223/93‑1 (MQ=255)

                                  aaCAGCAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACaaaaa                                                       >  1:1204429/1‑93 (MQ=255)

                                   aCAGCAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAAc                                                      >  2:1835097/1‑93 (MQ=255)

                                    cAGCAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTTGCACGGGCGCCGGGGCAACACGTAACGGCAAGCGTAAACAAAAACa                                                     >  2:3412135/1‑93 (MQ=255)

                                     aGCAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAg                                                    >  2:1771105/1‑93 (MQ=255)

                                     aGCAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAg                                                    >  2:1459986/1‑93 (MQ=255)

                                     aGCAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAg                                                    >  2:4604438/1‑93 (MQ=255)

                                     aGCAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAg                                                    <  1:1133925/93‑1 (MQ=255)

                                      gCAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGc                                                   <  2:3801006/93‑1 (MQ=255)

                                        aaaCGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAAGAAAAACAGCtt                                                 >  1:3435973/1‑93 (MQ=255)

                                        aaaCGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCtt                                                 >  1:3662980/1‑93 (MQ=255)

                                        aaaCGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCtt                                                 >  2:1851462/1‑93 (MQ=255)

                                         aaCGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTc                                                <  1:1080123/93‑1 (MQ=255)

                                          aCGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTcc                                               <  2:4577299/93‑1 (MQ=255)

                                          aCGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTcc                                               <  2:3446021/93‑1 (MQ=255)

                                          aCGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTcc                                               >  2:3164637/1‑93 (MQ=255)

                                          aCGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTcc                                               <  1:4727764/93‑1 (MQ=255)

                                           cGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCa                                              <  1:4532086/93‑1 (MQ=255)

                                            gTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCac                                             >  1:1487688/1‑92 (MQ=255)

                                            gTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCaa                                             <  1:2717920/93‑1 (MQ=255)

                                            gTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCaa                                             >  2:4242374/1‑93 (MQ=255)

                                            gTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCCCCGTGGCACCACGTAACGGCAAGCGTAAACAAAAACAGCTTCCac                                             >  1:1487687/1‑92 (MQ=255)

                                              aaaCCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAAcc                                           >  2:624214/1‑93 (MQ=255)

                                              aaaCCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAAcc                                           <  1:3622469/93‑1 (MQ=255)

                                              aaaCCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAAcc                                           <  2:3837077/93‑1 (MQ=255)

                                               aaCCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCt                                          >  2:1859599/1‑93 (MQ=255)

                                               aaCCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCt                                          >  1:319772/1‑93 (MQ=255)

                                               aaCCGTCAAGACGCACCGCCTGATTGACTGTTTGTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCt                                          <  1:2122521/93‑1 (MQ=255)

                                                aCCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTc                                         >  1:2751094/1‑93 (MQ=255)

                                                 ccGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTcc                                        <  2:3307113/93‑1 (MQ=255)

                                                  cGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCg                                       >  1:4314889/1‑93 (MQ=255)

                                                  cGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCg                                       <  2:4392970/93‑1 (MQ=255)

                                                  cGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCg                                       >  1:1998198/1‑93 (MQ=255)

                                                  cGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCg                                       >  1:606376/1‑93 (MQ=255)

                                                   gTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCgg                                      >  1:988407/1‑93 (MQ=255)

                                                   gTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCgg                                      <  2:573682/93‑1 (MQ=255)

                                                    tCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGc                                     >  2:2022268/1‑93 (MQ=255)

                                                    tCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGc                                     >  1:1093203/1‑93 (MQ=255)

                                                     cAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCg                                    <  1:2296370/93‑1 (MQ=255)

                                                      aaGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGAGTAAACAAAAACAGCTTCCAACCTCCGGCGt                                   >  1:4282383/1‑93 (MQ=255)

                                                       ataCGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTg                                  >  1:817210/3‑93 (MQ=255)

                                                       aGACGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTg                                  >  1:1156198/1‑93 (MQ=255)

                                                          cGCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCtg                               <  1:439173/93‑1 (MQ=255)

                                                           gCACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCtgg                              <  1:4050354/93‑1 (MQ=255)

                                                            cACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCtggt                             <  1:4147675/93‑1 (MQ=255)

                                                            cACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCtggt                             <  1:3395584/93‑1 (MQ=255)

                                                            cACCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCtggg                             >  1:2516456/1‑92 (MQ=255)

                                                             aCCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCtggtg                            <  2:1572794/93‑1 (MQ=255)

                                                             aCCGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCtggtg                            <  1:4035039/93‑1 (MQ=255)

                                                              ccGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCtggtgg                           >  2:3376430/1‑93 (MQ=255)

                                                              ccGCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCtggtgg                           <  2:3550656/93‑1 (MQ=255)

                                                                gCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGgcg                         >  2:3714735/1‑93 (MQ=255)

                                                                gCCTGATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGgcg                         >  2:3714736/1‑93 (MQ=255)

                                                                    gATTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCa                     >  1:4430847/1‑93 (MQ=255)

                                                                     aTTGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAg                    >  2:608159/1‑93 (MQ=255)

                                                                      ttGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGt                   >  1:4099610/1‑93 (MQ=255)

                                                                      ttGACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGt                   <  1:1970191/93‑1 (MQ=255)

                                                                        gACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTc                 >  2:4192323/1‑93 (MQ=255)

                                                                        gACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTc                 <  1:1878264/93‑1 (MQ=255)

                                                                        gACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTc                 >  1:4384509/1‑93 (MQ=255)

                                                                        gACTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTc                 <  2:2467557/93‑1 (MQ=255)

                                                                         aCTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTcg                >  2:1299971/1‑93 (MQ=255)

                                                                         aCTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTcg                <  1:1673593/93‑1 (MQ=255)

                                                                          cTGTTTTTGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTcgc               >  1:719421/1‑93 (MQ=255)

                                                                                ttGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCa         <  1:2520173/93‑1 (MQ=255)

                                                                                ttGGCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCa         <  1:711801/93‑1 (MQ=255)

                                                                                  ggCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATc       <  1:2653325/93‑1 (MQ=255)

                                                                                  ggCACGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATc       <  2:1073770/93‑1 (MQ=255)

                                                                                      cGGGCGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCt   >  1:2441739/1‑93 (MQ=255)

GGTAATTCGCTGGGTTAGCGGGTTATCTTCAATTAACAGCAAACGTAAACCGTCAAGACGCACCGCCTGATTGACTGTTTTTGGCACGGGTGCCGTGGCAACACGTAACGGCAAGCGTAAACAAAAACAGCTTCCAACCTCCGGCGTGCTGGTGGCGCTCAGTTCGCCGCCCATCGCCTG  >  NC_000913/1053311‑1053490

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