BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,333,152	A→G	L16P (CTA→CCA)	yjdB ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,333,152	0	A	G	99.8%	1263.4 / NA	403	L16P (CTA→CCA)	yjdB	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (201/200); total (203/200)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 424 total aligned reads displayed.

ACCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGT  >  NC_000913/4333062‑4333241

 ccTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTtggta                                                                                        <  1:633328/93‑1 (MQ=255)

   tGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAgg                                                                                      <  2:6471227/93‑1 (MQ=255)

   tGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAgg                                                                                      >  2:4122093/1‑93 (MQ=255)

    gTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGgc                                                                                     >  1:5330976/1‑93 (MQ=255)

     tttAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGgcg                                                                                    >  2:1697028/1‑93 (MQ=255)

      ttAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGgcgc                                                                                   >  2:6778399/1‑93 (MQ=255)

        aaaaaaGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCtt                                                                                 >  1:5941975/1‑93 (MQ=255)

          aaaaGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCa                                                                               <  1:61646/93‑1 (MQ=255)

          aaaaGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCa                                                                               >  1:7688246/1‑93 (MQ=255)

            aaGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAAc                                                                             >  2:397965/1‑93 (MQ=255)

            aaGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAAc                                                                             <  1:2834447/93‑1 (MQ=255)

            aaGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGACGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAAc                                                                             <  2:3263836/93‑1 (MQ=255)

               gCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATc                                                                          <  2:831762/93‑1 (MQ=255)

                cAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCa                                                                         >  1:2403887/1‑93 (MQ=255)

                cAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCa                                                                         <  1:3547198/93‑1 (MQ=255)

                cAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCa                                                                         <  1:515194/93‑1 (MQ=255)

                 aaTATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCAc                                                                        >  2:3537563/1‑93 (MQ=255)

                 aaTATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCAc                                                                        >  1:7440564/1‑93 (MQ=255)

                  atatTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACg                                                                       >  1:6463045/1‑93 (MQ=255)

                   tatTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACgg                                                                      <  2:5054157/93‑1 (MQ=255)

                   tatTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACgg                                                                      >  1:1287478/1‑93 (MQ=255)

                   tatTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACgg                                                                      >  1:5846641/1‑93 (MQ=255)

                    atTCAGGCAGATAGAGATATAAAAAGCGGTCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGgt                                                                     <  2:4876046/93‑1 (MQ=255)

                    atTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGgt                                                                     <  2:3784857/93‑1 (MQ=255)

                    atTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGgt                                                                     >  2:5743072/1‑93 (MQ=255)

                     ttCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGgtg                                                                    >  2:2680338/1‑93 (MQ=255)

                     ttCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGgtg                                                                    <  1:1933496/93‑1 (MQ=255)

                     ttCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGgtg                                                                    >  2:6388262/1‑93 (MQ=255)

                      tCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGgtgt                                                                   <  2:1115475/93‑1 (MQ=255)

                       cAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGtt                                                                  >  1:5618274/1‑93 (MQ=255)

                         ggCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTc                                                                >  2:1491468/1‑93 (MQ=255)

                         ggCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTc                                                                <  1:6853870/93‑1 (MQ=255)

                          gCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTcc                                                               <  2:5124574/93‑1 (MQ=255)

                          gCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTcc                                                               >  2:1545243/1‑93 (MQ=255)

                             gATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATc                                                            <  1:5786816/93‑1 (MQ=255)

                             gATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATc                                                            >  2:7791050/1‑93 (MQ=255)

                             gATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATc                                                            >  1:6527441/1‑93 (MQ=255)

                             gATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATc                                                            <  2:2737094/93‑1 (MQ=255)

                               tAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGa                                                          >  1:1888712/1‑93 (MQ=255)

                               tAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGa                                                          >  2:2321611/1‑93 (MQ=255)

                               tAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGa                                                          >  2:5927086/1‑93 (MQ=255)

                                agagATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGaa                                                         <  1:4043680/93‑1 (MQ=255)

                                 gagaTATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAAc                                                        >  2:5251963/1‑93 (MQ=255)

                                 gagaTATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAAc                                                        >  2:7534377/1‑93 (MQ=255)

                                 gagaTATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAAc                                                        >  1:3763688/1‑93 (MQ=255)

                                  agaTATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACa                                                       <  1:6792382/93‑1 (MQ=255)

                                     tataAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAg                                                    <  1:7116140/93‑1 (MQ=255)

                                     tataAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAg                                                    <  2:4377076/93‑1 (MQ=255)

                                     tataAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAg                                                    <  1:4342921/93‑1 (MQ=255)

                                     tataAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAg                                                    <  1:501342/93‑1 (MQ=255)

                                      ataAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGt                                                   >  2:6941140/1‑93 (MQ=255)

                                        aaaaaGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTgc                                                 >  2:6192635/1‑93 (MQ=255)

                                         aaaaGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTgcg                                                >  2:4279947/1‑93 (MQ=255)

                                          aaaGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTgcgc                                               <  2:4214335/93‑1 (MQ=255)

                                           aaGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTATGCGCTTCAACATCACGGTGTTTCCATCGATCATAGTGCGCa                                              >  2:7041614/1‑93 (MQ=255)

                                           aaGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCa                                              >  2:1216049/1‑93 (MQ=255)

                                           aaGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCa                                              >  2:1764356/1‑93 (MQ=255)

                                            aggagCCAACAATAGCCAGGCGAGGAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCat                                             <  2:5369423/89‑1 (MQ=255)

                                            aGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCat                                             >  1:989313/1‑93 (MQ=255)

                                             gCGGCAAACAATAGCCAGGCGAGTAAATTCAAAGAGGATCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCata                                            <  2:494001/93‑1 (MQ=255)

                                               ggCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATg                                          >  1:3962485/1‑93 (MQ=255)

                                                 ccAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCt                                        >  2:1947206/1‑93 (MQ=255)

                                                 ccAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCt                                        >  1:3296652/1‑93 (MQ=255)

                                                  caacaaTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTc                                       <  2:6860112/93‑1 (MQ=255)

                                                  caacaaTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTc                                       <  1:6030458/93‑1 (MQ=255)

                                                   aacaaTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCg                                      >  1:5415953/1‑93 (MQ=255)

                                                     caaTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgca                                    >  1:5228386/1‑93 (MQ=255)

                                                       aTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgcaaa                                  <  2:6246448/93‑1 (MQ=255)

                                                       aTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgcaaa                                  <  1:6681182/93‑1 (MQ=255)

                                                       aTAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgcaaa                                  >  1:2526821/1‑93 (MQ=255)

                                                        tAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgcaaag                                 >  2:4468988/1‑93 (MQ=255)

                                                        tAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgcaaag                                 <  2:5600910/93‑1 (MQ=255)

                                                         aGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCgcaaagc                                <  2:5848320/93‑1 (MQ=255)

                                                              ggCGAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAAt                           >  1:6970273/1‑93 (MQ=255)

                                                                 gAGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCACAGCAAAATTaa                        >  2:1342520/1‑93 (MQ=255)

                                                                  aGTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAg                       >  1:3148711/1‑93 (MQ=255)

                                                                   gTAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGc                      >  2:4732396/1‑93 (MQ=255)

                                                                   gTAAATTCAAAGAGGGTCTTTGTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGc                      <  1:5479383/93‑1 (MQ=255)

                                                                    tAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGcc                     >  2:3890348/1‑93 (MQ=255)

                                                                    tAAATTCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGcc                     <  2:3091164/93‑1 (MQ=255)

                                                                        ttCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAAcc                 <  1:3394143/93‑1 (MQ=255)

                                                                         tCAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCt                >  1:110144/1‑93 (MQ=255)

                                                                          cAAAGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCtt               <  1:1005342/93‑1 (MQ=255)

                                                                           aaaGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTa              >  2:744114/1‑93 (MQ=255)

                                                                            aaGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTaa             >  1:4600855/1‑93 (MQ=255)

                                                                            aaGAGGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTaa             >  1:1390491/1‑93 (MQ=255)

                                                                             agagGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAg            >  2:3177733/1‑93 (MQ=255)

                                                                              gagGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGt           <  2:6039384/93‑1 (MQ=255)

                                                                              gagGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGt           >  2:2098895/1‑93 (MQ=255)

                                                                               agGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGtt          >  1:2103837/1‑93 (MQ=255)

                                                                               agGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGtt          >  1:4445386/1‑93 (MQ=255)

                                                                               agGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGtt          >  2:2517297/1‑93 (MQ=255)

                                                                               agGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGtt          >  1:1762548/1‑93 (MQ=255)

                                                                               agGGTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGtt          <  1:1199739/93‑1 (MQ=255)

                                                                                gggTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTc         <  1:7845872/93‑1 (MQ=255)

                                                                                  gTCTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCtt       <  2:4568713/93‑1 (MQ=255)

                                                                                    cTTTTTGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTaa     <  1:2952140/93‑1 (MQ=255)

                                                                                     tttttGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAg    >  1:725271/1‑93 (MQ=255)

                                                                                     tttttGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAg    <  1:4771720/93‑1 (MQ=255)

                                                                                       tttGGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGt  <  1:2564266/93‑1 (MQ=255)

ACCTGTTTAAAAAAGGCAATATTCAGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGT  >  NC_000913/4333062‑4333241

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