BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,109,226	A→G	L190P (CTA→CCA)	yiiQ ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,109,226	0	A	G	85.3%	719.1 / 115.6	306	L190P (CTA→CCA)	yiiQ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (25/20); new base G (143/118); total (168/138)
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.68e-01

Only 100 of 308 total aligned reads displayed.

GCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAA  >  NC_000913/4109134‑4109318

  aGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGgc                                                                                            <  1:3572045/93‑1 (MQ=255)

    ggCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGgcgc                                                                                          >  1:3890727/1‑93 (MQ=255)

    ggCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAgcgc                                                                                          >  1:6271351/1‑93 (MQ=255)

     gcagcaGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGcc                                                                                         <  1:2606990/93‑1 (MQ=255)

      cagcagGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCa                                                                                        <  2:1910585/93‑1 (MQ=255)

      cagcagGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCa                                                                                        <  1:2304401/93‑1 (MQ=255)

       agcagGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAg                                                                                       >  2:3660140/1‑93 (MQ=255)

        gcagGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGc                                                                                      >  1:751294/1‑93 (MQ=255)

          agGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCtt                                                                                    >  2:3910331/1‑93 (MQ=255)

            gTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTaa                                                                                  <  1:3530564/93‑1 (MQ=255)

             tCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAAt                                                                                 <  1:5692549/93‑1 (MQ=255)

             tCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAAt                                                                                 <  1:2991590/93‑1 (MQ=255)

             tCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAAt                                                                                 <  1:6719113/93‑1 (MQ=255)

              cTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATc                                                                                <  1:2763120/93‑1 (MQ=255)

               tgtgAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCg                                                                               >  1:5472471/1‑93 (MQ=255)

               tgtgAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCg                                                                               >  1:3697571/1‑93 (MQ=255)

                 tgAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGt                                                                             <  2:4453090/93‑1 (MQ=255)

                 tgAAACAGTATAGAGATTCATCAGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGt                                                                             <  2:1077977/93‑1 (MQ=255)

                   aaaCAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTc                                                                           <  2:5055244/93‑1 (MQ=255)

                     aCAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCaa                                                                         <  2:2250108/93‑1 (MQ=255)

                        gTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAg                                                                      >  2:4721017/1‑93 (MQ=255)

                        gTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAg                                                                      >  2:3415296/1‑93 (MQ=255)

                           taGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGa                                                                   <  2:5745676/93‑1 (MQ=255)

                           taGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGa                                                                   >  2:183520/1‑93 (MQ=255)

                             gagaTTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAg                                                                 <  2:3998295/93‑1 (MQ=255)

                             gagaTTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAg                                                                 <  1:2206910/93‑1 (MQ=255)

                             gagaTTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAg                                                                 <  1:3309690/93‑1 (MQ=255)

                             gagaTTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAg                                                                 >  1:7182085/1‑93 (MQ=255)

                             gagaTTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAg                                                                 >  2:1858759/1‑93 (MQ=255)

                             gagaTTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAg                                                                 <  1:5425759/93‑1 (MQ=255)

                              agaTTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAgg                                                                <  1:2681832/93‑1 (MQ=255)

                               gaTTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGt                                                               <  2:7487091/93‑1 (MQ=255)

                                aTTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTc                                                              <  2:6912883/93‑1 (MQ=255)

                                aTTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTc                                                              >  2:2997165/1‑93 (MQ=255)

                                 ttCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCa                                                             <  1:575712/93‑1 (MQ=255)

                                   cATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGc                                                           <  2:2202096/93‑1 (MQ=255)

                                     tCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGccc                                                         <  2:7082106/93‑1 (MQ=255)

                                     tCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGccc                                                         >  2:6140156/1‑93 (MQ=255)

                                         cacaAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCttt                                                     <  1:5186044/93‑1 (MQ=255)

                                           caAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTc                                                   >  2:136031/1‑93 (MQ=255)

                                           caAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTc                                                   >  1:6540995/1‑93 (MQ=255)

                                           caAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTc                                                   >  1:1503739/1‑93 (MQ=255)

                                           caAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTc                                                   >  1:6020050/1‑93 (MQ=255)

                                            aaaGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCg                                                  >  2:6501685/1‑93 (MQ=255)

                                             aaGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGc                                                 >  2:1123302/1‑93 (MQ=255)

                                             aaGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGc                                                 >  2:2518558/1‑93 (MQ=255)

                                             aaGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGc                                                 >  2:679934/1‑93 (MQ=255)

                                              aGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGcc                                                >  2:1790743/1‑93 (MQ=255)

                                              aGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGcc                                                >  1:383037/1‑93 (MQ=255)

                                                gctttgctttTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGt                                              <  1:4527689/93‑1 (MQ=255)

                                                gctttgctttTTGTCATTTATTCAAACCTTCAAGAGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGt                                              <  2:5291422/93‑1 (MQ=255)

                                                 ctttgctttTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGtt                                             >  2:5515004/1‑93 (MQ=255)

                                                  tttgctttTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTg                                            >  2:7165626/1‑93 (MQ=255)

                                                  tttgctttTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTg                                            >  1:3963854/1‑93 (MQ=255)

                                                   ttgctttTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGt                                           <  1:4682536/93‑1 (MQ=255)

                                                    tgctttTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTc                                          <  2:6960984/93‑1 (MQ=255)

                                                     gctttTTGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTcc                                         >  1:1240328/1‑93 (MQ=255)

                                                       tttttGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCcgc                                       <  1:1654999/93‑1 (MQ=255)

                                                          ttGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGac                                    <  1:1073226/93‑1 (MQ=255)

                                                          ttGTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGac                                    >  2:4155639/1‑93 (MQ=255)

                                                            gTCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGacaa                                  >  2:3170993/1‑93 (MQ=255)

                                                             tCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGacaac                                 <  1:6135728/93‑1 (MQ=255)

                                                             tCATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGacaac                                 >  1:2495232/1‑93 (MQ=255)

                                                              cATTTATTCAAACCTTCAAGCGATTCAGATGGCGCGAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGacaaca                                >  2:4521750/1‑93 (MQ=255)

                                                              cATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGacaaca                                >  2:5977588/1‑93 (MQ=255)

                                                              cATTTATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGacaaca                                <  2:2443258/93‑1 (MQ=255)

                                                                tttATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATa                              <  2:511240/93‑1 (MQ=255)

                                                                  tATTCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAAc                            <  1:1182673/93‑1 (MQ=255)

                                                                    ttCAAACCTTCAAGGGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGc                          >  2:6449377/1‑93 (MQ=255)

                                                                    ttCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTTCGCGACAACATAACGc                          <  2:68435/93‑1 (MQ=255)

                                                                     tCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCa                         <  1:943165/93‑1 (MQ=255)

                                                                     tCAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCa                         >  2:2693760/1‑93 (MQ=255)

                                                                     tCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCa                         >  2:5615076/1‑93 (MQ=255)

                                                                      cAAACCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAg                        >  1:4131420/1‑93 (MQ=255)

                                                                       aaaCCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGt                       <  2:1445398/93‑1 (MQ=255)

                                                                       aaaCCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTACGCCGTTGTCCGCGACAACATAACGCAGt                       >  1:1390059/1‑93 (MQ=255)

                                                                       aaaCCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGt                       <  2:352/93‑1 (MQ=255)

                                                                        aaCCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTg                      <  2:2296863/93‑1 (MQ=255)

                                                                        aaCCTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTg                      >  2:5370110/1‑93 (MQ=255)

                                                                          ccTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCa                    <  2:4978585/93‑1 (MQ=255)

                                                                          ccTTCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCa                    >  2:7238026/1‑93 (MQ=255)

                                                                           cTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCAc                   <  2:6249471/93‑1 (MQ=255)

                                                                            ttCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCAcc                  >  1:4889241/1‑93 (MQ=255)

                                                                             tCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCt                 >  2:3773871/1‑93 (MQ=255)

                                                                             tCAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCt                 <  1:4206451/93‑1 (MQ=255)

                                                                              cAAGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCtt                >  2:5023585/1‑93 (MQ=255)

                                                                               aaGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTc               <  1:6864497/93‑1 (MQ=255)

                                                                               aaGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTc               <  1:730043/93‑1 (MQ=255)

                                                                               aaGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTc               <  1:6928607/93‑1 (MQ=255)

                                                                                aGCGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCt              >  1:7412774/1‑93 (MQ=255)

                                                                                 gCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTg             >  2:5188599/1‑93 (MQ=255)

                                                                                  cGATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGt            <  1:5068972/93‑1 (MQ=255)

                                                                                   gATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGtc           <  2:1499443/93‑1 (MQ=255)

                                                                                   gATTCAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGtc           >  1:2371437/1‑93 (MQ=255)

                                                                                     ttCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGtctc         >  1:5993877/1‑93 (MQ=255)

                                                                                      tCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCg        >  1:1745902/1‑93 (MQ=255)

                                                                                       cAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCgg       <  1:7092670/93‑1 (MQ=255)

                                                                                       cAGATGGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCgg       <  2:7362319/93‑1 (MQ=255)

                                                                                          aTAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGctgt    >  1:6158765/1‑93 (MQ=255)

                                                                                           tAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGta   >  2:4321765/1‑93 (MQ=255)

GCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAA  >  NC_000913/4109134‑4109318

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