BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,860,862	T→C	A123A (GCA→GCG)	b4553 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,860,862	0	T	C	99.7%	1220.6 / ‑5.5	387	A123A (GCA→GCG)	b4553	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (220/166); total (220/167)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

Only 100 of 399 total aligned reads displayed.

CAGCGCCGTCACAATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCTGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGTTCCAGGTCAGGAAACT  >  NC_000913/3860772‑3860951

cAGCGCCGTCACAATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGc                                                                                         >  1:6168828/1‑93 (MQ=255)

     ccGTCACAATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCAccc                                                                                    >  2:197362/1‑92 (MQ=255)

     ccGTCACAATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACcg                                                                                    <  1:2821920/93‑1 (MQ=255)

     ccGTCACAATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAACCCCGGCCATCATTGTCAGACCGCTACCCGCCACcg                                                                                    >  1:1706072/1‑93 (MQ=255)

       gTCACAATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACcgcg                                                                                  >  2:7924880/1‑93 (MQ=255)

        tCACAATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGt                                                                                 >  2:5986224/1‑93 (MQ=255)

        tCACAATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGt                                                                                 >  1:23212/1‑93 (MQ=255)

          acaATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTcc                                                                               >  2:6336387/1‑93 (MQ=255)

           caATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCt                                                                              <  2:2935460/93‑1 (MQ=255)

            aaTGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTg                                                                             >  1:1909027/1‑93 (MQ=255)

            aaTGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTg                                                                             >  1:5580471/1‑93 (MQ=255)

             aTGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGa                                                                            >  2:4169758/1‑93 (MQ=255)

             aTGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGa                                                                            <  1:7136954/93‑1 (MQ=255)

              tGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAg                                                                           <  2:4561076/93‑1 (MQ=255)

              agCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCAATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAg                                                                           <  2:1620928/92‑1 (MQ=255)

                 cGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGa                                                                        >  1:4149303/1‑93 (MQ=255)

                  ggAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGaa                                                                       >  1:5919100/1‑93 (MQ=255)

                  ggAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGGTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGaa                                                                       >  1:759340/1‑93 (MQ=255)

                   gAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGaaa                                                                      >  2:432560/1‑93 (MQ=255)

                    aaaGGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAAt                                                                     <  2:5842746/93‑1 (MQ=255)

                     aaTGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATc                                                                    >  1:4774555/1‑93 (MQ=255)

                     aaTGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATc                                                                    <  1:3228753/93‑1 (MQ=255)

                      aTGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAAtag                                                                   >  1:2140398/1‑91 (MQ=255)

                      aTGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCg                                                                   <  1:680323/93‑1 (MQ=255)

                       tGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGa                                                                  <  2:7297264/93‑1 (MQ=255)

                        gATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGAc                                                                 <  1:4330047/93‑1 (MQ=255)

                        gATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGAc                                                                 >  1:4993171/1‑93 (MQ=255)

                           aaTTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGcc                                                              <  2:3132227/93‑1 (MQ=255)

                           aaTTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGcc                                                              >  1:2407353/1‑93 (MQ=255)

                            aTTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCg                                                             >  1:7832017/1‑93 (MQ=255)

                              tGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTCCCCGCCACCGCGTCCTGAGTGAAATCGACGCCGaa                                                           <  2:6445668/93‑1 (MQ=255)

                              tGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGaa                                                           <  1:7726963/93‑1 (MQ=255)

                              tGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGaa                                                           <  2:224602/93‑1 (MQ=255)

                               gcgcCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAg                                                          >  2:4710646/1‑93 (MQ=255)

                                  ccGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAg                                                       <  1:249388/93‑1 (MQ=255)

                                  ccGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAg                                                       >  2:762695/1‑93 (MQ=255)

                                   cGATAATGCTGGTATCGAGGGGTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGc                                                      >  2:1141411/1‑93 (MQ=255)

                                    gATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCt                                                     <  2:6502182/93‑1 (MQ=255)

                                    gATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCt                                                     >  1:1786423/1‑93 (MQ=255)

                                    gATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGGAGCt                                                     >  2:6175246/1‑93 (MQ=255)

                                     aTAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCtt                                                    <  2:5319086/93‑1 (MQ=255)

                                     aTAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCtt                                                    >  2:7428975/1‑93 (MQ=255)

                                        aggCTCGTATCGAGGGTTATAATCCCGGCCATCATTGTCAGACCGCTACCCGCGCCCGCGTCGTGAGTGAAACCGACGCCGAAGTAGCTTccc                                                 <  2:5076804/91‑1 (MQ=255)

                                         tGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTcccc                                                <  1:1069002/93‑1 (MQ=255)

                                         tGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTcccc                                                >  1:1493482/1‑93 (MQ=255)

                                            tGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAgg                                             <  2:855046/93‑1 (MQ=255)

                                            tGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAgg                                             >  1:5738033/1‑93 (MQ=255)

                                                aTCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCAt                                         <  1:7393395/93‑1 (MQ=255)

                                                 tCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCAtt                                        >  1:4715483/1‑93 (MQ=255)

                                                  cGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTc                                       <  1:5094005/93‑1 (MQ=255)

                                                   gAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTcc                                      >  2:3335242/1‑93 (MQ=255)

                                                    aGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTccc                                     >  2:3762214/1‑93 (MQ=255)

                                                    aGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTccc                                     >  2:6937942/1‑93 (MQ=255)

                                                    aGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTccc                                     >  1:3740894/1‑93 (MQ=255)

                                                        ttttAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATcg                                 >  1:7272381/1‑93 (MQ=255)

                                                         tttAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATcgc                                >  2:3004431/1‑93 (MQ=255)

                                                           tAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGt                              <  1:3296667/93‑1 (MQ=255)

                                                           tAATCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGt                              >  2:2431800/1‑93 (MQ=255)

                                                           tAATCCCCGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGt                              <  1:808059/93‑1 (MQ=255)

                                                            aaTCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGtt                             <  1:5212250/93‑1 (MQ=255)

                                                             aTCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTtg                            >  1:6704576/1‑93 (MQ=255)

                                                             aTCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTtg                            >  1:1988278/1‑93 (MQ=255)

                                                             aTCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTtg                            <  1:3621119/93‑1 (MQ=255)

                                                             aTCCCGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTtg                            >  1:7888887/1‑93 (MQ=255)

                                                               cccGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTtgat                          >  1:42393/1‑93 (MQ=255)

                                                                ccGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTtgatg                         >  2:5423045/1‑93 (MQ=255)

                                                                 cGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTtgatga                        <  1:2222663/93‑1 (MQ=255)

                                                                 cGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTtgatga                        >  2:1266681/1‑93 (MQ=255)

                                                                 cGGCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTtgatga                        <  1:4479451/93‑1 (MQ=255)

                                                                  ggCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGACATCGACGCCGAAGTAACTTCCCCACGTCATTCCCATCCCGTTGATGaa                       >  1:2649356/1‑93 (MQ=255)

                                                                   gCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGaaa                      >  2:438769/1‑93 (MQ=255)

                                                                   gCCATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGaaa                      <  1:6312902/93‑1 (MQ=255)

                                                                    ccATCATTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAAt                     <  2:308676/93‑1 (MQ=255)

                                                                     catcatTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATa                    <  2:2657652/93‑1 (MQ=255)

                                                                     catcatTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATa                    <  1:2774537/93‑1 (MQ=255)

                                                                     catcatTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATa                    <  2:1780567/93‑1 (MQ=255)

                                                                      atcatTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAg                   >  2:3573569/1‑93 (MQ=255)

                                                                      atcatTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAg                   >  2:5690039/1‑93 (MQ=255)

                                                                       tcatTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGt                  >  2:2832128/1‑93 (MQ=255)

                                                                       tcatTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGt                  >  2:4627727/1‑93 (MQ=255)

                                                                       tcatTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGt                  >  1:5371395/1‑93 (MQ=255)

                                                                        catTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGtt                 <  2:3988871/93‑1 (MQ=255)

                                                                         atTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGTTc                >  1:6906267/1‑93 (MQ=255)

                                                                         atTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGTTc                <  1:2988225/93‑1 (MQ=255)

                                                                         atTGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGTTc                <  1:461295/93‑1 (MQ=255)

                                                                          ttGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGTTcc               >  1:1160742/1‑93 (MQ=255)

                                                                          ttGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGTTcc               >  2:7789126/1‑93 (MQ=255)

                                                                           tGTCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGTTCCa              <  1:6464927/93‑1 (MQ=255)

                                                                             tCAGACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGTTCCAgg            <  1:4040782/93‑1 (MQ=255)

                                                                                gACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGTTCCAGGTCa         >  1:965651/1‑93 (MQ=255)

                                                                                gACCGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGTTCCAGGTCa         <  1:6572070/93‑1 (MQ=255)

                                                                                  ccGCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGTTCCAGGTCAgg       <  2:1440411/93‑1 (MQ=255)

                                                                                    gCTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGTTCCAGGTCAGGaa     >  1:1005353/1‑93 (MQ=255)

                                                                                     cTACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGTTCCAGGTCAGGaaa    <  1:129113/93‑1 (MQ=255)

                                                                                      tACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGTTCCAGGTCAGGAAAc   >  1:7507775/1‑93 (MQ=255)

                                                                                      tACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGTTCCAGGTCAGGAAAc   <  2:70417/93‑1 (MQ=255)

                                                                                      tACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGTTCCAGGTCAGGAAAc   <  1:39181/93‑1 (MQ=255)

                                                                                      tACCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGTTCCAGGTCAGGAAAc   <  2:3172678/93‑1 (MQ=255)

                                                                                       aCCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGTTCCAGGTCCGGAAAcc  <  2:4798773/93‑2 (MQ=255)

                                                                                       aCCCGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGTTCCAGGTCAGGAAACt  <  2:6663034/93‑1 (MQ=255)

CAGCGCCGTCACAATGCCGGAAATGATAATTGCGCCGATAATGCTGGTATCGAGGGTTTTAATCCCGGCCATCATTGTCAGACCGCTACCTGCCACCGCGTCCTGAGTGAAATCGACGCCGAAGTAGCTTCCCCAGGTCATTCCCATCGCGTTGATGAAATAGTTCCAGGTCAGGAAACT  >  NC_000913/3860772‑3860951

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