BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,757,874	T→C	H225R (CAT→CGT)	b1680 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,757,874	0	T	C	99.1%	365.9 / ‑4.5	117	H225R (CAT→CGT)	b1680	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (73/42); total (74/43)
Fisher's exact test for biased strand distribution p-value = 3.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.50e-01

Only 100 of 117 total aligned reads displayed.

AGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATCCACCAGCACTTTT  >  NC_000913/1757786‑1757963

aGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGccc                                                                                       <  1:2612301/93‑1 (MQ=255)

 gAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCg                                                                                      >  1:2484576/1‑93 (MQ=255)

 gAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCg                                                                                      <  2:2806302/93‑1 (MQ=255)

    ccGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCggcg                                                                                   >  1:3486913/1‑91 (MQ=255)

    ccGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGgag                                                                                   <  2:3379683/93‑1 (MQ=255)

    ccGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGCGGGGCCATACAGTTTACGCCCGgag                                                                                   >  1:469561/1‑93 (MQ=255)

       cccccTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAac                                                                                >  1:552719/1‑93 (MQ=255)

        ccccTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAaca                                                                               >  1:735275/1‑93 (MQ=255)

        ccccTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAaca                                                                               <  1:3367712/93‑1 (MQ=255)

        ccccTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAaca                                                                               >  1:3905020/1‑93 (MQ=255)

         cccTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAacac                                                                              >  2:3268483/1‑93 (MQ=255)

         cccTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAacac                                                                              <  1:4221538/93‑1 (MQ=255)

         cccTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAacac                                                                              >  1:3470311/1‑93 (MQ=255)

          ccTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACg                                                                             >  1:1816119/1‑93 (MQ=255)

          ccTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACg                                                                             >  1:4085411/1‑93 (MQ=255)

          ccTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACg                                                                             >  1:118375/1‑93 (MQ=255)

            ttCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTa                                                                           >  2:2354764/1‑93 (MQ=255)

            ctCCCACGGCGGCATCTCCTGCAACAAGGCTTCTCTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTa                                                                           <  1:1307678/92‑1 (MQ=255)

              cccACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTaaa                                                                         >  2:3057180/1‑93 (MQ=255)

               ccACGGCGTCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTaaaa                                                                        <  1:4675506/93‑1 (MQ=255)

               ccACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTaaaa                                                                        >  1:542657/1‑93 (MQ=255)

               ccACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTaaaa                                                                        <  2:4220423/93‑1 (MQ=255)

                cACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAg                                                                       <  1:1678351/93‑1 (MQ=255)

                 aCGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGt                                                                      >  1:813094/1‑93 (MQ=255)

                 aCGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGt                                                                      >  2:904478/1‑93 (MQ=255)

                 aCGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGt                                                                      >  2:180256/1‑93 (MQ=255)

                  cggcggCATCTCCTGCAACAAGGCTTCTTTCACATGAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTc                                                                     >  2:4095789/1‑93 (MQ=255)

                    gcggcATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGc                                                                   <  2:493852/93‑1 (MQ=255)

                     cggcATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCa                                                                  <  1:1223986/93‑1 (MQ=255)

                     cggcATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCa                                                                  >  2:1812779/1‑93 (MQ=255)

                      ggcATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCaa                                                                 >  2:2028950/1‑93 (MQ=255)

                         agctcCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCGGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATcc                                                              <  2:4091020/91‑1 (MQ=255)

                         aTCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATcc                                                              >  1:1053407/1‑93 (MQ=255)

                         aTCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATcc                                                              <  2:806679/93‑1 (MQ=255)

                         aTCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATcc                                                              >  2:2151475/1‑93 (MQ=255)

                         aTCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATcc                                                              <  2:1542816/93‑1 (MQ=255)

                           ctcCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAg                                                            <  1:1220884/93‑1 (MQ=255)

                            tcCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAgc                                                           >  2:1806090/1‑93 (MQ=255)

                            tcCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAgc                                                           <  2:1217023/93‑1 (MQ=255)

                            tcCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAgc                                                           <  2:4538317/93‑1 (MQ=255)

                              cTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAgcgc                                                         >  1:1556386/1‑93 (MQ=255)

                              cTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAgcgc                                                         >  2:1066078/1‑93 (MQ=255)

                               tGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGcc                                                        >  1:2473943/1‑93 (MQ=255)

                                 caacaaGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTg                                                      >  2:4381735/1‑93 (MQ=255)

                                  aacaaGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGa                                                     >  1:3261880/1‑93 (MQ=255)

                                  aacaaGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGa                                                     >  2:1838269/1‑93 (MQ=255)

                                  aacaaGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGa                                                     <  1:16372/93‑1 (MQ=255)

                                   acaaGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGaa                                                    <  1:485846/93‑1 (MQ=255)

                                    caaGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAAc                                                   >  2:2500372/1‑93 (MQ=255)

                                    caaGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAAc                                                   >  1:37473/1‑93 (MQ=255)

                                     aaGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACa                                                  <  2:4182751/93‑1 (MQ=255)

                                        gCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATcc                                               <  1:153658/93‑1 (MQ=255)

                                        gCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATcc                                               >  1:807614/1‑93 (MQ=255)

                                        gCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATcc                                               >  2:3800311/1‑93 (MQ=255)

                                         cttcttTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACCTCCa                                              >  2:4509093/1‑93 (MQ=255)

                                         cttcttTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATccc                                              >  1:909952/1‑92 (MQ=255)

                                          ttcttTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCAc                                             >  1:3632109/1‑93 (MQ=255)

                                          ttcttTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCAc                                             >  1:3107839/1‑93 (MQ=255)

                                          ttcttTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCAc                                             >  1:3753750/1‑93 (MQ=255)

                                           tcttTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCAcc                                            <  2:673654/93‑1 (MQ=255)

                                           tcttTAACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCAcc                                            <  2:2596503/93‑1 (MQ=255)

                                             tttCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCgg                                          >  1:790113/1‑93 (MQ=255)

                                              ttCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGga                                         <  2:1788184/93‑1 (MQ=255)

                                                cacaTAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGgatg                                       >  2:2219234/1‑93 (MQ=255)

                                                 acaTAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGgatga                                      <  2:2864649/93‑1 (MQ=255)

                                                  caTAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGgatgat                                     <  2:4045180/93‑1 (MQ=255)

                                                  caTAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGgatgat                                     <  2:4634502/93‑1 (MQ=255)

                                                   aTAAAGAATGCCAATTCCGTTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGatgatg                                    >  2:4568641/1‑93 (MQ=255)

                                                   aTAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGatgatg                                    >  2:1030888/1‑93 (MQ=255)

                                                    tAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGc                                   <  1:1926269/93‑1 (MQ=255)

                                                      aaGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCAt                                 >  2:3496547/1‑93 (MQ=255)

                                                       aGAATGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATc                                <  1:4209217/93‑1 (MQ=255)

                                                          aTGCCAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCAcc                             >  2:1198934/1‑93 (MQ=255)

                                                            tccAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGc                           <  2:2599388/92‑1 (MQ=255)

                                                             ccAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGcc                          <  1:3760441/93‑1 (MQ=255)

                                                             ccAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGcc                          >  1:2949701/1‑93 (MQ=255)

                                                             ccAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGcc                          >  2:3646602/1‑93 (MQ=255)

                                                              cAATTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCGCCGGATGATGCAACGCCGCCt                         <  2:2731619/93‑1 (MQ=255)

                                                               aaTTCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTg                        >  1:4629364/1‑93 (MQ=255)

                                                                aTTCCGGTGGGGCCATACAGCTTAAGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCAGAACATCCACCGGATGATGCATCACCGCCtca                       >  2:741037/1‑91 (MQ=255)

                                                                 ttCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAg                      <  1:2025058/93‑1 (MQ=255)

                                                                  tCCGGTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAgc                     <  2:189248/93‑1 (MQ=255)

                                                                     ggTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGcc                  >  1:2835459/1‑93 (MQ=255)

                                                                      gTGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCa                 >  2:478147/1‑93 (MQ=255)

                                                                       tGGGGCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCAt                >  1:4605300/1‑93 (MQ=255)

                                                                         gggCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATcc              >  1:151922/1‑93 (MQ=255)

                                                                         gggCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATcc              >  1:2465795/1‑93 (MQ=255)

                                                                          ggCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATcca             >  1:3166474/1‑93 (MQ=255)

                                                                          ggCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATcca             >  1:2763855/1‑93 (MQ=255)

                                                                           gCCATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATccac            >  1:917046/1‑93 (MQ=255)

                                                                            ccATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATccacc           >  1:3293013/1‑93 (MQ=255)

                                                                            ccATACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATccacc           <  1:80559/93‑1 (MQ=255)

                                                                               tACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATCCACCAGc        >  1:2317627/1‑93 (MQ=255)

                                                                               tACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATCCACCAGc        >  1:4497346/1‑93 (MQ=255)

                                                                               tACAGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATCCACCAGc        >  1:1438474/1‑93 (MQ=255)

                                                                                aCAGTTTACGCCCGGAGAACTCGTAAAAGTCGCAATCCAGTGCCTTATCATCCACTTGATGATGTATCTCCGCCTGAGCGCCATCCACCAGCa       >  1:4678896/1‑93 (MQ=255)

                                                                                  aGTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATCCACCAGCACt     >  2:2139728/1‑93 (MQ=255)

                                                                                   gTTTACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATCCACCAGCACtt    >  2:623651/1‑93 (MQ=255)

                                                                                    tttACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATCCACCAGCACttt   >  1:3254164/1‑93 (MQ=255)

                                                                                     ttACGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATCCACCAGCACtttt  <  1:4231991/93‑1 (MQ=255)

AGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATCCACCAGCACTTTT  >  NC_000913/1757786‑1757963

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