BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,098,156	A→G	T19A (ACA→GCA)	ycdX →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,098,156	0	A	G	99.6%	794.4 / ‑4.4	227	T19A (ACA→GCA)	ycdX	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (129/97); total (129/98)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.91e-01

Only 100 of 233 total aligned reads displayed.

GCTATCCTTGGCGGGAATTGAATACAGGAGAGAGTTATGTATCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCACATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTGCGATCACCGATCATGGCCCGGATATGGAAGATGCGCCG  >  NC_000913/1098066‑1098245

  tATCCTTGGCGGGAATTGAATACAGGAGAGAGTTATGTATCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCAtt                                                                                       >  2:268489/1‑93 (MQ=255)

    tCCTTGGCGGGAATTGAATACAGGAGAGAGTTATGTATCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTaa                                                                                     <  2:6200542/93‑1 (MQ=255)

    tCCTTGGCGGGAATTGAATACAGGAGAGAGTTATGTATCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTaa                                                                                     <  2:5266114/93‑1 (MQ=255)

       ttGGCGGGAATTGAATACAGGAGAGAGTTATGTATCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTg                                                                                  >  2:4785265/1‑93 (MQ=255)

       ttGGCGGGAATTGAATACAGGAGAGAGTTATGTATCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTg                                                                                  >  2:1364637/1‑93 (MQ=255)

         ggCGGGAATTGAATACAGGAGAGAGTTATGTATCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGAt                                                                                >  1:4031566/1‑93 (MQ=255)

          gCGGGAATTGAATACAGGAGAGAGTTAGGTATCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGAtt                                                                               <  1:2468319/93‑1 (MQ=255)

            gggAATTGAATACAGGAGAGAGTTATGTATCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTAc                                                                             >  1:1416689/1‑93 (MQ=255)

             ggAATTGAATACAGGAGAGAGTTATGTATCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACa                                                                            >  2:1135027/1‑93 (MQ=255)

              gAATTGAATACAGGAGAGAGTTATGTATCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACAt                                                                           <  1:5423687/93‑1 (MQ=255)

                  tGAATACAGGAGAGAGTTATGTATCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGcc                                                                       <  2:3265036/93‑1 (MQ=255)

                   gAATACAGGAGAGAGTTATGTATCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGccc                                                                      >  2:2948758/1‑93 (MQ=255)

                    aaTACAGGAGAGAGTTATGTATCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCa                                                                     <  1:5103042/93‑1 (MQ=255)

                    aaTACAGGAGAGAGTTATGTATCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCa                                                                     <  2:2860303/93‑1 (MQ=255)

                     aTACAGGAGAGAGTTATGTATCCCGTCGACCTTCATATGCATACCGTTGGCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAg                                                                    >  1:6230612/1‑93 (MQ=255)

                     aTACAGGAGAGAGTTATGTATCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAg                                                                    >  1:4099688/1‑93 (MQ=255)

                     aTACAGGAGAGAGTTATGTATCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAg                                                                    >  2:3564574/1‑93 (MQ=255)

                      tACAGGAGAGAGTTATGTATCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAgg                                                                   <  1:1694975/93‑1 (MQ=255)

                       aCAGGAGAGAGTTATGTATCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGc                                                                  >  1:208975/1‑93 (MQ=255)

                       aCAGGAGAGAGTTATGTATCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGc                                                                  <  2:1698812/93‑1 (MQ=255)

                       aCAGGAGAGAGTTATGTATCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGc                                                                  <  2:435740/93‑1 (MQ=255)

                       aCAGGAGAGAGTTATGTATCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGc                                                                  >  2:2472681/1‑93 (MQ=255)

                         aGGAGAGAGTTATGTATCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCca                                                                >  1:2127563/1‑93 (MQ=255)

                          ggAGAGAGTTATGTATCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCcaa                                                               >  2:3216561/1‑93 (MQ=255)

                          ggAGAGAGTTATGTATCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCcaa                                                               >  2:2118583/1‑93 (MQ=255)

                          ggAGAGAGTTATGTATCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCcaa                                                               <  2:2209336/93‑1 (MQ=255)

                            agagagTTATGTATCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCcaaac                                                             >  1:718580/1‑93 (MQ=255)

                            agagagTTATGTATCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCcaaac                                                             >  2:5862049/1‑93 (MQ=255)

                             gagagTTATGTATCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCcaaaca                                                            >  2:3796741/1‑93 (MQ=255)

                             gagagTTATGTATCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCcaaaca                                                            >  1:3975271/1‑93 (MQ=255)

                             gagagTTATGTATCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCcaaaca                                                            <  2:5521119/93‑1 (MQ=255)

                             gagagTTATGTATCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCcaaaca                                                            >  1:2441274/1‑93 (MQ=255)

                              agagTTATGTATCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCcaaacaa                                                           <  2:5116320/93‑1 (MQ=255)

                              agagTTATGTATCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCcaaacaa                                                           >  1:2218732/1‑93 (MQ=255)

                               gagTTATGTATCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCcaaacaaa                                                          >  2:398297/1‑93 (MQ=255)

                                agTTATGTATCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACaaaa                                                         <  2:3680395/93‑1 (MQ=255)

                                agTTATGTATCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACaaaa                                                         <  2:1280413/93‑1 (MQ=255)

                                    aggTATCCCGTCGACCTTCATATGCATACCGTTGCCAGCACCCATGCATATAGCGCATTAAGTGATGACATTGCCCAGGCCAAACAAAAGGGc                                                     <  2:4961576/91‑1 (MQ=255)

                                    aTGTATCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGc                                                     >  2:5183325/1‑93 (MQ=255)

                                     tGTATCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCa                                                    <  1:3513909/93‑1 (MQ=255)

                                      gTATCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCAt                                                   >  1:4776214/1‑93 (MQ=255)

                                      gTATCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCAt                                                   <  1:1800114/93‑1 (MQ=255)

                                       tATCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCAtt                                                  >  1:2184257/1‑93 (MQ=255)

                                        aTCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTa                                                 >  1:553527/1‑93 (MQ=255)

                                        aTCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTa                                                 >  2:690176/1‑93 (MQ=255)

                                        aTCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTa                                                 >  1:2083564/1‑93 (MQ=255)

                                          cccGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTaaa                                               >  1:3747753/1‑93 (MQ=255)

                                           ccGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAAc                                              >  2:4902442/1‑93 (MQ=255)

                                           ccGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAAc                                              <  2:266611/93‑1 (MQ=255)

                                            cGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACt                                             >  1:4678320/1‑93 (MQ=255)

                                             gTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACtt                                            <  1:77694/93‑1 (MQ=255)

                                              tCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACttt                                           <  1:1000646/93‑1 (MQ=255)

                                               cGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACtttt                                          >  2:5871734/1‑93 (MQ=255)

                                               cGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACtttt                                          <  2:3374439/93‑1 (MQ=255)

                                                gACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACttttt                                         >  2:655846/1‑93 (MQ=255)

                                                gACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACttttt                                         >  1:2836491/1‑93 (MQ=255)

                                                gACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCACGCCAAACAAAAGGGCATTAAACttttt                                         >  2:5880038/1‑93 (MQ=255)

                                                 aCCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTg                                        <  1:4725806/93‑1 (MQ=255)

                                                 aCCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTg                                        >  1:4496806/1‑93 (MQ=255)

                                                 aCCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTg                                        <  1:3832002/93‑1 (MQ=255)

                                                  ccTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTGc                                       <  2:776056/93‑1 (MQ=255)

                                                  ccTTCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTGc                                       <  1:6020320/93‑1 (MQ=255)

                                                     tCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTGCGAt                                    >  1:2915121/1‑93 (MQ=255)

                                                     tCATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTGCGAt                                    >  1:5702164/1‑93 (MQ=255)

                                                      cATATGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTGCGATc                                   <  1:5044996/93‑1 (MQ=255)

                                                        tatGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTGCGATCAc                                 >  2:5632027/1‑93 (MQ=255)

                                                        tatGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTGCGATCAc                                 >  2:3885977/1‑93 (MQ=255)

                                                         atGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTGCGATCAcc                                <  2:4774567/93‑1 (MQ=255)

                                                          tGCATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTGCGATCACCg                               >  1:1548294/1‑93 (MQ=255)

                                                            cATACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTGCGATCACCGAt                             <  2:2897481/93‑1 (MQ=255)

                                                             aTACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTGCGATCACCGATc                            <  1:5453338/93‑1 (MQ=255)

                                                             aTACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTGCGATCACCGATc                            <  1:494038/93‑1 (MQ=255)

                                                              tACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTGCGATCACCGATCa                           >  1:6138302/1‑93 (MQ=255)

                                                              tACCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTGCGATCACCGATCa                           >  1:3790951/1‑93 (MQ=255)

                                                               aCCGTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTGCGATCACCGATCAt                          <  1:911497/93‑1 (MQ=255)

                                                                  gTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTGCGATCACCGATCATGGc                       >  2:3155086/1‑93 (MQ=255)

                                                                  gTTGCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTGCGATCACCGATCATGGc                       <  1:536063/93‑1 (MQ=255)

                                                                     gCCAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTGCGATCACCGATCATGGCCCg                    >  1:641537/1‑93 (MQ=255)

                                                                      ccAGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTGCGATCACCGATCATGGCCCgg                   >  1:3371352/1‑93 (MQ=255)

                                                                        aGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTGCGATCACCGATCATGGCCCGGat                 >  1:5338605/1‑93 (MQ=255)

                                                                        aGCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTGCGATCACCGATCATGGCCCGGat                 >  1:3159511/1‑93 (MQ=255)

                                                                         gCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTGCGATCACCGATCATGGCCCGGata                >  2:257207/1‑93 (MQ=255)

                                                                         gCACACATGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTGCGATCACCGATCATGGCCCGGata                <  2:4558458/93‑1 (MQ=255)

                                                                              caTGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTGCGATCACCGATCATGGCCCGGATATGGaa           <  2:6092740/93‑1 (MQ=255)

                                                                               aTGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTGCGATCACCGATCATGGCCCGGATATGGAAg          <  2:639299/93‑1 (MQ=255)

                                                                                tGCATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTGCGATCACCGATCATGGCCCGGATATGGAAGa         <  1:4305731/93‑1 (MQ=255)

                                                                                  cATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTGCGATCACCGATCATGGCCCGGATATGGAAGATg       <  2:5312534/93‑1 (MQ=255)

                                                                                  cATATAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTGCGATCACCGATCATGGCCCGGATATGGAAGATg       <  2:2606102/93‑1 (MQ=255)

                                                                                   atatAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTGCGATCACCGATCATGGCCCGGATATGGAAGATgc      >  2:4671365/1‑93 (MQ=255)

                                                                                   atatAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTGCGATCACCGATCATGGCCCGGATATGGAAGATgc      <  2:3478790/93‑1 (MQ=255)

                                                                                   atatAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTGCGATCACCGATCATGGCCCGGATATGGAAGATgc      >  1:2679960/1‑93 (MQ=255)

                                                                                   atatAGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTGCGATCACCGATCATGGCCCGGATATGGAAGATgc      >  1:1366405/1‑93 (MQ=255)

                                                                                    tataGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTGCGATCACCGATCATGGCCCGGATATGGAAGATgcg     <  2:1112885/93‑1 (MQ=255)

                                                                                    tataGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTGCGATCACCGATCATGGCCCGGATATGGAAGATgcg     <  2:2381165/93‑1 (MQ=255)

                                                                                     ataGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTGCGATCACCGATCATGGCCCGGATATGGAAGATgcgc    <  2:1974112/93‑1 (MQ=255)

                                                                                     ataGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTGCGATCACCGATCATGGCCCGGATATGGAAGATgcgc    >  1:4378054/1‑93 (MQ=255)

                                                                                     ataGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTGCGATCACCGATCATGGCCCGGATATGGAAGATgcgc    <  1:3023034/93‑1 (MQ=255)

                                                                                     ataGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTGCGATCACCGATCATGGCCCGGATATGGAAGATgcgc    >  1:2341278/1‑93 (MQ=255)

                                                                                       aGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTGCGATCACCGATCATGGCCCGGATATGGAAGATGcgccg  >  2:5450093/1‑93 (MQ=255)

                                                                                       aGCGCATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTGCGATCACCGATCATGGCCCGGATATGGAAGATGcgccg  >  1:1140745/1‑93 (MQ=255)

GCTATCCTTGGCGGGAATTGAATACAGGAGAGAGTTATGTATCCCGTCGACCTTCATATGCATACCGTTGCCAGCACACATGCATATAGCACATTAAGTGATTACATTGCCCAGGCCAAACAAAAGGGCATTAAACTTTTTGCGATCACCGATCATGGCCCGGATATGGAAGATGCGCCG  >  NC_000913/1098066‑1098245

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